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1 .  INTRODUCTION 


1.1  Program  Objectives 

The  goal  of  the  work  reported  herein  is  the  real  time  operation  of 
the  cell  detection  and  tracking  algorithms  previously  developed  by 
Environmental  Research  §  Technology,  Inc.  (ERT)  for  the  Air  Force  Geo¬ 
physics  Laboratory  (AFGL) .  Specifically,  the  work  included  the  follow¬ 
ing  tasks:  (1)  encoding  the  automatic  echo  assessment  and  forecasting 
Algorithms  developed  under  Contract  No.  F19628-77-C-0058  on  the  Echo 
Track  and  Significance  Estimator  (ETSE)  Interdata  Model  7-32  Computer; 

(2)  testing  the  algorithms  for  real  time  operation;  (3)  operating  the 
computer  and  writing  assembly  level  programs  on  the  Interdata  7-32  com¬ 
puter  during  the  Joint  Agency  Doppler  Technology  tests  at  the  National 
Severe  Storms  Laboratory,  Norman,  Oklahoma,  March  through  June  1978;  and 
(4)  refining  the  automatic  assessment  and  forecasting  algorithms  based 
upon  a  critical  radar  meteorological  analysis  of  the  use  of  the  algorithms. 

The  cell  detection'*  and  tracking  algorithms  were  developed  to  auto¬ 
matically  process  weather  radar  data  to  provide  real  time  identification 
of  severe  weather  and  short  range  (0-20  minute)  forecast  of  regions  of 
potential  hazard  to  aircraft  operation. 

1.2  Summary  of  Results 


The  cell  detection  and  tracking  algorithms  previously  developed 
under  Air  Force  Contract  (Crane,  1978;  Crane,  1977)  were  designed  for 
real  time  operation  on  a  medium  scale  computer  based  on  our  experience 
with  an  extremely  limited  data  sample  (five  consecutive  azimuth  scans) 
and  with  the  CDC  6600  computation  facility  at  AFGL.  The  subsequent 
development  of  a  viable  real  time  weather  radar  data  processing  system 
for  use  by  the  Air  Force,  however,  required  access  to  a  significantly 
larger  set  of  weather  radar  data. 

The  initial  development  of  the  cell  detection  technique  had  been 
undertaken  for  the  Federal  Aviation  Administration  (FAA)  (Crane,  1976) 
for  use  in  automatic  air  traffic  hazard  detection.  Their  continued 
interest  in  the  processing  scheme  and  its  evaluation  resulted  in  a  con¬ 
tract  between  ERT  and  the  FAA  to  process  several  hours  of  radar  observa¬ 
tions  made  simultaneously  with  aircraft  penetrations.  ERT  was  also  under 
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contract  with  the  Bureau  of  Reclamation,  Department  of  Interior  (BuRec) 
during  the  period  of  performance  of  this  contract  for  the  acquisition 
and  analysis  of  significant  amounts  of  radar  data  for  the  development 
of  precipitation  augmentation  strategies  for  the  High  Plains  (HIPLEX) . 

As  a  part  of  that  contract,  the  cell  detection  and  tracking  algorithms 
were  installed  on  the  CYBER-74  Computer  System  in  Denver  and  used  to 
obtain  statistical  data  on  the  spatial  organization  of  precipitation 
production  within  high  plains  storms. 

The  three  concurrent  programs,  the  development  of  real  time  tech¬ 
niques  for  the  Air  Force,  hazard  detection  algorithm  evaluation  for  the 
FAA,  and  the  spatial  organization  analysis  for  BuRec  provided  the  exper¬ 
ience  with  a  significantly  larger  data  base  needed  for  the  refinement  of 
the  tracking  algorithms  (Task  4)  and  the  development  of  a  viable  real 
time  processing  program. 

The  work  under  this  contract  was  organized  in  the  four  tasks  listed 
in  Section  1.1:  (1)  encoding  the  algorithms  for  real  time  operation  on 

the  Interdata  7-32  computer;  (2)  testing  the  algorithms  for  real  time 
operation;  (3)  operating  the  ETSE  in  Norman,  Oklahoma  during  the  Joint 
Agency  Doppler  Technology  tests;  and  (4)  refining  the  algorithms  for 
improved  tracking  and  short  range  forecasting.  ERT  provided  a  programmer 
operator  for  the  1978  Joint  Agency  Doppler  Technology  tests  in  fulfill¬ 
ment  of  Task  3.  The  detection  and  tracking  algorithms  refined  as  a 
result  of  Task  4  have  been  coded  and  installed  on  the  Interdata  7-32  com¬ 
puter  in  fulfillment  of  Task  1.  Their  description  is  the  subject  of 
this  report.  Program  listings  and  a  copy  of  the  operating  instructions 
are  included  in  Appendix  C  and.D.  The  algorithms  have  been  tested  in  com¬ 
pliance  with  Task  2  and,  as  coded  and  operated,  performed  in  real  time. 

The  algorithms  operate  in  real-time  on  the  Interdata  7-32  computer  as 
required.  Real-time  operation  on  the  Interdata  7-32  in  the  manner  used 
in  non-real-time  analysis,  however,  requires  additional  programs  to  fetch 
and  store  the  raw  radar  data.  For  real-time  operation  at  the  level  of 
performance  of  the  previously  developed  non-real -time  program,  a  new 
operating  system  is  required  for  the  Interdata  7-32  computer  which  util¬ 
izes  the  real-time  interrupt  capability  of  the  machine  to  run  the  cell 
detection  and  tracking  programs  as  background  programs  with  the  data 
averaging  and  storage  programs  in  foreground.  Such  an  undertaking  was 
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beyond  the  scope  of  this  contract,  and  a  simpler  program  has  been  speci¬ 
fied  which  allows  operation  in  real-time  on  alternate  scans;  one  to  fetch 
and  store  the  data,  the  second  to  process  the  data.  This  program,  a 
modification  to  the  TSE  program  provided  by  Raytheon  (Boak  et  al,  1977), 
is  listed  in  Appendix  C. 

1.3  Software  Development 

The  previously  developed  cell  detection  and  tracking  programs  were 
extensively  modified  for  use  on  the  Interdata  7-32  computer.  The  track¬ 
ing  program  described  in  the  final  report  of  the  previous  contract 
(Crane,  1978)  was  completely  rewritten  starting  with  a  new  set  of  algor¬ 
ithms.  Experience  with  the  larger  volume  of  radar  data  available  from 
BuRec  forced  the  program  revision.  Two  major  problems  existed  with  the 
initial  tracking  program;  excessive  computer  storage  requirements  for  the 
large  numbers  of  cells  encountered  in  practice  and  an  inherent  inability 
of  the  program  logic  to  separately  establish  individual  smoothed  track 
velocities  for  each  cell.  The  new  tracking  program  develops  the  volume 
cell  attributes  discussed  in  the  previous  contract  report,  establishes 
the  existence  of  cell  clusters,  provides  an  estimate  of  cell  significance, 
and  maintains  both  instantaneous  and  smoothed  velocity  estimates  for  each 
cell. 

The  cell  detection  program  encoded  for  use  on  the  Interdata  7-32 
computer  is  a  streamlined  version  of  the  original  cell  detection  program. 
The  fixed  contour  outlines  are  developed  not  as  line  segments  enclosing 
the  contoured  area  but  as  azimuth  strobes  within  the  echo  region  in  con¬ 
formity  with  the  ETSE  display  scheme.  Attributes  are  not  generated  for 
the  fixed  contours.  The  cell  detection  subroutine  operates  as  before. 
Addressing  in  the  arrays  used  in  the  subroutine  has  been  extensively 
altered  to  increase  operating  speed. 

The  programs  can  operate  over  a  wide  range  of  reflectivity  thresholds 
but  should  be  used  at  a  relatively  high  reflectivity  threshold,  process¬ 
ing  only  data  with  reflectivity  values  above,  say,  40  dBZ.  This  thres¬ 
hold  was  selected  to  reduce  the  number  of  cells  being  processed  in  the 
computer.  The  reduction  in  the  number  of  cells  being  processed  improves 
operating  speed  and  matches  the  output  requirements.  Experience  with 
storms  in  the  high  plains  indicates  that  severe  storms  produce  large 
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numbers  of  celts.  fell  counts  for  a  JO  dRZ  ref  loot  ivity  threshold  are 
over  ISO  during  the  active  period  of  a  storm;  the  total  number  of  separ¬ 
ate  cells  observed  during  a  storm  often  exceeds  several  thousand.  By 
way  of  contrast,  the  output  requirements  established  by  the  remote  dis¬ 
play  system  (communicated  by  MTU  personnel'),  are  for  not  more  than  1 J 
cells  at  any  one  time.  The  more  than  an  order  of  magnitude  reduct  ion  m 
active  cell  count  can  best  be  accomplished  by  increasing  the  threshold 
reflectivity  value. 

In  addition  to  reflectivity  threshold  selection,  significance  is 
established  using  the  integrated  tangential  shear  of  the  radial  velocity 
for  each  cell.  Tangential  shear  cell  detection  as  previous Iv  coded 
ii'rane,  19'")  is  not  attempted  in  the  streamlined  version  of  the  program. 
1'hc  tangential  shear  data  are  used  to  develop  a  shear  attribute  for  each 
cell  detected  on  the  basis  of  reflectivity  alone,  l'he  program  operates 
1"  detecting  all  cells  that  occur  above  the  processing  threshold  but 
saving  for  tracking  and  output  only  the  lt>  cells  having  the  highest 
reflectivity,  integrated  tangential  shear  product .  Internally,  the  pro¬ 
gram  processes  M  volume  cells  but  onlv  the  1 J  volume  cells  with  the 
highest  reflectivity,  shear  product  are  output  after  each  volume  scan. 

The  program  can  be  modified  to  process  and  output  data  for  more  cells  by 
changing  array  si:es  and  test  limits. 

lie  object  of  the  work  reported  in  this  contract  was  to  streamline 
the  original  version  of  the  cell  detection  and  tracking  programs  for  real¬ 
time  use  on  the  Interdata  ~-.'J  computer  with  the  operating  system  provided 
by  Interdata.  Many  of  the  features  of  the  original  program,  such  as  the 
generation  of  fixed  contour  attributes  and  the  independent  detection  of 
tangential  shear  peaks,  were  removed  to  establish  the  real-time  program. 
1'heso  features  may  be  recovered  onlv  if  the  original  version  of  the  pro¬ 
gram,  the  program  opevatmg  on  the  Mid  ilV-oeOd  computer,  is  installed 
on  the  Interdata  "-3J  for  noil-real-time  processing. 

1.4  Organi rat  ion  of  the  Report 

This  report  considers  onlv  the  software  developed  for  use  on  the 
1 TSV,  documentation  for  Tasks  l,  and  4.  Task  o  covered  the  programmer 
computer  operator  for  the  1 measurements  in  Norman .  Oklahoma.  The 
results  of  that  task  were  reported  in  the  quarterly  reports  and  will  not 
be  considered  further. 
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Background  material  and  algorithm  refinement  hast'd  upon  results 
from  the  FAA  and  RuRec  programs  are  considered  in  Section  Section  o 
documents  the  program  for  the  ITS!  ,  Section  I  suinmari  res  program  status 
and  makes  recommendations  for  future  work.  Program  listings,  flow  charts, 
variable  definition,  and  operating  instructions  will  he  found  in  the 
appendices. 
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2.  BACKGROUND 


2.1  Overview  of  Automatic  Processing  Scheme 

Conventional  weather  radars  produce  large  amounts  of  data  -  a  signi¬ 
ficant  fraction  of  which  is  highly  redundant.  Doppler  radars  produce 
even  larger  amounts  of  data.  Significant  weather  events  may  be  imbedded 
in  the  mass  of  redundant  data.  It  is  the  goal  of  the  automatic  process¬ 
ing  scheme  to  extract  the  relevant  information  from  the  mass  of  data  to 
(1)  reduce  the  data  transmission  requirements  for  the  communication  of 
weather  data  obtained  from  a  radar,  (2)  to  screen  the  data  prior  to 
display  to  meteorologists,  (3)  to  preprocess  the  data  for  automatic 
hazard  detection,  and  (4)  to  prepare  the  data  for  use  in  objective  short 
range  forecasting. 

The  processing  scheme  is  structured  to  use  the  cell  detection 
algorithm  in  on-radar- site  computers  to  perform  the  bulk  of  the  data 
reduction.  The  cell  data  are  then  communicated  to  regional  computers 
(or  to  a  second  program  in  a  stand  alone  radar  data  processor)  for  track¬ 
ing  and  interpretation.  For  a  national  network  of  weather  radars,  the 
tracking  program  would  accomplish  the  task  of  netting  different  radars 
and  developing  a  single,  best  estimate  description  of  the  current  weather 
for  use  in  displays,  hazard  detection  and  warning,  and  short  range  fore¬ 
cast. 

For  this  contract  the  cell  detection  and  tracking  algorithms  are 
operated  in  a  single  computer;  the  final  output  is  track  data  for  the 
12  most  significant  cells.  The  track  data  include  smoothed  cell  veloci¬ 
ties  which  are  used  in  the  tracking  program  for  data  association  and  may 
be  used  externally  for  short  range  forecast. 

The  output  cell  and  track  data  are  for  significant  features  in  the 
larger  mass  of  radar  observations.  Significance  is  defined  in  an  ad  hoc 
manner  using  cell  intensity,  area,  vertical  development,  and  tangential 
shear  data.  Parameters  that  are  intuitively  associated  with  significant 
events  such  as  severe  hail,  severe  thunderstorm  turbulence,  and  tornadoes 
have  been  elected  for  the  determination  of  a  significant  cell.  Opera¬ 
tional  experience  with  the  processing  algorithms  and  a  large  sample  of 
data  is  required  before  the  values  of  the  thresholds  used  to  establish 
cell  significance  can  be  refined.  The  current  algorithm  has  been 
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partially  tested  using  aircraft  penetration  data.  For  the  measurements 
currently  available,  a  positive  correlation  has  been  obtained  between  the 
location  of  significant  cells  and  aircraft  turbulence.  An  example  of 
this  association  is  presented  in  Figures  1  and  2  (output  obtained  from 
the  work  at  ERT  sponsored  by  the  FAA) .  The  significant  cells  are  indi¬ 
cated  by  the  tightly  clustered  symbols  for  reflectivity  values  greater 
than  40  dBZ.  The  time  marks  are  6  km  apart  along  the  aircraft  track. 

The  aircraft  was  within  a  typical  cell  diameter  of  two  significant  cells 
between  1640  and  1641,  a  time  marked  by  the  strongest  acceleration 
fluctuations  (turbulence)  during  the  penetration.  By  way  of  contrast, 
the  remainder  of  the  penetration  was  quiet  and  did  not  show  strong 
acceleration  fluctuations  and  was  not  in  the  vicinity  of  significant 
cells . 

The  display  in  Figure  1  immediately  identifies  the  locations  of  the 
significant  cells  and  graphically  presents  the  essential  data  contained 
in  the  radar  observations.  In  contrast,  a  section  of  a  conventional  con¬ 
tour  display  and  of  the  cell  display  are  presented  in  Figure  3.  The 
essence  of  the  data  is  immediately  evident  in  the  cell  display.  The 
important  45  dBZ  cell  that  results  from  the  strong  updraft  depict  'd 
between  1640  and  1641  MDT  in  the  aircraft  data  is  observed  in  the  cell 
display  but  not  in  the  contour  display.  This  cell  produced  the  strongest 
turbulence.  Even  with  a  color  display,  this  important  region  would  not 
be  evident  although  the  higher  reflectivity  contours  in  Figure  3a  would 
stand  out  more  vividly,  in  the  manner  of  the  darkened  symbols  in  Figure 
3b. 

A  display  such  as  Figure  3b  is  readily  interpreted  but  is  not  the 
normal  weather  radar  display  familiar  to  trained  radar  meteorologists. 

The  cell  display  provides  the  important  details  but  at  a  scale  that  is 
smaller  than  used  by  most  meteorologists  and  since  many  cells  may  be 
observed  at  any  single  time,  a  display  of  all  the  cells  may  prove  to  be 
confusing.  The  number  of  cells  active  at  any  one  time  with  reflectivi¬ 
ties  above  15  dBZ  and  the  number  of  significant  cells  for  a  set  of  obser¬ 
vations  of  showers  in  Kansas  (output  from  data  processed  by  ERT  for 
BuRec)  are  depicted  in  Figure  4  as  a  function  of  time  and  volume  scan 
sequence.  The  total  number  of  cells  present  within  a  25  to  150  km 
annulus  of  Goodland,  Kansas  approached  200  during  the  most  active  part 
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Figure  1  Aircraft  Penetration  22  .July  1976  as  Observed 
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of  the  storm.  The  number  of  significant  cells  was  as  high  as  3b.  A 
display  of  alt  the  cells  would  ho  very  difficult  to  interpret.  A  dis¬ 
play  of  only  the  significant  cells,  though  readily  interpreted  by  a 
computer,  may  still  he  difficult  for  a  meteorologist  to  interpet.  A 
further  reduction  in  the  number  of  significant  cells  is  required.  The 
program  developed  under  this  contract  utilizes  tangential  shear  data  in 
addition  to  reflectivity  information  to  further  reduce  the  number  of 
cells  for  display  as  the  most  significant  cells.  The  utility  of  this 
algorithm  for  significant  cell  selection  still  needs  verification. 

The  small  cells  detected  and  tracked  hy  the  algorithms  developed 
under  this  and  prior  contracts  are  well  behaved  in  time  and  space.  The 
cells  show  vertical  development,  persist,  and  have  average  velocities 
that  approximate  a  steering  level  w  ad.  Summarv  statistics  for  the 
dune  10,  197b  storm  displayed  in  figure  •!  are  presented  in  figures  8  to 
9.  These  data  provide  statistical  summaries  of  several  cell  character¬ 
istics  representing  either  the  values  averaged  over  the  lifetime  of  a 
cell  (average)  or  the  peak  value  obtained  during  the  cell  lifetime  (peak). 
Lifetime  data  are  presented  in  figure  10.  The  data  in  figure  8  depict 
the  statistics  of  the  highest  reflectivity  values  reached  hy  a  cell  dur¬ 
ing  its  lifetime.  The  data  are  drawn  from  a  sample  of  900  cells  whoso 
lifetimes  exceeded  1(>  minutes.  These  data  show  that  reflectivity  alone 
was  not  used  as  a  criterion  for  significance.  Over  cells  had  a  peak 
reflectivity  in  excess  of  80  dRZ  while  only  8  significant  cells  had  a 
peak  reflectivity  in  excess  of  f>0  dIC..  In  the  processing  used  to  obtain 
these  data,  significance  was  defined  based  upon  the  vertical  structure 
and  horizontal  dimension  of  the  cell  as  well  as  its  reflectivity.  A  high 
reflectivity  echo  that  was  observed  only  at  one  elevation  angle  was  not 
considered  to  ho  a  significant  cell. 

The  cells  detected  bv  the  processing  scheme  have  relatively  small 
areas.  The  average  cell  area  for  all  colls  and  for  significant  cells  do 
not  differ  as  indicated  in  figures  b  and  The  peak  cell  diameter  is 
of  the  order  of  A  km.  The  cells  tend  to  he  vertical  structures.  Statis 
ti catty,  the  area  of  a  cell  at  the  lowest  elevation  angle  does  not  differ 
from  the  area  of  the  cell  at  the  height  at  which  it  has  a  maximum  reflec¬ 
tivity  value.  Note  that  the  area  at  anv  height  is  defined  hy  a  contoured 
region  A  dR  below  its  peak  reflectivity  value  at  that  height.  At  another 
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height,  the  contour  used  to  define  its  area  will  have  a  different 
reflectivity  value. 

The  tracking  algorithms  track  each  cell  individually.  The  distri¬ 
butions  of  cell  velocity  averaged  over  the  cell  lifetime  ari  depicted  in 
Figures  8  and  9.  These  data  show  no  significant  differences  between  the 
velocities  of  all  the  cells  and  of  significant  cells.  The  individual 
track  velocities  may  differ  from  the  mean  (or  steering  level  wind)  velo¬ 
city  by  as  much  as  4  m/s  (rms  in  each  component  of  the  wind).  Observa¬ 
tions  show  that  the  deviations  are  not  entirely  random.  Larger  scale 
convergence  and  divergence  patterns  are  evident  in  the  cell  trajectories. 
On  a  smaller  scale,  it  is  evident  that  cells  affect  each  others  motion. 

A  tendency  has  been  observed  for  cells  to  follow  each  other  along  the 
same  track  even  though  they  may  not  have  initially  developed  along  the 
track . 

The  cells  persist  for  a  range  of  lifetimes.  On  average,  cells  with 
reflectivities  above  IS  dB  last  for  a  little  over  12  minutes.  Signifi¬ 
cant  cells  last  over  SO  minutes  on  average.  These  results  show  that 
ideally,  a  full  cycle  of  radar  observations  (a  volume  scan)  should  be 
acquired  in  between  four  and  five  minutes  to  get  more  than  two  observa¬ 
tions  of  a  cell  during  its  lifetime.  Practically,  a  longer  time,  six  to 
seven  minutes,  is  required  for  the  processing  algorithms  as  implemented 
under  this  contract.  Processing  speed  can  be  increased  by  preparing  a 
new  operating  system  for  the  Interdata  7-32  computer  but  this  was  beyond 
the  scope  of  the  contract.  The  motion  of  significant  cells,  if  detected 
early  during  its  development,  may  be  extrapolated  for  upwards  of  20 
minutes  before  they  disappear  indicating  that  short  range  forecasting  of 
cell  location  is  feasible. 

2.2  Cell  Detection 

The  cell  detection  algorithms  have  been  previously  defined  (Crane, 
1977,  1978)  and  will  not  be  detailed  again.  Flow  charts  for  the  process¬ 
ing  algorithm  are  presented  in  Appendix  B.  Briefly,  a  cell  is  a  region 
within  a  contour,  a  fixed  number  of  quantization  steps  below  a  local 
maxima  that  includes  no  other  cells.  For  most  observations,  a  quantiza¬ 
tion  step  of  1  dB  and  the  use  of  contours  3  dB  below  a  local  reflectivity 
maxima  seems  to  work  best.  The  quantization  step  and  contour  threshold 


were  empirically  established  by  t'rane  using  data  from  Wallops 

Island,  Virginia.  The  1  dB  step  and  5  dB  threshold  produced  a  detection 
probability  better  than  O.PS  on  these  scans  and  a  false  alarm  rate  of 
less  than  2  per  scan.  By  increasing  the  threshold,  the  false  alarm  rate 
was  reduced  but  at  the  expense  of  a  lowered  detection  probability. 

Cell  detection  is  performed  for  all  localized  reflectivity  maxima 
that  exceed  a  processing  threshold.  In  the  post  mission  processing  ver¬ 
sions  of  the  cell  detection  program,  attributes  are  obtained  for  all  the 
cells  detected  above  the  lowest  threshold  fixed  reflectivity  processing 
contour  (fixed  contour  or  echo  region).  The  streamlined  edition  of  the 
cell  detection  program  provides  output  only  for  the  most  significant 
cells. 

..3  Tracking 

An  entirely  new  tracking  algorithm  was  developed  during  the  period 
covered  hv  this  contract  primarily  for  application  to  the  post  mission 
processing  requirements  of  the  FAA  and  BuKoc  contracts.  This  program 
was  subsequently  modified  for  real  time  application  under  this  contract. 
The  real  time  version  generates  volume  cells  and  volume  cell  tracks. 

The  same  tracking  algorithm  is  used  to  generate  the  3-dimensional 
cell  from  the  successive  azimuth  scans  within  a  volume  scan  sequence  and 
to  track  the  volume  cell  in  time.  1  ach  new  volume  cell  is  identified 
using  cell  data  from  a  single  azimuth  scan  that  cannot  be  associated 
with  other,  previously  established  tracks.  Association  is  attempted  with 
the  newly  detected  volume  cell  on  each  successive  azimuth  scan  until  the 
track  is  terminated.  A  track  is  terminated  when  no  new  data  are  obtained 
for  a  volume  cell  during  a  complete  volume  scan  cycle. 

Association  is  established  using  the  location  of  the  cell  on  a 
azimuth  scan  as  compared  with  the  volume  cell  location  extrapolated  to 
the  rime  of  the  azimuth  scan  together  with  the  differences  m  the  heights 
of  the  last  observed  data  and  the  current  data,  differences  in  cell  areas, 
and  differences  in  ref leet i v i tes .  A  measure  of  the  goodness  of  an  asso¬ 
ciation  between  a  cell  and  track,  is  established  for  each  possible  track, 
cell  combination.  The  final  cell  track  pairings  are  those  that  minimize 
the  sum  of  the  measures  (maximize  the  goodness  of  association)  for  all 
the  cells  and  tracks  that  may  be  possibly  paired  during  the  association 


process.  The  set  of  colls  and  tracks  that  m.  v  possihl>  he  paired  are 
defined  as  a  cluster. 

The  volume  cell  track  is  the  priman  entitv  maintained  by  the 
tracking  program  and  successive  radar  obsorv;  t ions  are  associated  with 
the  expected  position  of  the  cell  along  its  rack  at  the  tine  of  the 
radar  observation.  The  tracking  algorithm  was  developed  in  this  nanner 
to  allow  the  use  of  data  from  more  than  one  radar  system  since  all  that 
is  needed  for  the  association  algorithm  is  the  location  and  time  of  the 
cell  centroid  as  reported  bv  each  radar  together  with  other  attributes 
such  as  reflectivity,  tangential  shear,  etc. 

Two  velocities  arc  maintained  for  each  cell  track  -  an  instantaneous 
velocity,  the  difference  in  cell  centroid  location  divided  by  the  time 
interval  between  observations,  and  a  smoothed  velocity  obtained  by  low 
pass  digital  filtering  of  the  raw  velocity  data.  The  coefficients  in 
the  filter  were  established  by  trial  and  error  using  a  large  data  sample. 
The  initialization  of  the  track  velocity  for  each  cell  is  important  due 
to  the  extrapolation  process  used  for  tracking.  Ixperiments  with  the 
tracking  program  show  that  adequate  results  are  obtained  if  a  zero  velo¬ 
city  is  used  for  the  initial  velocity  estimate  but  better  results  can  be 
obtained  if  an  estimate  of  the  steering  level  wind  is  used  for  the  initial 
tracking  velocity.  The  program  automatically  updates  the  initial  velocity 
estimate  after  processing  sufficient  data  to  establish  a  stable  estimate. 
The  measure  of  success  of  the  tracking  program  was  taken  as  t  10  rate 
required  for  the  smoothed  track  velocities  to  stabilize. 

The  cell  detection  and  tracking  programs  were  initially  developed  to 
process  a  large  number  of  cells,  ip  to  512  active  cell  tracks  at  any  one 
time  and  to  calculate  upwards  of  30  attributes  for  each  cell,  cluster, 
and  fixed  contour  (echo')  region.  Processing  this  amount  of  data  is  not 
possible  in  a  real  time  environment  with  storage  and  display  limitations. 
The  basic  algorithms  for  cell  detection  and  tracking  have  been  maintained. 
The  number  of  tracks  to  be  processed  lias  been  reduced  by  increasing  the 
reflectivity  threshold  for  processing  and  by  incorporating  the  tangential 
shear  information  in  the  decision  process  for  saving  the  most  important 
12  to  IP  cells  of  50  or  more  cells  that  exceed  the  reflectivity  process¬ 
ing  threshold.  Further  storage  savings  have  been  accomplished  by  reduc¬ 
ing  the  number  of  attributes  for  each  cell. 


The  cell  detection  aiul  tracking  programs  were  not  coded  in  effi¬ 
cient  manner  for  operation  tin  the  f.TSl  Interdata  7-32  computer.  Inten¬ 
sive  program  revision  was  performed  to  reduce  the  number  of  subroutine 
calls  and  to  revise  the  addressing  procedure  to  reduce  the  time  required 
to  fetch  or  store  a  variable.  The  result  is  a  streamlined  cell  detection 
and  tracking  program  that  will  handle  a  reas.  liable  number  of  active  cells 
(.luring  the  time  required  for  an  azimuth  scan.  Specifically,  over  100 
cells  can  be  detected  and  processed  in  less  than  32  seconds  using  the 
Interdata  7-32  programs.  The  processing  time  can  be  reduced  further  by 
dynamically  varying  the  reflectivity  processing  threshold  to  maintain 
fewer  than  say  JO  detected  cells  but  this  has  not  1il.ii  necessary. 

The  programs  to  fetch  the  raw  radar  data  and  prepare  the  data  for 
use  in  the  cell  detection  program  are  included  in  Appendix  (!  and  I).  The 
real  program  is  designed  to  permit  data  gathering  and  cell  detection  on 
alternate  scans.  Real-time  processing  using  the  cell  detection  algor¬ 
ithms  on  every  azimuth  scan  is  possible  with  the  Interdata  7-32  computer 
but  will  not  occur  unless  considerable  effort  is  expended  to  develop  a 
new  operating  system  for  the  computer  tailored  to  use  the  interrupt  and 
background/ foreground  processing  capabilities  of  the  computer  to  provide 
quasi -real-time  cell  detection  and  simultaneous  real-time  radar  data 
acquisition,  averaging  and  storage.  Sufficient  time  is  available  for 
all  the  programs  to  operate  on  all  the  data  from  a  volume  scan  within  the 
time  of  a  volume  scan  but  the  processing  of  the  data  from  the  lower  ele¬ 
vation  angles  will  lag  behind  data  acquisition  and  only  catch  up  on  the 
higher  elevation  angles. 
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1  I'rovoss  i  ng  Opt  ions 

Two  programs  ovist  tor  opor.it  ion  on  flu-  I  INI  ,  V  post  mission  pio 
oossing  program  is  available  valh\l  OKWI  that  roads  proviouslv  pvopaied 
data  from  tiio  disk  ami  pot  forms  .ell  ,lotoot  ion  ami  tr.uking.  fins  pro 
g ram  is  rot  intended  for  real  t  into  oporat  ion  ami  threshold  level-  ma>  !>o 
reduced  to  allow  detection  of  a  largo  number  of  rolls.  I  ho  program  is 
intended  for  post  mission  Jala  aimlisis  when  t  i mo  i-  not  at  a  ptominm  ami 
larger  amounts  of  output  ran  So  handled  h\  t  ho  usor.  The  input  data 
must  ho  prepared  for  storage  on  tho  disk  using  a  moditioation  to  tho  I  St 
programs  generated  S'  Kavtheon.  Iho  prop  ram  r.  called  I  III  and  is  listed 
in  Vppemliv  1'. 

l'ho  operational  program  uses  tho  same  .oil  dotootion  and  tracking 
subroutines  but  is  railed  from  a  modified  version  ot  tho  Uavthoou  pro 
v '  .led  1  SI  programs.  'hose  programs  store  tho  data  on  one  aimuth  n  an 
ami  process  the  data  on  tho  ne\t.  Pata  I  rots  tho  tii-.t,  third,  ami  tittli 
srans  are  processed  during  tho  second,  touvth.  and  sivth.  I'm  am:  the 
seventh  sran  the  displavs  are  prepared  and  tho  programs  ,lto  re  i  n  i  (  i  a  I  i  o 
for  the  no\<  ru  le  ot'  the  seven  s,  an  sequence. 

All  the  programs  generated  inulor  this  lontraet  tor  ir.o  on  the  Intor 
data  computer  are  listed  in  Vppemliv  f  and  P.  iKVM  is  the  post  mi 

main  program  which  rails  fi'MV'K  and  IKVfk.  tho  output  t  rout  this  pro. -ram 
is  stored  on  disk  tor  subsequent  listing.  Ibis  program  requires  the  u.o 
of  a  preprooessiug  program  ITKI'Nt'  (modified  tor  lUP  written  bv  VITI 
personnel  ami  listed  in  Vppendiv  P.  Vho  real  time  program  in.  hides 
modified  versions  of  the  VSf  programs  supplio.l  h\  Kav t boon .  Tho  "i.'.litio 
programs  are  TSIMAtV  and  the  subroutines  KMVT  ami  ITNI'Sv  .  Ilio  v'v'Mv'U 
and  TUATh  subroutines  are  vailed  from  UMVT.  Vboso  programs  are  present  e 
in  Vppemliv 

A.  '  Toll  I'otovt  ivMl 

l'ho  eel  1  dote,  tiv’n  subroutines  fONlA'H  and  I'l  Vkl'  wore  -  uV -t  anl  t  a  1  1  \ 
modified  to  reduce  proooss :  ng  time  spent  in  addressing  the  data  ari.n  s, 
l'ho  pvinvipal  modi  float  von  was  to  vhango  all  the  arrays  to  single  dvmew 


••  lona  I  arrav  s  and  to  explicit  Iv  perform  tho  address  calculations  in  the 
program.  In  this  wav,  multiple  references  to  the  same  array  location 
would  not  involve  the  time  consuming  recalvulat  ion  of  tho  address  tor 
each  arrav  ret'erence, 

fONTOK  was  extensively  modified  to  remove  the  fixed  contour  attri¬ 
bute  generation  algorithm,  contour  data  are  still  prepared  hut  in  the 
a.'imuth  strobe  format  used  hv  the  TS1  displav  program.  No  radial  to- 
radial  association  is  re>|uired  for  this  streamlined  version  of  the  pro¬ 
gram  s i gn i  f i cant  1  v  decreasing  the  length  and  complex!  tv  vi  1'  the  I'rogram. 

I'l  AM'  has  been  changed  onlv  to  accept  the  mod  it' ted  system  of 
addressing  and  to  select  onlv  the  lb  most  significant  cells  for  further 
processing.  Koforeneo  to  t  hi'  lived  contour  wlent i f ieat ion  tag  was  also 
removed  from  I'l  \KP  since  the  tags  were  produced  in  the  association, 
attribute  generation  logic  of  the  c'ONTOR  subroutine  which  was  removed 
tor  (lie  Interdata  version  of  the  program. 

'  '  I rack i ng 


fhe  now  tracking  prog'.im  consists  of  the  subroutine  I'RAi'K  which 
calls  fOMPtlv  to  perform  the  cell  to  t  rack  assoc  i  at  ions .  five  subroutine 
fi'MI'XK  searches  the  track  list  and  the  cell  I  i  t  from  the  last  scan,  and 
finds  ,iil  pos'.ible  pairs  for  which  the  goodness  of  fit  measure  does  not 
exceed  a  prose  I  ei  t  ed  threshold  lhe  measure  i-  given  hv 


iV,lrl-*n 
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where  M  is  the  measure, 
and  i  s  ret  1  *  t  i  v  i  I  v  . 
lhe  weights  lv  were  set 
list ed  i n  I ub 1 e  I  .  I  he 
pan  mgs  with  a  measure 


t  he  sul  s»  v  i  pt  s  f  refer  t 
\  area.  II  height  and  \ .  \ 
bv  trial  and  error.  lhe 
be-,  t  cell,  t  rac k  pa  i  r  l  ng 

great ei  than  M  are  not 

P 


o  cell  and  T  to  track 
are  centroid  locations, 
current  values  are 
has  the  lowest  measure; 
a  I  I  owed . 


••ever.tl  pairings  are  possible  in  a  cluster  of  cells.  Subroutine 
KlSi'IV  selects  the  best  set  of  cell,  track  pan's  in  a  cluster.  The 


attributes  are  updated  in  MKAk  which  is  called  from  filMlWR  if  there  is 
no  tell  .luster  or  from  111  SOI  V  if  a  fluster  exists.  fhe  subroutine  HTRAk 


is  also  called  to  store  cell  data  in  the  VR  array  each  time  a  track  is 
updated.  This  information  is  used  in  calculating  the  measure  M.  The  VR 
array  data  are  either  from  the  lowest  elevation  angle  on  which  the  track 
was  observed  or  from  the  last  elevation  angle  at  which  it  was  observed. 
The  measure  used  to  evaluate  the  cell,  track  pairing  is  the  minimum  mea¬ 
sure  obtained  using  either  the  last  or  lowest  elevation  angle  data. 

At  the  end  of  a  volume  scan  cycle.  STRAK  is  called  to  calculate  the 
attributes  and  to  output  the  track  data.  Only  12  tracks  are  output  from 
STRAK  in  the  operational  version  although  a  maximum  of  tra  ks  are 
maintained  at  any  one  time.  The  list  of  the  16  attributes  maintained 
for  each  track  is  given  in  Table  /. 

tabu:  1 

TRACK  1  NO,  Will  OUTS  ,  MAXIMUM 
MI  ASURi:  AND  VIII.OCITY  1- 1 1.THR 


••  wx 

1/10 

IdlC.')'1 

WA 

1/25 

(km)  ’ 

WH 

1/.' 

(km)  1 

\ 

1/S 

tkm) 

M 

I’ 
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Velocity  filter: 

-  - 

aV  *  * 

*  VN 

~  a 

S-J  +  h  iv)i  * 

'•  VT 

V  velocity:  N.l 

:  refer  to  Northward  and  1, 

ponent  s 

Ax, Ay  -  change  in  position  between  volume  scans 
At  =  time  between  scans 
V  -  average  velocity,  all  cells 


tabu; 


TRACK  ATTRIBUTES 


Observation  Time  (seconds  from  start  of  year) 

East  Location  (km) 

North  Location  (km) 

Average  Reflectivity  (JBZ) 

Volume  (km-') 

Peak  Reflectivity  (dRZ) 

Height  of  Reflectivity  Peak  (km) 

Reflectivity  at  Lowest  Elevation 

Angle  (JBZ) 

Area  of  Cell  on  Lowest  Elevation 

Angle  (km") 

Height  of  Cell  at  Lowest  Elevation 

Angle  (km) 

Reflectivity  at  Cell  Top  (dBZ) 

Height  of  Cell  Top  (km) 

Track  Identifier 

Smoothed  Track  Velocity  East  (m/s) 

Smoothed  Track  Velocity  North  (m/s) 

Integrated  Tangential  Shear  (m/s/km) 


4.  PROGRAM  STATUS  AND  RECOMMENDATIONS 


The  processing  programs  are  operational  on  the  Interdata  ?-5’  com¬ 
puter.  Experience  must  now  he  gained  in  using  the  system  for  the  obser¬ 
vation  of  weather.  A  number  of  parameters  (Table  11  were  set  in  the 
program  on  the  basis  of  our  experience  with  the  data  obtained  from  our 
contracts  with  the  EAA  and  BuRec,  It  is  anticipated  that  a  different 
radar  system  operating  in  a  different  environment  may  need  a  different 
set  of  parameters.  These  parameters,  such  as  the  tracking  measure 
weights  and  the  constants  in  the  velocity  smoothing  filter,  are  readily 
changed  in  the  program.  Experience  with  a  larger  data  set  is  required 
to  obtain  the  best  estimate  values  for  the  parameters. 

It  is  recommended  that  the  cell  detection  and  tracking  programs  be 
used  on  available  data  to  develop  the  required  operational  experience  to 
adequately  use  the  new  displays.  The  new  output  attributes  are  in  a 
form  that  may  be  readily  adapted  to  objective  warning  and  forecasting 
systems . 
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APPENDIX  A 

Operating  Instructions* 
CSS  Files* 

TET  Files* 

Task  Establishment  Maps* 
Definition  of  Variables 
List  of  Arrays 
List  of  Common  Blocks 


*Listings  for  Real-Time  and  Post-Mission  Programs 


A-l 


REAL-TIME  PROGRAM  OPERATING  INSTRUCTIONS 

TO  EXECUTE  REAL-TIME  CELL  DETECTION 
AND  TRACKING  PROGRAM  (ERT) 

1.  Compile  each  subroutine  file: 

OCOMPILF.  ERT1 : Fi le  Name* 

2.  Create  an  object  file  for  each  of  the  compiled  subroutine  fit 

DCAI.OBJ  F.RT1  :File  Name* 

3.  Delete  old  task  file: 

DF.  ERT1 : CRANE. TSK 

4.  Create  a  new  task  file: 

ESTAB  F.RT1 :  GARY 

5.  Dispose  I/O  devices  hy  editing  the  CSS  file: 

F.DITR  CRF.AI..CSS 
0.  Execute  the  program: 

CREAl. 


*l.ist  of  subroutine  file  names: 


TSF.MAIN 

ATRAK 

TSKDATA 

BTRAK 

REALTM 

COMPAR 

INPARM 

RESOLV 

CONTOR 

COMB IN 

PEAKD 

STRAK 

TRACK 

A- 


REAL  TIME  CSS  FILE 


* 

07/19X79  16:56:49 

***LISTING  FOR  CRERL  CSS 
SN 
SJOB 

L  !  BO,  ERT1 :  TSEPRRAM 
T  .  BG 

flS4/  ERT1 :  TSEPRRAM.  DAT 
ST  /  @1#  @2,  93,  94,  05/  96,  99 
L~Cre7ERTrrTS£RTC~RTi: 

SE  PR  2/0/ 1/145.  0 
L  CRfiNE/  ERT1 : CRRNE 
T  CRRNE 
RS5/  CON : 

RS5/  CON : 

fissrERTrrpcr - 

RS9,  ERT1 : DRT1024 
RS10,  PPRI : 

RS7/  PACK :  TSEPPI.  PAR 

ASS/ PR: 

SIFNNULL  07 

- f&5T97 - 

$ENDC 

SIFNULX07 - 

RS5/  ERT1 :  TSEDEF.  PRR 
*ENDC 

RS4,  ERT1 :  TSEPARRM.  DAT 

ST - 

♦EXIT 


REAL  TIME  TET  FILE 


07/19/79  16:55: 46 

♦♦♦LISTING  FOR  ERT1 : GfiRV.  TET 
*N 
LOG 

JOB  ERT 
REMOTE 
ES  TASK 
MXSPACE  2800 
OPTIONS  F 
GET  400 

PRIORITY  10# 10 
TCOM  TASK/2/RW 
TCOM  2ST0RE/3/RW 
TCOM  D0NE/4/RM 
TCOM  EXTRA/5/RW 
IN  ERT1 : TSEMAIN 
IN  ERT1 : TSEDATA 
IN  ERT1 : REALTM 
IN  ERT1 : INPARM 
IN  ERT1 : CONTOR 

TN~£RTT7FE~ftK5 - 

IN  ERT1: TRACK 
IN  ERT 1 : ATRAK 
IN  ERT1 : BTRAK 
IN  ERT 1 : COMPAR 
IN  ERT1: RESOL V 
IN  ERT17STRRR 
IN  ERT1 : COMB IN 
RESOLVE  ERT1 : TSERTL 
EDIT  FVRTL 
BU  TASK# ERT1 : CRANE 
MAP 
END 


i  ■  1 1  i_«v  -  i_-3  i  noi_  l  .sncie.i'*  I  UUMy  I’lHf"'  l\r./U,  !  1MJ:  IAi>K.  I:S  I  AlU.  ISHMLNT  MM1 

DATE  07/19/79  TIME :  16  52:01 
JOB : ERT 

****  CTOP=0242FE  UTOP=023E40  MIN  PARTITION-  144  75K  **♦* 


PROGRAM  SEGMENTS: 


SEG 

TVPE 

NAME 

SIZE 

0 

IMPU 

143.  75K 

3 

TCOM 

2ST0RE 

TCM 

4.  25K 

5 

TCOM 

EXTRA 

TCM 

1.  25K 

15 

RTL 

TSERTL 

RTL 

7.  75K 

PROGRAM  LABELS 


000100 

TSEMAIN 

000620 

TSEDATA 

003838 

ERTRMAP 

01OA80 

I CLOCK 

010B70  COS 

010CB8 

SIN 

01ODFB 

T 

010EEB 

R 

011000 

R 

011110 

SORT 

011210 

ALOG10 

0112D8 

ALOG 

011430 

EXP 

011588 

FLOAT 

0115A8 

IFIX 

0115F8 

M0D2 

011640 

ABS 

011678 

IABS2 

011698 

H 

011708 

.  S 

011798 

01 

0117A0 

.  0 

011D48 

.  MES 

TASK  ENTRV-POlNTS: 

000664 

FPRDSC 

003S7C 

RMAP 

004204 

CONTOR 

004 9R4 

PEAKD 

009BEC 

TRACK 

00A0A4 

ATRAK 

00A87C 

BTRAK 

00AAF4 

COMPAR 

00BFAC 

RESOLV 

00EFCC 

STRAK 

010764 

COMB IN 

010RS2 

1 CLOCK 

010B72 

COS 

010BSA  .  COS 

010CBA 

SIN 

010CD2 

SIN 

010DFA 

.  I 

010EEA 

.  A 

011002 

.  R 

011112 

SORT 

01112A 

.  SORT 

011212 

ALOG10 

01122A 

.  AL0G10 

0112DA 

ALOG 

0112F2 

ALOG 

01131A 

.  L0GRT2 

011432 

EXP 

01144A 

EXP 

01158A 

FLOAT 

0115AA 

IFIX 

0115FA 

M0D2 

011642 

ABS 

01167A 

IABS2 

01169A 

.  H 

01170A 

.  S 

0117R2 

0 

0II7B2  "TOI  0'1192E~7OI  011D4A  ME5 


LOCAL  COMMON  BLOCKS 

011D68  CALB 

011D78  MUSIG 

011D80 

READZ 

011D8S 

SECTOR 

~n0llB90~REFI. 

8129B0  CAL 

0129C0 

RUNSUM 

~  0139D0 

TLIS 

0139E8  CALR 

016208  SWITCH 

020098 

AZ2 

0200B0 

AZM 

0200C8  PNTRS 

0200D8  INTL 

0200E8 

ZLOOK 

020258 

ECONST 

020260  MAPPAR 

020278  CNT 

020288 

DATA1 

020610 

DATA2 

021198  DRTR3  ~~ 

0214A0  NVLIS 

0214B0 

FILTER 

0214C0 

KTR 

0214C8  CDRAVS 

0222D0  CONST 

022300 

VPARM 

02230S 

D'v’AL 

0223I0CNT  RS 

“022320TLGS 

~  02233D 

TMRX 

02233S 

PWORK 

0223E8  FIXED 
023AF0  RSLV 

LIBRARV  ENTRIES: 

0223F8  PRSTOR 
023C78  COMB 

023AE0 

THRESH 

023AE8 

CONPK 

0F0802  U 

0F0052  .  V 

0F006A 

P 

0F00FA 

Q 

0F0182  @R 

0F02C8  T?H 

0F17E2 

CONMSG 

0F1 86 A 

FL0RT2 

TASK  COMMON  BLOCKS. 

030000  ZSTORE  050000  EXTRA 

END 


POST-MISSION  PROGRAM  OPKRATINO  INSTRUCTIONS 


to  kxkcutk  post-mission  cku  pktk.ction 
anp  tracking  program  (CRANK! 


Compile  each  subroutine  file: 

PCOMPIl.K  KRT1 : Ki  le  Name* 

!.  Create  an  object  file  for  each  of  the  compiled  subroutine  fil 
OCAI.OW  TRT1  :Ki  le  N;une* 

>.  Pelete  old  task  file: 


PT  t-RT  1 : CRANK. TSK 


Create  a  new  task 


f  i  le : 


l-STAR  CRT  1 :  CRANK. 


S.  Pispose  I/O  devices  by  editing  the  CSS  file: 
KPITR  CRANK. CSS 
K.xceute  the  program: 


CRANK 


list  of  subroutine  file  names 


CRANK. 

I  NT  ARM 

CONTOR 

PKAKP 

TRACK 

ATRAK 

BTRAK 

COMPAR 


RKSOI.V 

COMBIN 

STRAK 


POST-MISSION  OSS  PI  If 


"•listing  of'  tin'  ORANP.OSS  Pile- 
$OOR 

sr  pa  i  i -jo.;  o..'  o 

!  I  1B.PRT1  rTSPRlT  .  KTI 
1  ORANP . PRl'l : ORANP 
T  ORANP 
\S0. MM  I  : 

ASS, PR: 

AS4.MU.1  : 

ASt'.FRT!  :Usor  IV  fined  Output  Pile* 

AS”.  I:RT  l  :Pser  Pofined  Pal  ibrat  ion  Pile*'* 
ASS ,  PRT1  :  User  Defined  Input  Pat  a  Pile*** 
ST 

$TPRM.TOR 
£p\  1  T 


File  to  which  program  output  is  to  he  sent 

'Pisk  Pile  containing  values  for  the  parameters:  10UT.  PBB,  MAW 
MANS,  Sl.OPF.  01.PATA* 


Pisk  file  generated  hy  program  "PRT” 


POST-MISSION  TET  FILE 


TET/22  R02-02 
JOB  ERT 
REMOTE 
ES  TASK 

r''  MXSPACE  2800 
OPTIONS  F 
GET  400 

r  PRIORITY  10,  10 
TCOM  TASK/2/RW 
TCOM TZST0RE/2/RW 
r  TCOM  D0NE/4/RW 
TCOM  EXTRfl/5/RW 
IN  ERT 1 : CRANE 
r-  IN  ERT1.1NPRRM 
IN  ERT1 : CONTOR 

/*»  IN  ERT1:  TRACK 
IN  ERT1.HTRRK 
IN  ERT 1 : BTRHK 

S'  IN  ERT1 :  COMPRR  — ~ - 

*  IN  ERT 1 : RESOL V 


f  IN  ERT1: COMBIN 

RESOLVE  ERT1TSERTL 
EDIT  FVRTL 

r  BITT ASK,  ERTl :  CRANE 
MAP 


r 

POST- 

MISSION  TASK 

TSTARI.l SIIMPNT  MAP 

- ,, 

0S32MT  TASK -ESTABLISHMENT  LORD  MRP 
DRTE  07/19/79  TIME  1310  27 

( 

JOB  ERT 

*+**  CTOP«01E5FE 

UTOP=01F1F8  MIN  PARTITION-  121 

50L  ♦*** 

r 

PROGRAM  SEGMENTS 

\  C 

SEG  TYPE  NRME 

SIZE 

|  , 

!  j 

0  IMPU 

120  50K 

\ 

3  TCOM  2ST0RE 

TCM 

4  25K 

c 

5  TCOM  EXTRA 

TCM 

2.  50K 

15  RTL  TSERTL 

RTL 

7  75K 

c 

PROGRAM  LABELS: 

h 

( 

000100  ERTRMAP 

00E850 

COS 

00E998 

SIN 

00EAD8 

i 

SOEBCSFr 

00ECE0 

R 

'  00EDFB 

SORT 

OOEErO 

ALOG 10  1 

c 

00EPB8  ALOG 

00F110 

EXF 

00F268 

FLOAT 

00F288 

IFIX 

00F2D8  ABS 

80F310 

IABS2 

00F330 

s 

O0F3C0 

0 

1 

00F968  MES 

i  . 

TBSKENTR\T-;POINTS 

-  . 

■ 

( 

001FD4  CONTOR 

002774 

PEAKD 

0079BC 

TRACK 

007F.74 

ATRAK  ; 

00864C  BTRAK 

0O88C4 

COMPAR 

009D7C 

RESOL V 

O0CD9C 

STROK 

< 

00E534  COMB IN 

00E852 

COS 

00E86A 

COS 

O0E99H 

SIN  H 

O0E9B2  SIN 

00ERDA 

.  I 

00EBCA 

A 

00ECE2 

R 

B&EDF5  SORT 

00EE0A 

SORT 

0OEEF2 

ALOG10 

00EFOA 

ALOG10  ' 

c 

00EFBA  ALOG 

80EFD2 

ALOG 

00EFFR 

LOGRT2 

00F112 

EXP 

00F12A  EXP 

OOF 26 A 

FLOAT 

00F28A 

IFIX 

0OF2DA 

ABS  f 

00F312  IRBS2 

00F332 

s 

00F3C2 

0 

O0F3D2 

Oi 

C 

00F54E  01 

0OF96A 

MES 

1; 

( 

LOCAL  COMMON  BLOCKS 

fj 

00F988  TLIS 

00F9A0 

CONST 

00F9D0 

SWITCH 

019860 

AZ2  j  ] 

( 

019878  A2M 

019890 

REEL 

01A4B0 

PNTR5 

01A4C0 

intl  j; 

01A4D0  ZLOOK 

01A640 

ECONST 

01A648 

MAPPAR 

O1A660 

CNT  1’ 

0iA9FS“ DATA2 

01B580 

DAT  A3 

01B888 

NVLIS  | 

| 

01B898  FILTER 

01B8A8 

KTA 

01B8B0 

CDRAYS 

01C6B8 

WARM  }, 

V 

01C6C0  DVAL 

01C6C8 

CNTRS 

01C6D8 

FLGS 

01C6L8 

T  MAX  i 

01C6F0  PWORK 

01C7A0 

FIXED 

01C7B0 

PRSTOR 

0 IDE 98 

THRI’SH  ,, 

( 

01DEA0  CONPK 

01DEA8 

RSLV 

01E030 

COMB 

I 

t  , 

« 

LIBRARY  ENTRIES 

1 

0F0002  .  U 

0F0052 

V 

0F006A 

P 

or OOF A 

Q 

1  1 

0F0182  9R 

0F02CB 

OF 186 A 

FL0AT2 

1 

1  ASK  COMMON  BLOCKS 

i 

030000  ZSTORE 

050000 

EXTRA 

.  0 

L 

1  ND 

i ... 

VAR  1  AH  1,1  1 . 1  ST 


('RANI:  f,  Ri-AI.TM 
+TI.1  -  lower  threshold  t *.1 H IT V 
♦TI.Z  -  higher  threshold  (dRZ) 

+RQIIANT  -  threshold  quantization  factor  (dBZ) 

♦  U  f, *Z.X  -  intermediate  values  in  dHZ  conversion 
+ 1 l Oil  -  reflectivity  offset  to  insure  positive  values 
* PB  -  calibration  constant 

*IUT\T'I,  -  constant  used  in  calculation  of  velocity 
*H1TVAR  -  constant  used  in  calculation  of  variance 
*1?4W0RP  -  square  of  variance 
*i  I.1VAT  -  raw  elevation  (encoder  units') 

*ll,i:VAA  -  elevation  in  degrees 

+ 1 IU.SN  •  elevation  in  degrees 

+T  -  time  (seconds  from  start  of  year) 

♦AZ  -  raw  azimuth  (deg) 

*AZT  -  azimuth  in  degrees 
*.\  -  a,  iinuth  in  radians 
*AZCHK  -  azimuth  in  deg.  ♦  ,%S‘> 

+K  -  radial  counter 
*R(1NA  -  start  azimuth  (deg) 

*I:.NPA  -  stop  azimuth  (deg) 

*S1NA  -  sin(A) 

*COSA  -  cos (A) 

+NA  -  check  on  first  azimuth  of  new  scan 

♦NAC  -  offset  multiplier  for  arrays  in  t'ONTOR  !i  I'UAKP 

♦I PAY  -♦PAY  -  Julian  date  data  collected 

*1  HR  -  hour 

+IMIN  -  minute  • time  scan  begins 
+  TSI-C!  -  second 

♦pri.TAZ  -  ar, imutn  increment  (radians) 

*Plll  -  elevation  angle  (radians) 

VOSPIIl  -  cos  (PHI) 

♦S1NPHI  -  sin(Plll) 

*1  .AR I'll  -  (>.48571-  -5  (km'  ) 

VOSPIIZ  -  (:0SI’1U-'*1'AUTH 

*AZI,AST  -  azimuth  of  previous  radial  (deg) 


+  -  INTldbR*: 
*  -  RP.Al. 


A-  l  0 


CONTOR 


♦  I FLAG  -  intermediate  print  flag 
*TI.S  -  lower  threshold  fTl.T  in  CRAM-) 
*TATRMN  -  test  on  area 

♦NHMC  -  array  addressing  offset 

♦NCI:!.  -  cell  counter 

♦IRM/VX,  +JMAX  -  array  limits 

♦NF.M1  -  address  variable 

♦NCI.  -  maximum  number  of  positions 

♦TI.l  -  lower  threshold 

♦  TF.M  -  event  counter  on  first  threshold 
♦II:M2  -  event  counter  on  second  threshold 
♦1PR  -  peak  start  location 

♦IF  -  peak  counter 
♦IFVFNT  -  event  number 


+  -  INTFGF.R*’ 
_  RliAL— . .  . 


\-U 


PEAKP 


+\BAPR,+NCAPR,+NBKA,+NCKA  -  address  variables 
+NAX.+NA  -  radial  counter 
+  LM  -  number  of  variables  in  UP  (Array) 

+  LMDP , +NAN , +NAN 1 , +LMM , + I DX , +NCLM , +  LDBM , +LPX , +NPDP , + 1 P  2  -  address  variables 

*FQUANT  -  threshold  quantisation  factor 

+  KOFST ,  +  I.IMT , +N IPP , +MXTR , +  KMAX  -  Array  limits 

♦ IF  -  event  number 

+THM  -  number  of  events  on  radial 

+ IEA,+KIE ,+KIHM  -  address  variables  on  IEM 

+ICEST  -  event  start  position 

+ICESP  -  event  end  position 

+.TER,  +  KA,  +  KB  -  address  variable  on  event  and  peak 

♦IPL  -  peak  start  location 

+IP  -  peak  stop  location 

+NTHRES  -  threshold  counter  for  peak 

+  l,PB  -  dB  below  peak  value  used  to  define  peak 

+ 1R1  -  range  to  peak 

♦III  -  dBZ  value  at  IR1 

+IT  -  dBZ  above  threshold  at  IR1  for  I.PB  thresholds 
+.JMXPM  -  limit  on  IT 

♦  KA  -  address  variable  on  peak  thresholds 

+  IPT  -  number  of  thresholds  associated  with  peak 
+dR  -  limit  on  TPT 
+  IRGN  -  first  position  in  event 
+  1ND  -  last  position  in  event 
+  111  -  reflectivity  values  within  event 

+KA  -  address  variable  on  peak  threshold  values  within  event 
+ TMXJMX  -  limit  on  number  of  contours  per  radial 
+  IRF.G  -  start  or  stop  range  of  contour 
+ 1  Pl:  -  contour  counter 

+  fAPDR  l'i +IF.QL  -  address  variables  on  contour  threshold  and  number 

+KC,+KA,+KZ  -  address  variables  on  contour  thresholds 

+ TCV 1 , ,  +TC VM , + TCV l.B  -  threshold  values  on  a  peak 

+NPC  -  number  of  contours  by  threshold  on  this  radial 

+NP1.  -  contour  counter  by  threshold 

+ 1 1  IBM  -  start  range  of  contour 

♦ IHB  -  IHBM  +  1 

+  II1P  -  stop  range  of  contour 

+K,+KY,+KZ  -  address  variables  on  next  threshold 

+  1.PF.  -  number  of  contours,  this  thrchold,  on  radial 

+  I.PI,  -  contour  counter,  this  threshold 

+MPCFI.  -  IP  number  for  possible  cell,  this  threshold 

+TATC  -  IP  TATR(NIYFE)  points  to 

+.IE1  -  first  event  previous  radial 

+.IE2  -  last  event  previous  radial 

+.IF.M  -  address  variable  on  events,  previous  radial 
♦IPB  -  number  of  peak  thresholds 

♦  KB, +KBR,  +  KBA,  +  KRC  -  address  variables  on  i ontcur  thresholds 
+TBVI.  -  threshold  value  on  a  peak 

♦  NP2  -  numb  t  of  contours  on  a  threshold,  prei ioio  radial 
♦NP1  -  cont  >ur  counter,  bv  threshold,  previous  radial 

♦  -  INTEGER  *2  *  -  RFAI. 


1 

U 


VI. 


PEAKD  (continued) 


♦LPCEL  -  ID  number  for  possible  cell,  previous  radial 
♦IEQL  -  ID  TATR(LPCEL)  points  to 
♦MPC  -  equals  NPCEI.  if  associated 

♦  JEQl,  -  next  higher  threshold 

+JM  ,  +JN2,  +.TX3 ,  +.l\’  l ,  t.TNS  -  address  variable  on  peak 

♦1ST  -  start  range  of  contour 

♦ISP  -  stop  range  of  contour 

*R  -  area  per  azimuth  degree  at  peak  range 

♦  III  -  reflectivity  at  that  range 
*RU  -  reflectivity  weighted  area 
*SAZ  -  sin(nzimuth) 

*CAZ  -  cos (azimuth) 

♦  KNN.  +  KN  -  address  variable  on  NPCEL 

♦IMPX  -  address  variable  on  second  threshold  of  peak 

♦  IND,+  INDX, ♦LNX  -  address  variable  on  I.PCEI, 

♦IN.+  1EQL  '  address  variable  on  LPCEL 

♦JEQL  -  area  address  of  NPCEI.  pointed  to  by  area  of  I.PCEI. 
♦IPTT  -  number  of  peak  thresholds 

♦KTl .♦KTA,*KTR  -  address  variables  on  contour  threshold 

+NPCT  -  number  of  contours,  this  threshold 

♦TEQL  -  threshold  value 

♦INDXT  -  address  variable  on  NPCEI. 

♦NIMN  -  number  of  possible  cells  tested  so  far 
+NIDP  -  limit  on  number  of  possible  cells 
♦IE  -  event  number 

♦IPT  -  number  of  peak  thresholds  in  event 
♦NPC  -  number  of  thresholds  this  peak 
+NCVM  -  threshold  value  on  peak 

♦  NP1,  -  contour  counter 


+  -  INTEGER* 2 
*  -  REAL 


\  - 1 


I’lAKI*  (continued) 


♦I  -  each  position  out  the  radial 

+  IA  -  address  variable  on  1 

+IFQL  -  ID  value  at  each  position 

♦J  -  event  number  on  previous  radial 

♦JA  -  address  variable  on  J 

+  IPB  -  number  of  peak  thresholds  in  event 

♦KA,+KAP,+KAM  -  address  variables  on  event  and  threshold 

+NP  -  number  of  contours,  this  threshold,  on  previous  radial 

♦ITHRM  -  code  for  eliminating  possible  cell 

♦MG  -  address  in  UP  (array)  for  measure  of  cell  significance 

+I.MT  -  limit  on  number  of  cells  to  be  carried  in  order  of  significance 

+JK1  -  address  variable  on  previous  radial 

+.1KI,  -  address  variable  on  current  radial 


+  -  INTF.GF.R*2 


1 


TRACK 


*VKM  -  cos  (elevation')  x  conversion  m  to  km 
*SAVKM  -  unit  area  convert  from  m~  to  krn^ 

+MA  -  address  in  l:.CL  (array)  of  measure  of  cell  significance 
+M  -  cell  counter 

+Ml,+M2,+M3,-*-M4,+M5,+M6,+M7  -  address  of  each  attribute  in  ECL  (array)  by  cell 
+K0FST,+NAN2  -  array  addressing  offset 

*VKMF,  -  convert  reflectivity  weighted  line  of  sight  distance  to  horirontal  diM. 
* FNSN , +NSCAN  -  scan  counters 
+KTL  -  time 

♦JPAY ,+IPAY  -  Julian  date 
♦JHR, +IHR 

+JMIN.+IMIN  start  time  of  first  scan  in  sequence 
+JSEC  ,  +  ISF.C 

♦NC  -  cell  number  (current  scan) 

+NCG,+NCBG,+NCB  -  NC  to  pass  through  common 

+NVMIN,+NYMX,+NCMX  -  number  of  significant  cells  detected  this  >can 
+IELSN  -  current  elevation  angle 
+IESNL  -  elevation  angle  last  scan 

ATRAK 

+NCEC  -  address  offset  on  NC 
+NVYC  -  address  offset  on  NV 
♦NC  -  cell  counter  for  ECL  (array) 

♦NV  -  cell  counter  for  VCL  (array) 

♦NCA  -  address  variable  on  NC 
+NYA  -  address  variable  on  NY 
♦IZ  -  cell  reflectivity 
*X  -  cell  position  east  [(-)west] 

*Y  -  cell  position  north  [(-)south] 

*H  -  cell  height 

+  1ZL  -  offsetted  reflectivity 

*Z  -  reflectivity 

*HI.  -  height  of  last  cell 

+  IZP  -  peak  reflectivity 

RTRAK 

♦NCEC  -  address  offset  on  NC 
+NWR  -  address  offset  on  NY 
♦NC  -  cell  counter  for  ECL  (array) 

+NV  -  cell  counter  for  VR  (array) 

+NYA  -  address  variable  on  NY 
♦NCA  -  address  variable  on  NC 

♦  -  INTEGER*  2 

*  -  REAL 


A- lb 


COMPAR 


+NCMX  -  number  of  cells  detected  this  scan 
+IM.+JM  -  array  limits 
+NC  -  cell  numbers  detected  this  scan 
♦NCEC  -  address  offset  on  NC 

♦  (.NCI  to  NC6)  -  address  variables  on  NC 

♦NV  -  cell  numbers  tracked  from  previous  scans 
♦NVVC  -  address  offset  on  NV  for  VCL  (array-) 

+NVYR  -  address  offset  on  NV  for  VR  (array-) 

+NLR  -  address  variable  on  NV 

+ATEST  -  estimate  of  cell  NVs  movement  from  last  scan  to  this  scan 
*nri:r  -  time  since  last  scan 

*DELX,  DELY  -  distance  on  X  f,  Y  coordinates  between  cell  NY  and  NC 
*l)ELW  -  a  measure  of  the  association  between  NV  and  NC  using  reflectivity, 
location,  area,  height 
+.10  -  overflow  of  0  and  11)  (arrays) 

+NSCAN  -  scan  number 
+  10  -  overflow  of  IllV  (arrays) 

*PX  -  DKLW  of  a  previous  association  with  this  NC  cell 
♦WT  -  NV  previously  associated  with  this  NC  cell 

♦NCT  -  NC  previously  associated  with  this  NV  cell 

*PX  -  DELW  of  a  previous  association  with  this  NV  cell 
♦NVMX  -  number  of  active  cell  tracks 
+NCR  -  NC  cell  to  test  NV  cells  against 

♦  M  B  -  NV  . 

+NCB  -  NC  1  asSOOKUed 

*HCT  -  update  height 


+  -  INTI  Cl  R*: 
*  -  RUAI. 


A- 10 


RF.SOt.V 


•NYT  -  NY  associated  with  Nl'T 

♦I  -  number  of  NY  associations  on  this  NCT 

♦JX  -  number  of  Nl'  cells  associated  with  1  NY  cells 

*NCT  -  NY  associated  with  NYT 

♦I  -  number  of  NY  associations  on  this  NYT 

♦JX  -  number  of  NY  cells  associated  with  1  NY  cells 

#NVT  -  NY  cell  that  NY  cell  is  associated  with  (highest  HI  IK 

♦.TV  -  number  of  cells  associated  to  this  NY  cell 

*kY  -  maximum  number  of  associations  to  one  NY  cell 

•IYS  -  NY  cell  association  on 

♦111  -  counts  NY  cells  checked 

•NYT  -  s:une  as  NAT  but  NY  on  NY 

+dt'  -  same  as  J\  lust  on  an  NY  cell 

*kl'  -  same  as  KY  but  on  an  NC  cell 

♦  Il'S  -  N('  cell  association  on 

♦  l.V  -  counts  NY  cells  checked 

♦NY  -  N\'  cell  that  Nl'  cell  is  associated  with 
♦KA,  +k\J ,  ♦k.-M  -  address  variables  on  kl' 

♦NC  -  NY  cell  that  NY  cell  is  associated  with 

♦IMSM  -  co*v.\kS  associations  on  NY  other  than  NO 

*P1:.UV  -  measure  of  compare 

*PWT  -  minimum  measure  of  compare 

+N\'  -  cell  with  PWT 

♦KNY  -  counts  cells  with  too  many  associations 
*P1T.T  -  test  on  cell  velocity 

♦NYT l  -  next  NY  cell  associated  with  Nl'  (in  order  of  P1T.K1 

*PWTl  -  minimum  measure  on  NYT1 

♦NY  1  -  cell  with  PKTl 

*Pl:l.K'l  -  test  on  NY1 

*Pl:.l,W:  -  test  on  NX 

*HT1-  -  update  height 

♦  -  TNThdl'R-: 


»  -  INT1Y1:R*4 


STRAK 


♦MV A  -  address  offset  on  NV 

♦  (NV1  to  NV41)  -  address  variables  on  NV  for  VCI.  (array) 

♦  NVR  -  address  offset  on  NV 

♦  (NR1  to  NRM  -  address  variables  on  NV  for  VR  (array) 

♦  NV'  -  cell  counter 

*VXT  -  velocity  east  [(-)west] 

*VYT  -  velocity  north  [(-) south] 

*DELTM  -  time  since  last  scan 

♦  IZV'AL  -  reflectivity 

♦  IDTC  -  percent  of  scans  cell  was  detected 
*VXC  -  sum  of  eastward  velocity  components 
*VYC  -  sum  of  northward  velocity  components 
+NFN  -  number  of  scans  processed 

♦NSN  -  number  of  velocity  values  summed  in  VXC  and  VYC 
*VX  -  average  eastward  velocity  of  all  cells 
*YY  -  average  northward  velocity  of  all  cells 
♦NVSCN  -  number  of  volume  scans 
+KTL  -  time  of  last  scan 


♦  -  INTEGER* 2 

*  -  REAL 
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LIST  OF  ARRAYS  -  THEIR  SI  IK  AND  CONTENT 


In  CRANE 


ANC(1028)  -  header  information  -  INTEGER* 4 
ZEE (1024)  -  raw  data  -  INTEGER*4 
IREF (1024)  -  reflectivity  -  INTE(IER*2 
IVEL(1024)  -  radial  velocity  -  INTEGER*! 

IVAR( 1024)  -  variance  -  INTEGER*! 

RE (1025)  -  decoded  reflectivity  -  INTEGER*! 
IVELL(!56)  -  decoded  velocity  -  INTEGER*! 
IDVEL(256)  -  tangential  shear  -  INTEGER*! 
IARY(91)  -  convert  rcfelctivitv  to  dBZ  -  REAL 


In  CONTOR 


ICL(44)*  -  start  position  of  event  -  INTEGER*! 

ICI(44)*  -  stop  position  of  event  -  INTEGER*! 

IOC(22)  -  number  of  peaks  in  each  event  -  INTEGER*! 

IPRNG(34)  -  location  of  peaks  -  INTEGER*! 

IC.21(22)  -  start  position  of  event  on  second  threshold  -  INTEGER*! 
rC 2 2(22)  -  stop  position  of  event  on  second  threshold  -  INTEGER*! 


*ai  ray  contains  indicated  parameter's)  on  the  current  radial,  offsetted 
from  the  same  parameter! si  on  the  previous  radial 


\- 1 9 


In  PEAKD 


T(80)  -  all  possible  thresholds  a  peak  may  have  -  INTEGF.R*2 

TC(1980)*  -  thresh  of  each  peak  -  INTEGER* 2 

IPTC(44)*  -  number  of  thresholds  in  each  event  -  INTEGER*2 

IPCNT(1980)*  -  contour  counter  -  INTEGER*2 

IPCI (54001*  -  start  range  of  contour  segment  -  INTEGER*2 

IPC2(5400)*  -  end  range  of  contour  segment  -  INTEGER* 2 

IPC3(5400)*  -  number  of  peaks  within  the  segment  -  INTEGER* 2 

TATR(1400)  -  temporary  attribute  array  -  stores  peak  attributes  until  a 
cell  is  detected  or  peak  discarded  -  REAL 
1ACT(701  -  overflow  (too  many  peaks)  -  INTEGER*2 

UP(6)  -  cell  attributes  -  REAL 

1  -  area 

2  -  reflectivity 

3  -  location  in  km  east  of  radar 

4  -  location  in  km  north  of  radar 

5  -  tangential  shear 

6  -  a  measure  of  relating  cell  significance  =  const  *  area  +  tan  shear 


*arrny  contains  indicated  parameter(s)  on  the  current  radial,  offsetted 
from  the  same  parameter(s)  on  the  previous  radial 
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PEAKD  f,  TRACK:  ECL  (7x16x21 


( 1  ,  N ,  N'AN  1 1 


I 

variable 

1 

area  surf 

rea  1 

■» 

dBZ 

real 

3 

east 

rea  1 

4 

north 

rea  1 

5 

range 

rea  1 

0 

height 

rea  1 

defined  in  TRACK 

t  ang  shear 

real 

ATRAK : 

YCl.  (23x321 

( 1  ,\’l 

I 

variable 

— 

1 

east  (XI 

real 

7> 

north  (Y1 

real 

3 

dBZ 

int*4 

1 

4 

area  surf 

rea  l 

set  in 

ATRAK 

lowest  elevation 

5 

t  into 

int*4 

def  in 

CRANE 

angle  onlv 

6 

height 

real 

'l 

*"* 

range 

real 

— 

t 

1 

S 

NT;  track  10 

int*4 

def  in 

ATRAK 

_ 

1 

9 

.7  # scans  each  dot. 

int  *4 

I 

10 

V  "* 

real 

i 

I 

11 

ZZX 

rea  1 

i 

12 

>:zy 

real 

updated  each 

13 

7Z(H-HL1 xArea 

rea  1 

scan 

14 

11  summit  (HT  lastl 

real 

15 

dBZ  summit 

int*4 

i 

10 

dBZ  peak 

int  *4 

l~ 

ll  peak 

real 

— 

— - 

18 

X 

rea  l 

def  in 

STRAK 

19 

Y 

rea  1 

| 

undated  each 

20 

time 

intM 

1 

vol  scan 

21 

Yel  x 

real 

i 

j 

■>  T 

Yel  v 

real 

1 

i 

23 

Z  tang  shear 

rea  1 

def  in 

ATRAK 

updated  each  sea 

STRAK*  TCI.  (21) 


( n 


1 

variable 

1 

t  imo 

int  *4 

•» 

X 

rea  1 

3 

Y 

rea  l 

•i 

♦« 

rea  l 

!> 

(ll-lll.)  xArea 

rea  1 

o 

(not  used)* 

7 

dB2  peak 

int  *4 

8 

II  peak 

rea  l 

0 

(not  used) 

10 

die. 

int  *4 

11 

area  surf 

rea  l 

12 

height 

rea  1 

13 

(not  used) 

14 

(not  used) 

IS 

die  summit 

int  *4 

16 

II  summit 

real 

17 

ABS  (NT)  I  RAK  III 

int  *4 

IS 

(nor  used) 

10 

(not  used) 

20 

AX/ AT 

rea  1 

21 

AY/  AT 

rea  1 

•array  locations  not  used  in  the  current  version  of  the  program 
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In  aiMPAII 

track  colls  from  previous  scan  (NY)  to  colls  in  current  scan  t Nil t 

lUCl(K')  -  points  to  NY  each  Nf  is  associate.!  with  -  1NTIV1!R*2 
i  .  o .  1  VC!  1  ( N(!)  *  NY 

ll'(l<>)  -  stores  the  "measure"  of  the  above  association  -  INTl.iil  U*2 
I  IkT  J*  t  It')  -  if  one  N'C  is  associated  with  more  than  one  NY  then  IVt'.'lNi'l  * 
lllVh(32)  -  points  to  NV  each  N\'  is  associated  with  lYl'llUK*’ 
tV(32)  -  stores  the  "measure"  of  the  above  association  INVIVIKV 
1 1 IV 2 C .> 2 >  -  if  one  NV  is  associated  with  more  than  one  NT  then  10Y2(NV1  i 

II! (32x101  -  if  IIKYCNin  *  0  -  INTIVIKV 

li!('!.ll  =  number  of  conflicts 
ll!(",2-10)  '  the  NY's  associated  with 

l! ( 32x0)  -  stores  the  measures  of  each  NC-NY  association  -  Uf.Al 

10(32x10)  -  same  as  the  ICiarrayl  for  conflicts  on  NY  -  1 NT1V1  11'* 2 

0(32x0)  -  stores  the  measures  of  each  NY -NY  association  KI'AI 


hi  RliSCH.V 

IV (32,")  -  ordered  U'. (array)  or  lO(arrav)  -  INTVVdRV 
V( 3SI)  -  stores  the  measures  of  the  tested  associations  - 


-  IN  VI  til  U 


-  IMIVIK 


KIM 


LIST  OF  COMMON  BLOCKS  AND  TIIFIR  ASSOCIATED  ROUT! NFS 


BLOCK 


RLOCK  NAME/SIZE* 

CRANE 

DATA 

CONTOR 

PEAKD 

TRACK 

ATRAK 

BTRAK 

COMPAR 

RESOLV 

STRAK  COMB  IN 

AZM 

i4.ro 

X 

X 

X 

X 

M2 

(M 

X 

X 

X 

X 

X 

C DRAYS 

(897) 

X 

X 

X 

CLST 

CS7) 

X 

X 

X 

X 

X 

X 

X 

CNT 

(31 

X 

X 

X 

X 

CNTRS 

(3) 

X 

CONI’ IC 

m 

X 

X 

CONST 

U  2) 

X 

X 

X 

X 

DATA  1 

c:< o 

X 

X 

X 

X 

X 

X 

X 

X 

data: 

(7371 

X 

X 

X 

X 

X 

X 

DATA  3 

(15)31 

X 

X 

X 

X 

X 

DVAI. 

(I) 

X 

X 

ECONST 

Cl 

X 

X 

X 

EXTRA 

( 35)7 .  SI 

X 

F I LTFR 

(31 

X 

X 

X 

X 

FIXED 

(-’.M 

X 

X 

X 

FLOS 

(3.51 

X 

X 

X 

X 

X 

INTI. 

(31 

X 

X 

X 

X 

X 

X 

X 

KTA 

Cl 

X 

X 

X 

X 

KTB 

HI 

X 

X 

X 

X 

MAPI’ A  R 

(4.51 

X 

NVLIS 

(41 

X 

X 

X 

X 

X 

X 

X 

X 

l’NI'RS 

(41 

X 

X 

X 

X 

X 

X 

X 

I’RSrOR 

(14051 

X 

X 

X 

t’WORK 

(43.51 

X 

X 

X 

REEL 

(  7"T0 

X 

X 

X 

X 

RSLV 

( 7 : .  5 1 

X 

X 

SWITCH 

(1047.51 

V 

X 

X 

thresh 

(  .51 

X 

X 

X 

IT.  IS 

(SI 

X 

X 

X 

X 

X 

X 

'I’M  AX 

(.51 

X 

X 

VPARM 

Cl 

X 

X 

X 

ZLOOK 

(5)1  .51 

X 

X 

X 

ZSTORF. 

( 10281 

X 

*sCo  is  number  of  3Z  hit  words 


AC4 

i 


Loop  on  Volume  Scan 


Calibration  Data 


i  on ) 


chant; 


START 


Event  Identification 

Subroutine  u'M'i'K 


loop  on  ThrathoM 


Loop  ivi  Saamant 


Loor  on  Paal,  wthm  f  v*nt 
U-ow*\t  1hr*»hok1) 


Loop  on  Ranga  w\Kn  t  wm 


Loop  on  ThrwhoM 


Loop  <M1  Snjtnant  naw  M,0na<  1hi*,t,old 
C  Radial 


Loop  on  TVathold  R  Papol 


loop  on  Saijmanti  B  Radial 


A»toC>ait  B 
to  C  Radial 


Oatact  Call  .t 
not  Updatad  th.«  Radial 
and  C  Radi#'  vaUm*  aia 
than  Cal  l  TVatppij 


Peak  Detection 

Suhr.'ut  ■ ’>«'  \KP 


Event  Association 


INPUT  AFTER 
EACH  SCAN 


07/19/79  16:31: 48 

♦♦♦LISTING  FOR  ERT1 :  TSEMflIN.  FTN 

♦N 

♦ASSM 

NORX3 

TSEMflIN  PROG 
♦FORT 

♦TITL  FILE  TSEMflIN  -  MAIN  PROGRAM  FOR  ETSE 
IMPLICIT  INTEGERS  <R-2) 

REAL  DB,  B I TVEL 

INTEGERS  RHO,  TFSO,  GRND,  1/  ZERO,  T,  TSV,  STOP,  ZTH,  TFS,  K,  BETA 
INTEGERS  INDCTR,  ANGC7),  DISPLA,  MMU,  SIGMA,  ELEVC2) 

INTEGERS  TIME,  RSAVE < 16),  SECOND,  NRCEAD 

INTEGERS  NEXT  I M,  I T I  ME 

REAL  PCTMIN,  RRAREA 

COMMON  /CALB/  DB, B I TVEL,  NRC 

COMMON  /MUS I G/PCTM I N,  MMU,  SIGMA 

COMMON  /READZ/NRCEAD,  OLD 

COMMON  /SECTOR/ INDCTR 

COMMON  /EXTRA/RHO, GRND, ZTH,  BETA, K, RRAREA C 12, 24), RH02, ZTH2, BETA2, 
1PCT2,  MMU2,  S I GMA2,  M I NUT2,  ANG2  <  6 ) ,  CHANGE 
EQU I VALENCE  <  D I SPLA,  ANG  < 1 ) ) , <  DEMAND, ANG  <  6 ) > ,  <  ELE V  < 1 ) ,  ANG  C  4  )  ) 
EQUIVALENCE  < ANGLE, ANG<2) ),  <OLDATA, ANG<7) ) 

DATA  DISPLA/4/,  NLEVEL/3/ 

REWIND  4 
INDCTR=0 
10000  CONTINUE 
♦ASSM 

FREZE 

COPV  SVC1. 

STM  0, RSAVE 

READ  SVC  1, READBLOK  READ  IN  DISPLAY,  ANGLES 

LH  R0, READBLOK+SVC1.  STfl  READ  IN  DEVICE  STATUS 
BNZ  IOERR  BRANCH  IF  NOT  ZERO 

LH  R0,  ANGLE  LOAD  IN  BEGINNING  ANGLE 

CHI  R0,  360  LESS  THAN  360? 

BNL  ALLD  NO,  FULL  SCAN 

LIS  R0,  1 

STH  P0,  INDCTR  INDCTR=1 

ALLD  LM  0, RSAVE 

♦FORT 

READ<7,  333)  I  OUT,  DBB,  MAXV,  MAXS,  SLOPE 
333  FORMAT <  13) 

BITVEL=MAXV 
DB=65  28 

READ <5,  110)RHO,  ZTH, BETA,  IPCTMN,  MMU, SIGMA,  MINUTE 
110  FORMAT  < 13) 

CHANGE*© 

PCTMIN* I PCTMN+0.  01 
SECOND*M I NUTE*60 
OLD=OLDATA 
50  CONTINUE 

CALL  I CLOCK C 2,  TIME) 

NEXTIM*TIM£+SECOND 
SVELE*RNG  <  5 ) 

ANG<5)*ANGC4>+5 
DO  51  N=l,  NLEVEL 
CALL  PPRDSC  < ANG ) 

CALL  RMRP 
RNG<5)-SVELE 

SI  CONTINUE  c-2 

CALL  PPRDSC < ANG) 


C  CfiLL  TSEPLT (DISPLA) 

I F  <  CHANGE  ) 55,  35,  52 
CHANGE «0 

52  IF (RH02.  NE.  (-1) >RH0=RH02 
IF<PCT2.  NE.  (-1>>IPCTMN=PCT2 
I F  (  MMU2.  NE  <-i>>NMU=MMU2 

IF (SIGMA2  NE.  (-1) )SIGMA=-SIGMA2 
IFC2TH2.  NE  ( -1 > )ZTH^ZTH2 
IF (BETA2.  NE.  ( -1 ) >BETA^BETA2 
IF<MINUT2.  NE.  (-1) )MINUTE=MINUT2 
DO  34  I«l,  6 

IF( ANG2( I >.  NE  (-1)  )ANG( I )=RNG2( I ) 
54  CONTINUE 

PCTMIN«IPCTMN*0.  01 
SECOND=60*M I NUTE 

53  IF (DEMAND- 1)56,  75,  30 

56  CALL  I CLOCK <2, I TIME) 


IF<  I  TIME.  GT.  TIME)GO  TO  58 
IF (NEXT I M.  LE  86400) GO  TO  50 
NEXT I M=NEXT I M-86400 


58  IF  (NEXT  IM.  LT  I T I  ME  ) G 

60  CONTINUE 

♦ASSM 

R0 

EQU 

0 

R1 

EQU 

1 

R2 

EQU 

2 

R3 

EQU 

3 

R4 

EQU 

4 

STM 

R0, RSAVE 

- -  . 

L 

Rl,  NEXT  I M 

ST 

Rl, SVCTIM 

SVC 

2. TIMBLK 

LM 

R0,  RSAVE 

♦FORT 

GO 

TO  50 

75  PAUSE  0 

GO 

TO  50 

♦RSSM 

ALIGN 

4 

IOERR 

EXBR 

Rl,  R0 

NHI 

Rl, X '00FF ' 

CHI 

Rl,  X  '  88  ' 

BE 

ALLD 

CHI 

Rl,  X'90 

BE 

ALLD 

SVC 

2, ERRCODE 

SVC 

2, ERRBLOK 

SVC 

2, PAUSE 

B 

READ 

ALIGN 

4 

READBLOK 

EQU 

* 

DB 

X'58' 

DB 

4 

— 

DC 

H'0' 

DC 

A(DISPLA) 

DC 

A(OLDATA+l> 

DSF 

3 

ALIGN 

4 

ERRCODE 

EQU 

>•< 

‘  -  -  -• 

DB 

0 

DB 

6 

DC 

H'0' 

DC 

A(ERCD) 

ALIGN 

4 

ERRBLOK 

EQU 

♦ 

STORE  TIME  OF  DRV 
TIME  OF  DRV  WAIT 
RETURN  TO  FORTRRN 


END  OF  FILE? 

VES,  FULL  SCAN,  DEFAULT  DISPLAY 
END  OF  MEDIUM 

VES, FULL  SCAN, DEFAULT  DISPLAY 

PRINT  OUT  ERR  MSG  TO  CONSOLE 

PAUSE 

TRY  AGAIN 


READ  BINARY  AND  WAIT 

LU  4 

STATUS 

BEGINNING  OF  BUFFER 
END  OF  BUFFER 


4  CHRS  OF  ASCII  HEX 


TO  ERROR  MSG 


r-.' 


30  CHK'S 


DB 

7 

DC 

H'30' 

30  CHRS 

ERCD 

DC 

C' 

DC 

C'  I/O  ERROR 

IN  RTN  TSEMRIN 

RLIGN 

4 

PRUSE 

EQU 

* 

DB 

0,  1 

RLIGN 

4 

TIM6LK 

EQU 

* 

DB 

0 

DB 

10 

CODE  10 

DB 

0,  0 

SVCTIN 

DC 

F  ■'O' 

TIME  OF 

♦FORT 


TIME  OF  DRV 


* 

07/19/79  17:07: 53 

♦♦♦LISTING  FOR  ERT1 : TSEDRTR.  FTN 

*N 

♦ASSM 

SCRRT 
SQUE2  3 
TSEDRTR  PROG 
*FORT 

$TITL  FILE  T'SEDATR  -  DATA  INPUT  SUBROUTINE  FOR  ETSE 
SUBROUTINE  PPRDSC<RNG) 

IMPLICIT  INTEGERS  <R-Z) 

I NTEGER^4  PPR I < 102S ) , PPR 1 2  < 102S  > ,  RSRVE < 16  > 

+, RAW,  DECOD <1028),  DECOD2<1028> 

RERL  DB, BITVEL 

INTEGER+2  I REF <1024), IVEL<1024>, IVELL<1024), 

+  IREF2<1024),  IVEL2<1024> 

INTEGERS  PPRRNG 

INTEGER+2  RNG<6),  RE <1025),  HR <258), TL1,  TL2,  RQURNT 
EQU I VRLENCE  <  PPRRNG, PPR I <  3 ) ) ,  <  PP,  PPR I <  2  > ) 
EQUIVALENCE  <DEC0D<5),  IREF<1)>,  <DEC0D<517),  IVEL<1)), 
+<DEC0D2<5),  IREF2<1) ),  <DEC0D2<517),  IVEL2<1)) 

COMMON  /REFL/  RE,  HR,  NCL,  NID,  NIDP,  INCL,  I  MX, 

+  IMN,  TL1, TL2, RQURNT,  IDVEL<258) 

COMMON  /CAL/DB,  BITVEL,  NRC 
COMMON  /SECTOR/ I NDCTR 
COMMON  /ZSTORE/PPR I 
'  '  COMMON  /RUNSUM/PPRI2 

CALL  C0NMSG<6, "PPRDSC") 

2  REWIND  9 

NRZ=0 

ELEV1=RNG<4)+11.  37778 
ELEV2=RNG  <  5 ) *11.  37778 
IF<INDCTR.  EQ.  2)  I  NDCTR  =  0 
IF  < I NDCTR  EQ  0)GO  TO  1 
CW  =  ANG<3)  -  RNG<2) 

IF  <CW  .  GT.  180)  CW  =  CW  -  360 

IF  <CW  .  LT  <-180)>  CW  =  360  +  CW 

IF  <CW  .  LT.  0)  CW  =  0 

IF  <CW  .  GT.  0)  CW  =  15 

BGNR  =  11.  37778  ♦  RNG<2> 

ENDR  =  11.  37778  ♦  RNG  <  3 ) 

BGNA=MOD  <  BGNR,  4096) 

ENDfl=MOD  <  ENDR,  4096) 

IF < ENDR.  GT.  40S6)ENDA=0 

STPFLG  =  0 

NITG=4 

NRC=768 

IF  <CW  .  GT.  0  .  AND.  ENDR  .  LT.  BGNR)  STPFLG  =  1 

IF  <CW  .  EQ.  0  .  AND.  ENDR  .  GT.  BGNR)  STPFLG  =  1 

IF<BGNR  GT.  10  .  AND  BGNR.  LT.  4086 )GO  TO  i 


BGNA=0 

1 

STPFLG=2 

CONTINUE 

10000 

CONTINUE 

♦ASSM 

FREZE 

CROSS 

COPV 

SVC1 

R0 

EQU 

0 

R1 

EQU 

1 

R2 

EQU 

2 

R3 

EQU 

3 

< 

STM 

0,  RSRVE 

L 

R4,  WAITRERD+SVC1 

SRD 

RIS 

R4,  15 

ST 

R4,  WRITRERD+SVC1  ERD 

BflL 

R13,  WREAD 

LIS 

Rl,  4 

L 

R0, PPRKR1) 

NHI 

R0,  3 

RIS 

R0,  1 

SLLS 

R0,  8 

RIS 

R0,  4 

SLLS 

R0,  2 

SIS 

R0,  1 

L 

R4,  PPRIBLK+SVC1. 

SRD  GET  BEGINNING  ADDRESS 

RR 

R4,  R0 

COMPUTE  END  ADDRESS 

ST 

R4, PPRIBLK+SVC1 

ERD  STORE  IN  ERD 

L 

R4,  Wfl I TRERD+SVC1.  SRD 

RR 

R4,  R0 

ST 

R4,  WAITREAD+SVC1 

ERD  SAME  NUMBER 

L 

R4, OUTBLK+SVC1.  SRD 

RHI 

R4, 1028*4-1 

ST 

R4, OUTBLK+SVC1  ERD 

L 

R4, PPRIBLK2+SVC1 

SRD  GET  NEXT  BEGINNING  RDDR 

RR 

R4,  R0 

COMPUTE  END  RDDRESS 

ST 

R4, PPRIBLK2+SVC1 

ERD  STORE  IN  SVC  BLOCK 

€ 

L 

R4, 0UTBLK2+SVC1 

SRD 

RHI 

R4,  1028*4-1 

ST 

R4, 0UTBLK2+SVC1 

ERD 

LDRI 

R15, PPRI 

** 

LDRI 

R14, PPRI 2 

LCS 

R6,  1 

SET  R6  TO  -1  FOR  COUNTER 

DETECT 

BRL 

R13,  WRERD 

READ  IN  RN  AZIMUTH 

LH 

Rl, 8<R15> 

LORD  IN  AZIMUTH  DATA 

NHI 

Rl,  X'FFF' 

RND  OUT  UNWANTED  BITS 

BflL 

R13,  WRERD 

GET  ANOTHER  AZIMUTH 

LH 

R2, 8<R15> 

GET  AZIMUTH  DRTR 

NHI 

R2, X'FFF' 

AND  OUT  UNWANTED  BITS 

CR 

Rl,  R2 

COMPARE  TWO  AZIMUTHS 

BL 

CWISE 

IF  RKR2,  RADAR  IS  GOING  CWISE 

LIS 

R3,  0 

DIRECTION  FLAG 

B 

WHflTIZIT 

CONTINUE 

CWISE 

LIS 

R3,  15 

DIR  15LAG  =  CW 

WHATIZIT 

LH 

R4,  INDCTR 

SECTOR  SCAN  OR  FULL  CIRCLE? 

BZ 

EDETECT 

FULL  CIRCLE 

CH 

R3,  CW 

IS  DIRECTION  OF  ROTATION  CORRECT 

BNE 

DETECT 

WRONG  DIRECTION,  WAIT 

OR 

R3,  R3 

WHICH  DIRECTION  IS  IT? 

BZ 

CCW 

COUNTER  CLOCKWIZE 

LH 

R5, STPFLG 

THI 

RS,  2 

CASE  2? 

BNZ 

CWCASE2 

CW  CASE  2 

< 

CMF*1 

CH 

R2, BGNR 

ANGLE  <  BGNA? 

BL 

WRIT 

VES,  GET  RERDV 

B 

DETECT 

NO,  TRY  RGRIN 

# 

CWCASE2 

CHI 

R2,  X '  800  ' 

ANGLE  >  180? 

BL 

CMP1 

NO,  ALL  OK 

SHI 

R2,  X'1000 ' 

VES,  SUBTRACT  360 

0 

B 

CMP1 

CCW 

LH 

R*>,  STPFLG 

(■-<> 

BNZ 

CCWCASE2 

CMP2 

CH 

R2,  BGNA 

BL 

DETECT 

WAIT 

BAL 

R13,  WREAD 

LR 

Rl,  R2 

LH 

R2,  8<R15> 

NHI 

R2,  X'FFF' 

CR 

Rl,  R2 

BL 

CW3 

LIS 

R3,  0 

CMP6 

CH 

R3,  CW 

BNE 

DETECT 

OR 

R3,  R3 

BZ 

CCW2 

THI 

R5,  2 

BNZ 

WCWCASE2 

CMP3 

CH 

R2,  BGNA 

BL 

WAIT 

READ1 

LH 

R2, 12<R15> 

NHI 

R2,  X'FFF' 

CHI 

R2,  681 

BP 

ZERO 

ELDET 

CH 

R2, ELEV1 

BM 

DETECT 

CH 

R2,  ELEV2 

BP 

DETECT 

AIS 

R6,  1 

BNP 

DETECT 

READ 

BAL 

R13,  GOREAD 

LH 

R2,  PPRANG 

NHI 

R2,  X'FFF' 

OR 

R3,  R3 

BZ 

CCWl 

THI 

R5, 1 

BN2 

CWCASE1 

CH 

R2,  ENDA 

BNL 

DONE 

B 

READ 

CW3 

LIS 

R3,  15 

B 

CMP6 

ZERO 

LIS 

R2,  0 

B 

ELDET 

CWCASE1 

CH 

R2,  BGNA 

BNL 

READ 

CH 

R2,  ENDA 

BNL 

DONE 

B 

READ 

CCWl 

THI 

R5,  1 

BNZ 

CCWCASE1 

CH 

R2,  ENDA 

BL 

DONE 

B 

READ 

CCWCASE1 

CH 

R2,  BGNA 

BL  ~ 

READ 

CH 

R2,  ENDA 

BL 

DONE 

B 

READ 

CCU2 

THI 

R5,  2 

BNZ 

WWCASE2 

CMP4 - 

CH  ' 

R2,  BGNA 

BNL 

WAIT 

B 

RERD1 

WWCASE2 

CHI 

R2,  X  ' 800  ' 

BNL 

CPIP4 

AH  I 

R2, X ' 1000 ' 

ANGLE  >  BGNA? 

NO,  TRV  AGAIN 
VES,  GET  READV 

GET  NEXT  AZIMUTH 
AND  OUT  UNWANTED  BITS 


CASE  2? 

VES,  BRANCH 
ANGLE  >  BGNA? 

NO,  KEEP  WAITING 


INCREMENT  COUNTER 
DO  IT  TWICE  TO  BE  SURE! 
VES,  START  READING 


BRANCH  IF  COUNTER  CLOCKWISE 
CASE  1? 

VES,  BRRNCH 
ANGLE  >  ENDA? 

VES,  ALL  DONE 
NO,  KEEP  READING 


59.  5  DEG  =  0 

ANGLE  <  BGNA? 

VES,  KEEP  READING 
NO,  ARE  WE  DONE  VET? 

VES 

NO,  KEEP  READING 
CHECK  FOR  CASE  1 

ANGLE  <  ENDA? 

VES,  ALL  DONE 
NO,  CONTINUE 
ANGLE  >  BGNA? 

NO,  KEEP  READING 

VES,  CHECK  FOR  FINISHED 


OK 

NOT  VET 


C-7 


C'  CCWCRSE2  CHI 
BNL 
AHI 

r  b 

WCWCRSE2  CHI 
®L 

c  SHI 

x  B 

EDETECT  BRL 

f  LH 

NHI 
CHI 

f  BP 

CPEV  CH 

BL 

C  CH 

BP 
LH 
NHI 
CHI 
BL 

READ2  BRL 

LH 
NHI 

NXT  OR 

BZ 
CR 
BL 
B 

2ER01  LIS 

B 

C  FUDGE  LIS 

B 

NXT2  BRL 

LH 
NHI 
OR  R3 
BZ 
CR 
BP 
B 

CCL0K4  CR 

ft  BP 

B 

CCW4  CR 

»  BL 

B 

OONE  LH 

i  NHI 

STH 
SVC 

I  LH 

BNZ 
LM 

|  #FORT 

RETURN 


R2, X'800' 

CMP2 

R2, X'lOOO' 

CMP  2 

R2»  X  " 800  " 

CMP3 

R2, X'1000' 

CMP3 

R13- WRERD 
R2> 12<R13> 

R 2, X'FFF' 

F2,  681 
TER01 
R2-  ELEV1 
EDETECT 
R 2-  ELEV2 
EDETECT 
R2,  PPRRNG 
R2,  X'FFF' 

R2,  6 
FUDGE 

R13, GORERD 
R4, PPRRNG 
R4, X'FFF' 

R3,  R3 
CCL0K4 
R4,  R2 
NXT2 
RERD2 
R2,  0 
CPEV 
R2,  6 
RERD2 
R13,  GORERD 
R4, PPRRNG 
R4,  X'FFF' 

.  R3 
CCW4 
R4,  R2 
DONE 
NXT2 
R4,  R2 
NXT2 
RERD2 
R4,  R2 
DONE 
NXT2 
R15,  PPRI 
R13, X'OFFF' 

R15,  PPRI 
1.  OUTBLK 

R0,  OUTBLK ♦SS'Cl  S’ 1 

ERROR 

0/  RSRVE 


RNGLE  <  180? 

NO/  RLL  OK 

VES<  ROD  360 

NOW  CHECK  FOR  BEGINNING 

RNGLE  >  180? 

NO* RLL  OK 

VES,  SUBTRRCK  360 

RERD  IN  NEW  AZIMUTH 
RERD  IN  ELEVATION 
AND  OUT  UNWRNTED  BITS 


WITHIN  RANGE? 
VES/  GET  AZIMUTH 

START  READING 

COUNTERCLOCKWISE 


IF  O.  MAKE  IT  6 

GET  NEW  AZIMUTH 
AND  OUT  UNWRNTED  BITS 


ANG  >  ENDR? 

VES/  FINISHED 
NO, KEEP  READING 


♦RSSM 


1, WA I TREAD 

f^ITRER04SVC1  st*  rerd  Returned  STATUS 
J^°R  IF  NOT  ZERO,  ERROR 

0  RSRVELK2  RERD  IN  °NE  **1"UTH 

c-.s 


Mil 


10  DECOD< I >“PPRI < I > 

JSIZ-NRC/NITG 

IMX»JSI2+1 

K«5 

DO  101  I»i, JSIZ 

REF-0 

VEL-0 

DO  20  J«l,  NITG 
RAW*PPRI<K> 

K«K+1 

♦ASSM 

ST  "  R0/RSAVE 

L  RO,  RAW  GET  PACKED  SOURCE  WORD 

EXHR  R0,  R0  SHIFT  POWER  TO  LOW  HALF 

NHI  R0, X'lFF'  NINE  BITS 

AHM  R0, REF  INTEGRAE  REFL 

L  R0,  RAW  GET  LOW  HALF 

SRHA  R0/  8  SHIFT  DOPPLER  TO  LOW  ORDER 

AHM  R0,  VEL  INTEGRATE  DOPPLER 

L  R0,  RSAVE  RESTORE  REGISTER 

♦F  ORT 

20  CONTINUE 

IFCREF.  LE.  0>GO  TO  111 
REF*R£F*.  0390625-DB 
IFCREF.  LT.  -39  >  REF=REF+100 
1REF  < I )=REF 
VEL“V£L*6  I  TVEL/-128. 

I VEL  < I > “VEL~ I VELL  < I > 

IVELLCD-VEL 
GO  TO  101 
111  I REF  < I ) “0 
1VELL< I >“0 
IVELa>“0 
101  CONTINUE 
♦ASSM 

STM  O, RSAVE 

SVC  1,  OUTBLK  OUTPUT  LAST  AZIMUTH 

SVC  1/  WAITBLK  WAIT  FOR  FIRST  READ  TO  FINISH 

SVC  If PFRIBLK  READ  NEXT  AZIMUTH 

LM  0.  RSAVE 

♦FORT 

DO  40  I«l,  4 
40  DEC0D2C I >“PPRI2< I > 

JS I Z“NRC/N I TG 

IMX*JSIZ+1 

K«5 

DO  201  1-1. JSIZ 

REF“0 

VEL«0 

DO  40O  J-l,  NITG 
RAW-PPRI2(K> 

K»K+1 

♦ASSM 


ST 

L 

R0,  RSAVE 

R0.  RAW 

GET  PACKED  SOURCE  WORD 

EXHR 

R 0,  R0 

SHIFT  POWER  TO  LOW  HALF 

NHI 

R0, X'lFF  ' 

NINE  BITS 

AHM 

R0»  REF 

INTEGRATE  REFL 

L 

R0,  RAW 

GET  LOW  HALF 

SRHA 

R8,  8 

SHIFT  DOPPLER  TO  LOW  ORDER 

AHM 

R0.  VEL 

INTEGRATE  DOPPLER 

L 

R0,  RSAVE 

RESTORE  REGISTER 

♦FORT 

400  CONTINUE 

IFCREF.  LE.  0>  GO  TO  411 


.  •  \«ni  WV>4>»*  L-'U 

IF  <REF.  LT.  -33>  REF  “REF +100 
IREF2< I >“RFF 
VEL“VEL*B I TVEL/12S. 

IVEL2C I >“VEL-IVELL< I ) 

I VELLC  I  >**VEL 
GO  TO  201 
411  IREF2C I >“0 
IVEL2C I >“0 
IVELL< I >“0 
201  CONTINUE 
♦RSSM 

STM  0,  RSAVE 

SVC  1, 0UTBLK2  OUTPUT  LAST  AZIMUTH 

SVC  1/ WRITBLK  WRIT  FOR  READ  TO  FINISH 

LH  R0, PPR I BLK2+SVC1.  STR  RERD  STRTUS 

BNZ  ERROR  IF  NOT  ZERO,  ERROR 

LH  R0, PPRIBLK+SVC1  STA 

BNZ  ERROR 

LH  R0,  0UTBLK+SVC1  STR 

BNZ  ERROR 

LH  R0,  0UTBLK2+SVC1  STR 


BNZ 

ERROR 

LH 

R6, NRZ 

RIS 

R 6,  2 

STH 

R6>  NRZ 

CHI 

R6,  446 

BLR 

F13 

LIS 

P6,  2 

STH 

R6, INDCTR 

B 

DONE 

ERROR 

SVC 

2, ERRCODE 

SVC 

2. ERRBLOK 

SVC 

2, PAUSE 

LM 

0, RSAVE 

++**P2  IS  FORTRAN  STMT  NO 

B 

*P2 

RLIGN 

4 

PAUSE 

EQU 

* 

DB 

0,  1 

RLIGN 

4 

ERRCODE 

EQU 

+ 

OPT 

DB 

0 

DB 

6 

DC 

H'0' 

DC 

ACERCD} 

RLIGN 

4 

ERRBLOK 

EQU 

+ 

DB 

0,  7 

DC 

H'15' 

ERCD 

DC 

C' 

DC 

CM/O  ERROR 

RLIGN 

4 

WRITBLK 

EQU 

* 

DB 

X'OS' 

DB 

10 

DB 

0,  0 

DSF 

5  • 

RLIGN 

4 

WR I TREAD 

EQU 

♦ 

- ~~ 

DB 

DB 

10 

DB 

0,0 

DC 

A<PPRI> 

DC 

ACPPRO 

DSF 

3 

RDD  TWO  TO  RIIMUTH  CTR 

TOO  MRNV  AZIMUTHS? 

NO,  KEEP  GOING 
VES,  INDCTR“2 

QUIT 

CONVERT  ERROR  CODE 
OUTPUT  MSG  TO  CONSOLE 
TASK  PAUSED 

RESTORE  FORTRAN  REGISTERS 
START  OVER 


PAUSE 


CODE  6 

DESTINATION 

PRINT  CONSOLE  MSG 
CODE  7 

PRINT  15  CHRS 
ERROR  CODE 


WAIT  ONLV 
LU  10 


READ  AND  WRIT 
LU  10 


C-10 


ALIGN 

^  PPRIBLK2  EQU 
DB 
DB 

C  DB 

DC 

JT  DC- 

C  DSF 

ALIGN 

0UTBLK2  EQU 
DB 
DB 
DB 
DC 
DC 
DSF 
ALIGN 

PPRIBLK  EQU 

DB 
DB 
DB 
DC 
DC 
DSF 
ALIGN 

OUTBLK  EQU 


4 

* 

X'51' 

18 

0,  8 

A<PPRI2> 

A<PPRI2> 

2 

4 

* 

X'31' 

9 

0,0 

ACDEC0D2) 

A<DEC0D2> 

3 

4 
* 

X'51" 

10 

0,  0 

A  <  PPR I ) 
A<PPRI> 


READ 
LU  10 


LU  9 


READ 
LU  10 


OUTST 


4F0RT 


X'31' 

9 

0,  0 

A<DECOD> 
A < DECOD) 
3 


WRITE 
LU  9 


RETURN 

END 


C- 1 1 
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♦♦♦LISTING  FOR  ERT1 : REALTM  FTN 

*N 

♦RSSM 

ERTRMAP  PROG 
♦FORT 

♦TITL  FILE  TSERMRP  -  -  PRINT  OUT  DATR  F I ELDS-CHRNGED  FOR  ERT  RERD  BV  CLB 
SUBROUTINE  RMRP 
IMPLICIT  INTEGERS  <I-N> 

INTEGERS  NRC,  VAR,  STORE <10),  PLACE,  OLDRTA 
INTEGER+4  RNC<1028>,  I DTI ME,  I2F,  ITZ,  I  ZB,  I2S 
+, SRC  <  1028  >,  DST (1028) 

INTEGER*2  RH02,  2TH2,  BETR2,  PCT2,  SIGMR2, ANG2, CHRNGE 
INTEGERS  TWENTY, ELEVEN,  DRV,  HOUR,  MINUTE,  SECOND,  TP,  ELEVRT ,  R2 
INTEGER*2  RZIM,  TC,  TR,  T 

INTEGERS  V.  I,  THETR1,  NRC1,  II,  THETR,  RHO,  STOP,  GRND 

INTEGERS  MERN,  POWER,  SIGMR,  TP2,  TP3,  SEGNO,  Q,  J,  2TH,  MMU 

INTEGERS  BEGIN,  SUM1,  SUM2,  JMIN,  M,  K,  L,  I 

INTEGERS  VI,  ZERO,  TWO,  BETA 

INTEGERS  RSAVE<16>, R1SRV 

INTEGERS  JREFC1024),  JDVEL<1024> 

INTEGER+2  RE,  HR,  TL1,  TL2,  RQURNT 
REAL  RNN,  RRRRER 

REAL  PCTMIN, RZT, ELEVHR,  BGNR,  R2CHK,  ENDR 
INTEGERS  FLAG,  IQ,  IB,  IE,  I M 
INTEGERS  DBB,  SLOPE 
INTEGER+4  NRCERD 

COMMON  /TLIS/  T,  Tfi,  JDRV,  JHR,  JMIN,  JSEC,  IDRV,  I  HR,  IMIN,  ISEC 
COMMON  /CflLR/  SRC, DST,  I REF C 1024 > 

COMMON  /SWITCH/  IC1<44>,  IC2C44>,  TC<1980>,  IPC1<5400>, 

♦  IPC2C5400),  IPC3<5400>,  IPCNT<1980>,  IPTCC44>, 

+  NEMB, NEMC,  NRC 

C0MM0N/R22/S I NR, COSR,  DELTR2,  ISCRNF, NEL, RI, SR 
C0MM0N/R2M/R2T, R2LRST ,  RZSTRR,  NR,  ELEVRR 
COMMON/REEL/RE  < 1025  > ,  HR<258>,  NCL,  NID,  NIDP,  INCL 
X,  I  MX,  IMN,  TL1,  TL2,  RQURNT,  IDVEL<25S> 

COMMON  /PNTRS/  NCMX,  NVMIN,  NUMX,  IELSN, NSCRN,  IESNL, NVSCN,  NT 
COMMON  / 1 NTL/  MHSN,  MNSN,  HM,  FNSN 

COMMON/EXTRfl/RHO, GRND,  2TH,  BETR,  K,  RRRRERC12, 24 >, RH02,  2TH2. BETR2, 

1PCT2,  MMU2,  SIGMR2,  MINUT2,  RNG2<6>,  CHRNGE 
COMMON/ZST  ORE/RNC 
COMMON  /ZIOOK/  I20FF, 2RRV<91> 

COMMON  /ECONST/  EARTH, VMK 

COMMON  /MRPPRR/  DRV,  HOUR,  MINUTE,  SECOND, DBB, MRXV, MRXS,  IOUT,  SLOPE 

COMMON  /CNT/  COSPHI, SINPHI,  C0SPH2 

COMMON  /DRTR1/  ECL<224>,  NOFST,  KOFST,  ICLRD,  NRN1 

COMMON  /DRTR2/  VCL < 736 >,  MXVC,  NVC1 

COMMON  /DATA!/  VR  < 192  > , MXVR,  NVR1 

COMMON  /NVLIS/  WARM,  NCRRM,  NVO,  ICO,  10,  JO,  JVR,  KTL 
COMMON  /FILTER/  TRTRMN, RRERMN,  DR2M 
EQUIVALENCE <RNC<3>,  JREF<2>>,  <RNC<317>,  JDVEL<2>> 

DRTR  TWENTV/18/ 

DATA  TWO/2/ 

DATA  PCTMIN/.  05/, MMU/0/,  SIGMfl/0/,  OLDRTR/0/ 

DATA  ZERO/0/ 

C 

C  *******  INITIALIZE  ARRAY 
C 

T-0 

AZLAST—999 
00  901  J-l,  MXVC 

901  VCL<J>«0  i-i: 


UU  J-i<  I'l.-AVIV 

902  VR<J>=0. 

RE<1>=0 

ICC=0 

C  CALL  C0NMSGC7,  'TSERMAPO 

TL1*TL1/RQUANT 
TL2=TL2/RQUANT 
REMIND  4 
REMIND  6 
DO  3  IX«1,  91 
I I=IX-IZOFF 
ZX=FLOAT< 1 1 >/10. 

3  ZARVC IX) =10.  **ZX 

00001  CONTINUE 

00750  FORMAT <1X>  'ENTER  PRF,  0=768,  1=922,  2=1075,  3=12290 
RCKM=.  075 
FLAG=0 

333  FORMAT  < 1 3  > 

READ<7,  333>0LDATA 

01002  FORMRTC/V"  AZIM  ELEV  RAN  ',8C'  REF  VEL  STDO> 
REMIND  9 
K=0 

10000  CONTINUE 

♦ASSM 


R0 

EQU 

0 

R1 

EQU 

1 

R2 

EQU 

2 

R3 

EQU 

3 

R4 

EQU 

4 

R5 

EQU 

5 

R6 

EQU 

6 

R7 

EQU 

7 

R8 

EQU 

8 

R9 

EQU 

9 

R10 

EQU 

10 

Rll 

EQU 

11 

R12" 

EQU 

12 

R13 

EQU 

13 

R14 

EQU 

14 

R15 

EQU 

15 

STM 

R0, RSAVE 

FREZE 

COPY 

SVC1. 

L 

R3,  ANCBLK+SVC1.  SAD 

ANCBLK  START 

AIS 

R3»  15 

READ  ANCILLARY  ONLY 

ST 

R3,  ANCBLK+SVC1.  EAD 

SVC 

1, ANCBLK 

THE  RERD 

LH 

R0, ANCBLK+SVC1.  STA 

BNZ 

ERROR 

LIS 

Rl,  0 

STH 

Rl, STOP 

LIS 

R5,  0 

MORKING  REGISTER 

LIS 

R9,  10 

MULTIPLICAND 

LHI 

Rll,  100 

MULTIPLICAND 

L 

R0, ANC<R1V 

DATA 

STBR 

R0,  R5 

HOUR 

LIS 

R10,  15 

MASK 

NR 

R10. R5 

1  HOUR 

SRLS 

R5,  4 

10  HOUR 

NHI 

R5,  3 

MHR 

R5»  R9  - 

C*10> 

AHR 

R5,  R10 

TOTAL 

STH 

R5,  HOUR 

HOURS 

EXBR 

R0,  R0 

STBR 

R0,  R5 

MINUTES 

LIS 

R10,  15 

MASK  C- 13 

< 


tw  r\xc>,  KO 

SRLS  R5/  4 
NHI  RS/  7 
MHR  R5/  R9 
RHR  R5/ R10 

STH  R5, MINUTE 

EXHR  RO, R® 

STBR  R®, R5 

LIS  Rl®, 15 

NR  Rl®/  R5 

SRi  S  R5/  4 
MHR  R5,  R9 
HHR  R5,  Rl® 

SRLS  R0>  8 
NHI  R0, X  F' 

MHR  R0. Rll 

RHR  R5, R0 

STH  R5/  DHV 

RIS  Rl. 4 

L  R0»  RNC  <  Rl > 

LHI  Rl®/  256 
STH  RIO. NRCl 

LHI  10/  1028*4-1 
R  10/  NRCBLK+SVC1. 

ST  10/  NRCBLK+SVC1 
ST  10, NRCERD 

NHI  R3/  X'C®' 

SRLS  R5,  6 
STH  R5,  TP 
EXHR  R0,  R0 
STBR  R0, R5 

LIS  Rif*.  15 

NR  RIO,  R5 

SRLS  R5, 4 

NHI  R5,  7 
MHR  R9,  R5 
RHR  R9, R10 

STH  R9, SECOND 

RIS  Rl,  8 

LHL  R0, RNC  <  Rl > 

NHI  F  0,  X'FFF' 

STH  F  0, ELEVRT 

LM  R0, RSRVE 

♦FORT 

NRC=NRC1 
GRND- FCTM I N*NRC 
REWIND  9 
17  ELEVEN*© 

TR*T 

DR=  <  2**TP ) *RCKM 
K-X+l 
NRl 1=NRC 
10001  CONTINUE 
♦RSSM 

STM  R®,  RSRVE 

SVC  1. NRCBLK 

LH  R®,  NRCBLK+SVC1 

BNZ  ERROR 

LIS  Rl, 0 

LH  R®, RNC  <  Rl > 

BM  MINUS 

LIS  Rl,  1 

STH  Rl,  STOP 

MINUS  LIS  Rl,  B 


■±  mwtt 
10  MINUTE 
MASK 
<*10> 

TOTRL  MINUTES 


DRVS 
MRSK 
1  DRV 
1®  DRV 
<*10> 

10®  DRVS 
MRSK 
<*100> 
TOTRL  DRVS 
DRVS 

NEW  DRTR 


STORE 

SRD 

EflD  STORE  END  RDDRESS  FOR  RERD 
STORE  RWfiV  FOR  SECOND  READ 
MRSK  TP 

STORE < UNFIXED) 

SECONDS 
MASK 
1  SECOND 

MRSK  10  SECONDS 
<*10> 

TOTRL  SECONDS 

STORE 

MORE  DRTR 

ELEVATION 

ANGLE 

STORE  ANGLE 


RERD  IN  DRTR 
STR  LORD  IN  STATUS 

IF  NOT  ZERO,  I/O  ERROR 


IF  STOP*® 
STORE  1 

GET  RZIMUTH 


LH  R®,  ANCCRl)  ANGLE 

NHI  R®,  X'FFF'  t:*14  MRSK 


J  I  ORE. 


» i  pi  *.w»  rw. 

AIS  Rl,  4 

LHL  R0>  RNC<R1>  ELEVATION 

NHI  R0,  X'FFF'  ANGLE 

STH  R0, ELEVAT  STORE  IT 

♦FORT 

DO  101  l  =  IMN,  I MX 
RE  < I >  * JREF  < I > 

IDVELC I >=JDVEL< I > 

101  CONTINUE 
00500  CONTINUE 
ICC=ICC+1 

ELEVAA=ELEVAT/11.  37776 
I ELSN= 1 F I X  <  ELEVAA > 

T* <  <DAV*24+H0UR  >  *60+MINUTE  >  +60+SECOND 
AZT=A2/11.  37778 
FLAG-FLAG+1 
A=AZT*.  01743 
RZCHK=R2T 

IF<K.  NE.  1)  GO  TO  105 

BGNA=A2T 

ENDA= A2T +359. 

105  CONTINUE 

IFCK.  LT.  180 >  GO  TO  106 
IF<A2T.  GT.  180  >  GO  TO  106 
AZCHK*=A2T +359. 

106  CONTINUE 

IF < A2CHK.  GT.  ENDA>  K=1 
SINA=SIN<A> 

COSA*COSCA> 

DELTR2=0.  0191987 
NA=K 

R£<258>=0 

IF  <NA.  EQ  1)NAC=1 

NAC=NAC+1- 

IFCNAC  GT.  1>NAC=0 

IF<NA.  NE.  1.  OR.  A2LAST.  LT  -990  >  GO  TO  224 

IDRV=DAV 

IHR=HOUR 

IMIN=MINUTE 

ISEC=SECOND 

PH I =ELEVAR*DRZM 

COSPH I -COS  <  PH I > 

SINPHI=SIN<PHI> 

COSPH2»COSPN I *COSPH I *EARTH 
CALL  TRACK 
224  CALL  CONTOR 
RZLAST=A2T 

IF<STOP.  NE.  1)  GO  TO  17 
800  CALL  STRAK 
STOP 


♦ASSM 

ERROR 

ALIGN 

SVC 

4 

2. ERRCODE 

DECODE 

ERROR  BITS 

SVC 

2, ERRBLOK 

OUTPUT 

ERROR  MSG  TO  CONSOLE 

SVC 

LM 

2.  PAUSE 

0, RSAVE 

RETURN 

TO  FORTRAN 

♦FORT 

IF<K.  EQ.  0>  GO  TO  10000 
IF<K.  EQ.  1)  GO  TO  10001 

-  GO  TO  1 

10003  CONTINUE 
♦ASSM 

ALIGN  4 
ERRCOPE  EQU  * 

DB  0 


('-!? 


DB 

6 

DC 

H'Q' 

DC 

A<ERCD> 

ALIGN 

4 

ERRBLOK 

EQU 

* 

DB 

0 

.  -  • 

DB 

7 

DC 

H'28' 

ERCD 

DC 

C' 

DC 

C'l/O  ERROR 

IN  RTN  RANGE-' 

ALIGN 

4 

PAUSE 

EQU 

♦ 

DB 

0, 1 

ALIGN 

4 

ANCBLK 

DB 

X'59' 

READ  BLOCK  FOR  LITTLE  I 

DB 

9 

LU 

DB 

0,  0 

DC 

A<ANC> 

START  ADDRESS 

~ 

DC 

ACANC) 

END  ADDRESS 

DSF 

3 

NRCBLK 

DB 

X'59' 

READ  BLOCK  FOR  BIG  READ 

DB 

9 

LU 

DB 

3,  0 

STATUS 

DC 

A< ANC> 

START  ADDRESS 

DC 

A<ANC> 

END  ADDRESS 

DSF 

3 

♦FORT 

END 


i 


aiti-nih  x  n 

POST- MI  SSI ON  ANALYSIS  YLRSION 

OuMo  BlOCk  PAT  A ,  I'ONTOK ,  PL.AkP. 
ll'Al  k.  ATRAk,  RTRAk ,  CONI' Mv 
RLSOI.V ,  COMB  IN  ;nul  STRAK  an’ 
Mention  l  in  Both  versions) 


♦RSSM 


SCRRT 
SQUEZ  3 


'IS  i:  DATA  0 


TSEDRTR  PROG 
tFORT 

♦  TITL  FILE  TSEDRTR  -  DRTR  INPUT  SUBROUTINE  FOR  ETSE 
SUBROUTINE  PPRDSC<,RNG> 

IMPLICIT  INTEGERS  <R-Z> 

1NTEGER*4  PPRI<lu28>,  PPR 1 2 < 1028 > ,  RSRVE < 16 > 
INTEGERS  OUT <260 >,  0UT2<260> 

INTEGER*4  PPRRNG 
INTEGERS  RNG<6> 

EQU I VRLENCE  <  PPRRNG,  PPR I <  3  >  > ,  <  PP, PPR I <  2  >  > 

EQUI VRLENCE  <0UT<1>,  PPR1 <1> >,  <0UT2<1>,  PPRI2<1> > 
COMMON  /SECTOR/ I NDCTR 
COMMON  / Z S T 0 RE /PPR I 
COMMON  /RUNSUM/PPR I 2 
RVEN=3 

CRLL  C0NMSG<6,  "PPRDSC  > 

2  REWIND  9 
NRZ=0 

ELEV1=RNG<4 >*11.  37778 
ELEV2=RNG<  5  >  *11.  37778 
IF< I NDCTR  EQ  2>INDCTR  =  0 
IF  <  I  NDCTR  EQ.  LOGO  TO  1 
C W  =  RNG <  3  >  —  RNG <  2  > 

IF  <CW  GT.  J 80 >  CW  —  04  ~  360 

IF  <CW  LT  < -180 >  >  CW  =  360  +  CW 

IF  <CW  LT  O >  CW  =  O 

IF  <CW  GT.  0 >  CW  *  15 
BGNfl  =  11.  37778  *  RNG>.2> 

ENDR  =  11  3777S  *  RNG<3> 

BGNR=MOO  <.  BGNR, 4096 > 

ENDR=MODvENDR,  4096 > 

IF<ENDR  GT  40S6 >ENDH=0 
STPFLG  =  0 

IF  <CW  GT  0  HND  ENDR  LT  BGNR >  STPFLG  =  1 

IF  <CW  EQ  0  AND  ENDR  GT  BGNR>  STPFLG  =  1 

I F < BGNR  GT.  10  RND.  BGNR  LT  4086>GO  TO  1 
BGNR=0 


1 

STPFLG=2 

CONTINUE 

10000 

CONTINUE 

♦ASSM 

FREZE 

CROSS 

COPY 

SVt 

RO 

EQU 

0 

Ri 

EQU 

1 

R2 

EQU 

2 

*3 

EQU 

3 

R4 

EQU 

4 

R5 

EQU 

5 

R6 

EQU 

6 

R14 

EQU 

14 

R15 

EQU 

15 

*13 

EQU 

13 

R7 

EQU 

7 

R8 

EQU 

8 

R9 

EQU  9 

R10 

EQU 

10 

Rll 

EQU 

11 

P4  -> 

trnt  i 

A  O 

0,  RSRVE 

R4,  WA1TREAD+SVC1.  SAD 
R4,  15 

R4, Wfl I TREAD+SVC1  ERD 
R13, WREAD 
Rl,  4 

R0, RPR I <R1> 

R0,  3 
R0,  1 
R0,  8 
RO,  4 
R0,  2 
RO,  1 

R4, PPRIBLK+SVC1.  SflD  GET  BEGINNING  ADDRESS 
R4, R0  COMPUTE  END  ADDRESS 

R4,  PPRIBLK+SV'Cl  ERD  STORE  IN  ERD 
R4, WAITPEAD+SVCl  ERD  SAME  NUMBER 


LETS  ONLY  256  WORDS  OF  VIDEO  OUT 


STM 

L 

RIS 

ST 

BRL 

LIS 

L 

NHI 

RIS 

SLLS 

RIS 

SLLS 

SIS 

L 


SI  R4,  2( 
ST  F 

L  F 

RR  F 

ST  F 

LETS  ONLY 


?04S 

R4, OUTBLK+SVC1  ERD 

R4, PPRIBLK2+SVC1.  SRD  GET  NEXT  BEGINNING  RDDR 
R4, R0  COMPUTE  END  ADDRESS 

R4, PPRIBLK2+SVC1  ERD  STORE  IN  SVC  BLOCK 


WORDS 


VIDEO  OUT 


— 

ST 

R4, 0UTBLK2 

LDA I 

R15,  RPR  I 

LDAI 

R14,  RPR  1 2 

LCS 

R6,  1 

DETECT 

BAL 

R13, WREAD 

LH 

Rl,  8  <  RICO 

—  _ 

NHI 

Rl,  X  FFF 

BAL 

R13, WREAD 

LH 

R2, 8<R15> 

NHI 

R2,  X  FF  F' 

CR 

Rl ,  R2 

BL 

CW  I AE 

LIS 

R3,  0 

B 

WHATIZIT 

CWISE 

LIS 

R3,  15 

WHATIZIT 

LH 

P4,  INDCTR 

BZ 

EDETECT 

i~H 

R3,  CW 

BNE 

DE  TECT 

OR 

R3,  Rl 

BZ 

CCW 

LH 

R5,  STPF’LG 

THI 

R5,  2 

BNZ 

CWCASE2 

£MPi 

CH 

R2, BGNA 

BL 

WAIT 

B 

DETECT 

CWCASE2 

CHI 

R2. X  800 

BL 

CMP1 

SHI 

R2,  X  10O0 

_ - _  _ 

B 

CMP1 

CCW 

LH 

P5, STPFLG 

THI 

R5,  2 

BNZ 

CCWCASE2 

01P2 

CH 

R2,  BGNA 

BL 

DETECT 

WO  T  T 

POI 

pi  ?.  i.|pppr> 

- - — - - 

SET  R6  TO  -1  FOR  COUNTER 
READ  IN  RN  AZIMUTH 
LORD  IN  AZIMUTH  DATA 
AND  OUT  UNWANTED  BITS 
GET  ANOTHER  AZIMUTH 
GET  AZIMUTH  DATA 
AND  OUT  UNWANTED  BITS 
COMPARE  TWO  AZIMUTHS 
IF  R1CR2,  RADAR  IS  GOING  CWISE 
DIRECTION  FLAG 
CONTINUE 
DIR  15LAG  =  CW 
SECTOR  SCAN  OR  FULL  CIRCLE 71 
FULL  CIRCLE 

IS  DIRECTION  OF  ROTATION  CORRECT 
WRONG  DIRECTION.  WAIT 
WHICH  DIRECTION  IS  IT? 

COUNTER  CLOCKW I ZE 

CASE  2? 

CW  CASE  2 
ANGLE  <  BGNA? 

YES.  GET  READY 
NO,  TRY  AGAIN 
ANGLE  >  180  -* 

NO, ALL  Ok 

YES,  SUBTRACT  360 


CHECK  FOR  CASE 

ANGLE  >  BGNA? 
NO,  TRY  AGAIN 
vpc.  opt  pconu 


LH 

l  -.X  •  1'  C- 

R2  SCR 15 > 

NHI 

R2- X'FFF" 

CR 

Rl,  R2 

BL 

CW3 

LIS 

R3,  0 

CMP6 

CH 

R3,  CW 

BNE 

DETECT 

OR 

R3,  R3 

BZ 

CCW2 

THI 

R5,  2 

BNZ 

WCWCASE2 

CMP3 

CH 

R2,  BGNA 

BL 

WAIT 

READ1 

LH 

R2, 12<R15> 

NHI 

R2,  X  'FFF' 

CHI 

R2,  681 

. 

BP 

ZERO 

ELDET 

CH 

R2,  ELEV1 

BM 

DETECT 

CH 

R2,  ELEV2 

BP 

DETECT 

AIS 

R6,  1 

BNP 

DETECT 

READ 

BRL 

R13, GOREAD 

LH 

R2,  PPRANG 

NHI 

R2,  X'  FFF' 

OR 

R3,  R3 

BZ 

CCWI 

_ 

THI 

R5,  1 

BNZ 

CWCASE1 

*  WANT  SEVERAL  FPIS 

B 

READ 

CH 

R2,  ENDA 

BNL 

DONE 

- 

B 

READ 

CW3 

LIS 

R3,  15 

B 

CMP6 

ZERO 

LIS 

R2,  O 

B 

ELDET 

CWCASE1 

CH 

R2,  BGNA 

BNL 

READ 

*  WANT 

SEVERAL  PPIS 

B 

READ 

CH 

R2,  ENDA 

BNL 

DONE 

B 

READ 

CCWI 

THI 

R5,  1 

BNZ 

CCWCASE1 

CH 

R2,  ENDA 

BL 

DONE 

B 

READ 

CCWCASE1 

CH 

R2,  BGNA 

—  - 

BL  - 

-READ 

CH 

R2, ENDA 

BL 

DONE 

B 

READ 

CCW2 

THI 

R5,  2 

BNZ 

WWCASE2 

-CMP4 

CH 

R2,  BGNA 

BNL 

WAIT 

B 

READ1 

MWCASE2 

CHI 

R2, X' 800 ' 

BNL 

CMP4 

AH  I 

R2, X'1000' 

D 

CMoa 

GET  NEXT  AZIMUTH 
FIND  OUT  UNWANTED  BITS 


C-fiSE  2? 

VES,  BRANCH 
ANGLE  >  BGNA? 

NO,  KEEP  WAITING 


INCREMENT  COUNTER 
DO  IT  TWICE  TO  BE  SURE! 
VES, START  READING 


BRANCH  IF  COUNTER  CLOCKWISE 
CASE  1? 

VES,  BRANCH 


ANGLE  >  ENDA? 
VES,  ALL  DONE 
NO,  KEEP  READING 


59.  5  DEG  =  0 

ANGLE  <  BGNA? 

VES,  KEEP  READING 


NO,  ARE  WE  DONE  VET? 

VES 

NO,  KEEP  READING 
CHECK  FOR  CASE  1 

ANGLE  <  ENDA? 

VES,  ALL  DONE 
NO,  CONTINUE 
ANGLE  >  BGNA? 

NO,  KEEP  READING 

VES,  CHECK  FOR  FINISHED 


OK 

NOT  VET 


D-4 


CCWCASE2 

CHI 

R2,  X'  800 

ANGLE  <  180? 

BNL 

CMP2 

NO,  ALL  OK 

AH  I 

R2,  X  / 1000  *** 

VES,  ADD  360 

B 

CMP2 

NOW  CHECK  FOR  BEGINNING 

WCWCASE2 

CHI 

R2,  X  '800" 

ANGLE  >  180? 

ESL 

CMP3 

NO, ALL  OK 

SHI 

R2, X- 1000" 

VES,  SUBTRACK  360 

B 

CNP3 

EDETECT 

BAL 

R13, WREAD 

READ  IN  NEW  AZIMUTH 

LH 

R2,  12<R15> 

READ  IN  ELEVATION 

NHI 

R2,  X  '  FFF ' 

AND  OUT  UNWANTED  BITS 

— 

CHI 

R2,  681 

BP 

ZEROl 

CPEV 

CH 

R2, ELEV1 

BL 

EDETECT 

CH 

R2, ELEV2 

BP 

EDETECT 

WITHIN  RANGE? 

LH 

R2, PPRANG 

VES,  GET  AZIMUTH 

NHI 

R2,  X  FFF  -' 

CHI 

R2,  6 

BL 

FUDGE 

READ2 

BAL 

R13, GOREAD 

START  READING 

LH 

R4, PPRANG 

- -  . 

NHI 

R4,  X  FFF 

NXT 

OR 

R3-  R3 

BZ 

CCL0K4 

COUNTERCLOCKW I SE 

CR 

R4,  R2 

BL 

NXT2 

B 

READ2 

•ZERO! 

LIS 

R2,  0 

B 

CPEV 

FUDGE 

LIS 

R2,  6 

IF  0,  MAKE  IT  6 

B 

READ2 

NXT2 

BAL 

R13, GOREAD 

LH 

R4, PPRANG 

GET  NEW  AZIMUTH 

_ _ 

NHI 

R4,  X'  FFF  " 

AND  OUT  UNWANTED  BITS 

OR  R3, 

R3 

BZ 

CCW4 

WANT  SEVERAL  PPIS 


CR 

R4,  R2 

ANG  >  ENDA? 

.. 

BP 

DONE 

VES,  FINISHED 

B 

NXT  2 

NO, KEEP  READING 

CCL0K4 

CR 

R4,  R2 

BP 

NXT2 

B 

READ2 

CCW4 

CR 

R4,  R2 

-  . 

BL 

DONE 

B 

NXT2 

DONE 

LH 

R15, PPRI 

NHI 

R15, X ' 0FFF " 

STH 

R15,  PPRI 

SVC 

1,  OUTBLK 

.  . 

LH 

R0,  OUTBLK+SVC1.  STA 

BNZ 

ERROR 

LM 

0, RSAVE 

♦FORT 

RETURN 

♦ASSM 

-  .... 

.  ALIGN 

4 

WREAD 

SVC 

1, WA I TREAD 

LH 

R0,  WAITREAD+SVC1.  STA  READ  RETURNED  STATUS 

BNZ 

ERROR 

IF  NOT  ZERO,  ERROR 

__ 

BR 

R13 

GOREAD 

SVC 

1, PPRIBLK2 

READ  IN  ONE  AZIMUTH 

Ti  »tr  OI  irno^r^-  mrr*  r>r» 


n« 


L13  ii-  CUUNIER  FOR  MAIN  BUFFER 

LHI  R4,  16  COUNTER  FOR  NEW  BUFFER 


AVE 

LIS 

R8,  O  ZERO 

POWER  ACCUMULATOR 

LIS 

RIO, 0  ZERO 

VELOCITV  ACCUNLTR 

LIS 

Rl, 3  AVERAGING  COUNTER 

AVER 

LHL 

R7, PPRICR4) 

LOAD  POWER  HALFWORD 

NHI 

R7,  X'lFF' 

GET  ONLY  THE  POWER 

AR 

RS,  R7 

ADD  TO  ACCUMLTR 

AIS 

R4,  2 

INCREMENT  ARRAY  POINTER 

LHL 

R7,  PPRI  <R4> 

LOAD  VELZVAR  HALFWORD 

NHI 

R7,  X'FFFF' 

AR 

R10,  R7 

INCREMENT  FOR  NEXT  HALFWORD 

AIS 

R4,  2 

INCREMENT  FOR  NEXT  HALFWORD 

SIS 

Rl,  1 

SUBTRACT  FOR  INNER  LOOP 

BNZ 

AVER 

BRANCH  IF  NOT  3  ADDED  UP 

DH 

R8,  AVEN 

DIVIDE  POWER 

NHI 

R9,  X  IFF 

SLA 

R9,  16 

SHIFT  POWER  TO  LEFT 

SRA 

RIO,  8 

DH 

RIO, AVEN 

DIVIDE  VELOCITY 

SLA 

Rll,  8 

NHI 

Rll, X'FFFF' 

GET  ONLY  LOWER  HALFWORD 

AR 

R9,  Rll 

PUT  TOGETHER 

- 

AIS 

R12,  4 

INCREMENT  NEW  ARRAY  POINTER 

ST  R9,  PPRI<R12>  PUT  INTO  ARRAY 
CLHI  R4, 3072 
BP  AVEEND 

B  AVE 


* 


AVEEND 

SVC 

1, OUTBLK 

OUTPUT  LA  T  AZIMUTH 

SVC 

3 . WAITBLK 

WAIT  FOR  FIRST  READ  TO  FINISH 

SVC 

1, PPRIBLK 

READ  NEXT  AZIMUTH 

SVC 

1 • 0UTBLK2 

OUTPUT  LAST  AZIMUTH 

SVC 

1. WAITBLK 

WAIT  FOR  READ  TO  FINISH 

LH 

PO, PPRIBLK 

2+SVC1.  STA  READ  STATUS 

BNZ 

ERROR 

IF  NOT  ZERO,  ERROR 

LH 

RU,  PPRIBLK  +SVC1.  STA 

BNZ 

ERROR 

LH 

PO,  0UTBLK+SVC1.  STA 

BNZ 

ERROR 

LH 

PO, 0UTBLK2 

+SVC1.  STA 

BNZ 

ERROR 

LH 

R6,  NAZ 

AIS 

R6,  2 

ADD  TWO  TO  AZIMUTH  CTR 

STH 

R6,  NAZ 

CHI 

P6,  1760 

TOO  MANY  AZIMUTHS? 

BLR 

113 

NO,  KEEP  GOING 

LIS 

F6,  2 

YES,  INDCTR=2 

STH 

F  6,  INDCTR 

B 

CONE 

QUIT 

ERROR 

SVC 

1 , ERRCODE 

CONVERT  ERROR  CODE 

SVC 

2,  ERRBLOK 

OUTPUT  MSG  TO  CONSOLE 

SVC 

2,  PAUSE 

TASK  PAUSED 

.  . 

LM 

0,  RSAVE 

RESTORE  FORTRAN  REGISTERS 

****P2 

IS  FORTRAN  STMT  NO 

.  2 

B 

*P2 

START  OVER 

ALIGN 

4 

PAUSE 

EQU 

* 

DB 

0,  1 

PAUSE 

_ 

ALIGN 

4 

ERRCODE 

EQU 

* 

OPT 

DB 

0 

DB 

6 

CODE  6 

. 

DC 

H'0 

DC 

A<ERCD) 

DESTINATION 

A 

r>-(i 

c 

t 

< 


ERRBLOk 

EQU 

♦ 

DB 

0,  7 

DC 

H '  15  ' 

ERCD 

DC 

C' 

DC 

C  I/O  ERROR 

RLIGN 

4 

WRITBLk 

EQU 

* 

DB 

X' Os¬ 

DB 

lo 

DB 

0,0 

DSF 

5 

RLIGN 

4 

» 

I 

I 

I 

I 

I 

> 

C 

> 

i 


WRITRERO  EQU  * 


DB 

X '  59  •" 

DB 

10 

DB 

0,  0 

DC 

RCPPRI > 

-  DC 

RCPPRI > 

DSF 

3 

RLIGN  4 

PPRIBLK2  EQU 

♦ 

DB 

X  51- 

DB 

10 

DB 

0,  0 

DC 

RCPPRI2: 

DC 

RCPPRI2' 

DSF 

3 

RLIGN 

1  4 

0UTBLR2  EQU 

♦ 

-  -  OB 

X  -  31  - 

DB 

9 

DB 

0,  0 

DC 

R  <  OUT  2  > 

DC 

R<0UT2> 

DSF 

—  RLIGN 

4 

PPRIBLK  EQU 

♦ 

DB 

X '  51  - 

DB 

10 

DB 

0,  0 

DC 

rcppri  :» 

-  DC 

RCPPRI,'- 

DSF 

3 

RLIGN 

4 

OUTBLK  EQU 

* 

DB 

X  -  31  - 

DB 

9 

OUTST  DB 

0,0 

DC 

RCOUT) 

DC 

fl<OUT> 

DSF 

3 

PRINT  CONSOLE 
CODE  ? 

PRINT  15  CHRS 
ERROR  CODE 


WRIT  ONLV 
LU  10 


RERD  RND  WRIT 
LU  10 


RERD 
LU  10 


LU  9 


RERD 
LU  10 


WRITE 
LU  9 


♦FORT 


RETURN 

END 


♦ 

07/19/79  12  49 : 01 

♦♦♦LISTING  FOR  EfiTl :  CRANE.  FTN 
*N 

*ASSM 

ERTRMAP  PROG 
*FORT 

♦TITL  FILE  TSERMAF  -  PRINT  OUT  DATA  FIELDS-CHANGED  FOR  ERT  READ  BV  CLB 

IMPLICIT  INTEGER+2  <I-N> 

INTEGER*2  NRC,  VAR,  STORE < 10 >,  PLACE,  OLDATA 
INTEGERS  ZEE < 1024 >, ANC< 1626 >,  IDT  I  ME,  IZF,  ITZ,  IZB,  IZS 
INTEGERS  TWENTV,  ELEVEN,  DAV,  HOUR,  MINUTE,  SECOND,  TP,  ELEVAT,  AZ 
INTEGERS  AZIM,  TC,  TA,  T 

INTEGERS  V,  I,  THETA1,  NRC1,  II,  THETA,  RHO,  STOP,  GRND 
INTEGERS  MEAN,  POWER,  SIGMA,  TP2,  TPS,  SEGNO,  Q,  J,  ZTH,  MMU 
INTEGERS  BEGIN,  SUM1,  SUM2,  JMIN,  M,  K,  L,  I 
INTEGERS  VI,  ZERO,  TWO,  BETA 
INTEGERS  RSAVE<16>,  R1SAV 

INTEGERS  IREF<1024>,  I VEL<1024>,  I  VAR  <1024),  I VELLC256) 

+  , JDVELC256) 

INTEGER*2  RE, HR, TL1,  TL2,  RQUANT 
REAL  RNN, RRAREA,  AXX 

REAL  FCTMIN, AZT,  ELEVAA,  BGNA,  AZCHK,  ENDA 
INTEGER+2  FLAG, IQ, IB, IE, I M 
INT£GER*2I0UT,  DBB,  MAXV,  MAXS,  SLOPE 
INTEGERS  NRCEAD 

COMMON  /TLIS/  T,  TA,  JDAV,  JHR,  JMIN,  JSEC,  I  DAV,  I  HR,  I  MIN,  I  SEC 

COMMON  /CONST/  VMISW,  DIV,  VMAG,  VMISWM,  ZDIV,  AD IV,  HI,  A2,  A3,  Bl,  B2,  HDIV 

COMMON  /SWITCH/  IC1<44>,  IC2<44>,  TC<19S0),  IPCK5400), 

+  IPC2C5400),  IPC3C5400),  IPCNT<19S0>,  IPTC<44>, 

♦  NEMB, NEMC,  NAC 

COMMON/ AZ2/S I NA,  COSA,  DELTAS,  ISCANF- NEL, RI,  SA 
CGMMON/AZM/AZT,  RZLAST ,  AZSTAR,  NA,  ELEVAA 
COMMON/REFL/RE  < 1025 ) ,  HR <258),  NCL,  NID,  NIDP,  INCL 
X,  I  MX,  IMN,  TL.,  TL2,  RQUANT,  IDVEL<258> 

COMMON  /PNTfS/  NCMX,  NVMIN,  NUMX,  IELSN, NSCAN,  IESNL,  NVSCN,  NT 
COMMON  /INTI  /  MHSN,  MNSN,  HM,  FNSN 

COMMON/EXTRA/RHO,  GRND,  ZTH,  BETA,  K,  RRAREA <28S>, RH02<16>,  AXX<28S>, 
1IPRFM 

COMMON/ZSTORE/ANC 

COMMON  /ZLOOK/  IZOFF, 2ARV<91> 

COMMON  /ECONST/  EARTH,  VMK 

COMMON  /MAPPAR/  DAV,  HOUR, MINUTE, SECOND,  DBB, MAXV, MAXS,  IOUT,  SLOPE 

COMMON  /CNT /  COSPHI ,  SINPHI ,  C0SPH2 

COMMON  /DATA1/  ECL<224),  NOFST, KOFST,  I CLAD,  NANI 

COMMON  /DA rA2/  VCL<736>, MXVC, NVC1 

COMMON  /DATA3/  VR <  192),  MXVR,  NVR1 

COMMON  /NVI.IS/  NVARM,  NCARM,  NVO,  ICO,  10,  JO,  JVR,  KTL 

COMMON  /FILTER/  TATRMN,  AREAMN,  DAZM 

EQU I VALENCE < RE < 2 > ,  IREF<1>>,  < IDVELC2), JDVEL<1>  > 

EQU I  VALENCE  <  ANC  <  5  > ,  ZEE  <  1  >  > 

DATA  TWENTV/18/ 

DATA  TWO/2/ 

DATA  PCTM1N/  05/, MMU/O/,  SIGMA/0/,  OLDATA/0/ 

DATA  ZERO/0/ 

C 

C  *******  INITIALIZE  ARRAV 

c - 

T-0 

RZLAST— 999 
DO  901  J-l,  MXVC 

“901  VCL<J>-0  P-8 

DO  902  J-l,  MXVR 


1 


VM I SUM* < VM I SM-1.  >/DIV 

RE<1>=0 

JDVEL<l>-0 

JDVEL  <  256  >  =0 

ICC-0 

C  CALL  C0NMSGC7,  'TSERMRPO 
TL1-TL1/RQUANT 
TL2-TL2/RQUANT 
REMIND  4 
REMIND  8 
DO  2  1-1,1024 
IREF<n«0~ 

I VEL< I )*0 
2  I VRR< I >— 0 
DO  3  IX-1,  91 
I I— IX-I20FF 
2X=FLORTCII>/10. 

3  2RRVC  1X5=10.  +*ZX 
00001  CONTINUE 
DB-65.  28 

750  FORMRT < IX,  'ENTER  PRF,  0=768,  1=922,  2=1075, 3=1229 '  > 
RCKM=  075 
FLRG=0 

-  READ  <  7,  333  >  I  OUT,  DBB,  MRXV,  MRXS,  SLOPE 

333  FORMRT  < I3> 

READ <7, 333>0LDATA 

01002  FORMRT <.//'  R2IM  ELEV  RAN  REF  VEL  STD'>> 

BITVEL-MRXV 
B I TVRR=MRXS 

-  REMIND  9 

K=0 


10000 

CONTINUE 

1ASSM 

R0 

EQU 

0 

Rl 

EQU 

1 

R2 

EQU 

2 

R3 

EQU 

3 

R4 

EQU 

4 

R5 

EQU 

5 

R6 

EQU 

6 

R7 

EQU 

7 

R8 

EQU 

8 

R9 

EQU 

9 

RIO 

EQU 

10 

Rll 

EQU 

11 

R12 

EQU 

12 

R13 

EQU 

13 

R14 

EQU 

14 

R15 

EQU 

15 

STM 

R0,  RSAVE 

FREZE 

COPV 

SVC1. 

L  R3, ANCBLK+SVC1.  SAD  ANCBLK  START 

AIS  '  R3, 15  READ  RNCILLRRV  ONLV 

ST  R3,  ANCBLK+SVC1.  EAD 

SVC  1,  ANCBLK  THE  READ 

LH  R0, ANCBLK+SVC1.  STB 
BNZ  ERROR 

LIS  Rl,  0 

STH  Rl, STOP 

LIS  R5, 0  MORKING  REGISTER 

LIS  R9, 10  MULTIPLICAND 

LH1  Rll,  100  MULTIPLICAND 

L  RO,  ANC<R1>  DATA 

STBR  R0,  R5  HOUR 


KJLO/  AO 

NR 

R10,  R5 

1  HOUR 

SRLS 

R5,  4 

10  HOUR 

NHI 

R5,  3 

MHR 

R5,  R9 

<*10> 

AHR 

R5, R10 

TOTAL 

STH 

R5,  HOUR 

HOURS  .  " . .  " 

EXBR 

R0,  R0 

STBR 

R0,  R5 

MINUTES 

LIS 

R10,  15 

MASK 

NR 

R19,  R5 

1  MINUTE 

SRLS 

R5,  4 

10  MINUTE 

NHI 

R5/7 

MASK  '  ‘ 

MHR 

R5,  R9 

<*10> 

flHR 

R5,  R10 

TOTAL  MINUTES 

STH 

R5, MINUTE 

EXHR 

R0,  R0 

STBR 

P0,  R5 

DAVS 

LIS 

R10,  15  " 

MASK 

NR 

R10,  R5 

1  DAV 

SRLS 

R5,  4 

10  DAV 

MHR 

R5,  R9 

<*10> 

flHR 

R5,  R10 

SRLS 

R0,  8 

100  DAVS 

NHI 

R0,  X'F  ' 

MASK 

MHR 

R0,  Rll 

<*100> 

flHR 

R5,  R0 

TOTAL  DAVS 

STH 

R5,  DRV 

DAVS 

flIS 

Rl,  4 

NEW  DATA 

L 

R0, ANCCR1) 

LIS 

R5,  0 

JUST  TO  BE  SURE 

STBR 

R0,  R5 

NRC  AND  TP 

LIS 

R10,  2 

MflSK 

NR 

R10,  R5 

NRC 

flIS 

R10, 1 

LHI 

R12,  256 

MHR 

R10,  R12 

NRC=  <  NRC+1 >  *256 

LHI  R10, 256 

STH 

R10,  NRC1 

STORE 

SLHLS 

10,  2 

MULTIPLV  BV  4  FOR  BVTE  COUNT 

flIS 

10, 15 

fl 

10, NRCBLK+SVC1.  SAD 

ST 

10,  NRCBLK+SVC1.  EAD 

STORE  END  RDDRESS  FOR  READ 

ST 

10, NRCEflD 

STORE  AWRV  FOR  SECOND  READ 

NHI 

R5,  X'C0" 

MASK  TP 

SRLS 

R5,  6 

STH 

R5,  TP 

STORE (UNFIXED) 

EXHR 

RO,  R0 

STBR 

R0,  R5 

SECONDS 

LIS 

R10, 15 

MASK 

NR 

R10,  R5 

1  SECOND 

SRLS 

R5,  4 

NHI 

F5,  7 

MASK  10  SECONDS 

MHR 

*9,  R5 

<*10> 

flHR 

R9,  R10 

TOTAL  SECONDS 

STH 

R9,  SECOND 

STORE 

RIS 

Rl,  e 

MORE  DRTR 

LHL 

R0,  RNCCRl) 

ELEVATION 

NHI 

R0,  X'FFF' 

ANGLE 

STH 

R0, ELEVRT 

STORE  ANGLE 

LM 

R0,  RSAVE  - 

~  .  -  • 

♦FORT 

NRC-NRC1 

GRND-PCTMIN*NRC 

- REWIND 

17  ELEVEN-0 


D-10 


I 

f  DR=<2**TP>*RCKM 

K=K+1 

NRC1=NRC 
f  10001  CONTINUE 

1.RSSM 

1 

STM 

R0,  RSflVE 

fa 

ff 

SVC 

1-  NRCBLK 

READ  IN  DATA  1 

LH 

R0, NRCBLK+SVC1 

STA  LOAD  IN  STATUS  1 

BNZ 

ERROR 

IF  NOT  ZERO,  I/O  ERROR  1 

f 

LIS 

Rl,  0 

1 1 

LH 

R0, ANC<Rl> 

1 

BM 

MINUS 

IF  STO P=0 

LIS 

Rl,  1 

STORE  1  , 

• 

STH 

Rl, STOP 

1  5 

MINUS 

LIS 

Rl,  8 

GET  AZIMUTH 

LH 

R0,  rnc<ri> 

ANGLE  ; 

i 

NHI 

R0,  X'FFF' 

MASK 

b 

STH 

R0,  flZ 

STORE 

* 

AIS  Rl, 

4 

LHL  R0, 

flNCCRl? 

ELEVATION 

: 

NHI  R0, 

X'FFF' 

ANGLE 

a 

STH  R0, 

ELES'AT 

STORE  IT 

t  * 

LHL 

RQ,  NRC1 

MAX.  NO.  OF  CELLS 

1 

LHI 

R2,  32767 

MAN.  HALFWORD 

I  ft 

LDRI 

R15, GRND 

1 

LHL 

Rl,  0  C R15  > 

CELL  COUNTER  | 

I 

LIS 

R2>  0 

SUM  ZR  1 

1  f 

LIS 

R4,  7 

FOR  1 

1 

LIS 

R5,  0 

SHIFTING  1 

t 

LIS 

R6,  0 

SUM  ZR2/128  I 

K  £ 

LDfll 

R15,  ZTH 

1 

K  £* 

LHL 

R?,  0<R15) 

MIN.  POWER  1 

1  C 

LIS 

R9,0 

CELLS  PER  SEG.  9 

LDfll 

R15, RHO 

H 

■ 

LHL 

R10, 0<R15> 

MIN.  SEGMENT  SIZE  jj 

■ 

STH 

R10,  PLACE 

INITIALIZE  PLACE  | 

■  # 

SIS 

R10,  1 

E 

i  ^ 

LHL 

Rll, MMU 

MINIMUM  MEAN  | 

■ 

LHL 

R12,  SIGMA 

MINIMUM  S'AR I ANCE  P 

I  • 

B 

POW 

j 

MAIN  LOOP 

1  •  LOOP 

flR 

Rl,  R10 

i 

K  NLOOP 

flIS 

Rl,  1 

NEXT  CELL  I 

■  ft 

CR 

Rl,  R8 

IF  DONE, 

■ 

BNL 

RERLLV 

LEAVE 

■  POU 

LR 

R15,  Rl 

■  ft 

SLLS 

R15,  2 

B 

L 

R5, ZEE<R15> 

RAW  DATA 

B 

EXHR 

R0,  R5 

POWER 

1  0 

LH 

BNZ 

R2,  OLDflTfl 

NOS 

B 

NHI 

R0, X'lFF' 

1 

■  ftft- 

STH 

R0, IREF<R1> 

STORE  REFLECTIVITY 

B  CONT 

EXHR 

R15,  R3 

PUT  MEAN  IN  TOP  OF  R15 

B 

SRfl 

R15,  24 

SHIFT  DOWN  24  WITH  SIGN  EXT. 

1  < 

BNMS 

POSITIVE 

ABSOLUTE  VALUE 

XHI 

R15, X'FFFF' 

COMPLEMENT 

B  —  •■ 

AIS 

R15, 1 

TWO'S  COMPLEMENT 

B  ft 

STH 

R15, IVEL<R1> 

STORE  VELOCITY 

■  POSITIVE  NHI 

R3, X'FF' 

VARIANCE 

m 

STH 

R5, IVAR<R1> 

STORE  VARIANCE 

1  • 

B 

NLOOP 

„ 

■  N09 

NHI 

__J 

R0,  X'FF' 

n- 1 1 

WMi 


€ 


* 


»nc‘>  i^r  s.  ivi.  / 

B  CONT 


c 


*  DATA  STORING  ROUTINE 

* 

REALLV  LM  R0, RSAVE 
♦FORT 

DO  101  1=1,236 

IF  CIREFCO.  LE.  0>  GO  TO  111 

IR£F<I)  =  II>EF<I>*100.  /256.  -DB 

IF<  IREF  <  I  LT.  -39  >  I  REF  <  I )  =  I  REF  <  I  >  +100 

I VEL  < I >  =  I VEL < I >  *6 1 TVEL/128. 

IDVEL< I >=IVEL< I >-IVELL< I > 

IVELL< I >=IVEL< I > 

C  F4W0RD- 1 VAR  < I > ♦B I TVAR/256. 

C  IVAR< I >  =  SORT  <F4W0RD> 

GO  TO  101 
00111  IREF  < I >=0 
IVEL< I >=0 
I VAR  < I ) =0 
101  CONTINUE 
00500  CONTINUE 
ICC=ICC+1 

ELEVAA=ELEVAT/11  3???S 
I ELSN= I F I X < ELEV A A > 

T=<  < DAV*24+H0UR > *60+M I NUTE ) *60+SECOND 
A2T=AZ/11.  37778 
FLAG=FLAG+1 
A=A2T*.  01743 
A2CHK=AZT 

IF<K.  NE.  1>  GO  TO  105 

BGNA=A2T 

ENDA=A2T +359. 

105  CONTINUE 

IF<K  LT.  180)  GO  TO  106 
IF <AZT.  GT.  180.  >  GO  TO  106 
A2CHK=AZT +359. 

106  CONTINUE 

IFCA2CHK.  GT.  ENDA>  K=1 
SINA=SIN< A) 

COSA=COS<A> 

DELTA2=0.  0191987 


NA=K 


» 


RE  < 258 ) =0 


IF<NA.  EQ.  1)NAC=1 

NAC=NAC+1 

IFCNAC.  GT.  1>NAC=0 

IFCNA.  NE.  1.  OR.  A2LAST.  LT.  -990.  >  GO  TO  224 

IDAV=DAV 

I HR=HOUR 

IMIN=MINUTE 

ISEC-SECONC 

PH I *ELEVAA*DAZM 

COSPHI=COS<PHI> 

SINPHI=SIN<PHI> 

C0SPH2-C0SPH I *COSPH I *E ARTH 
CALL  TRACK 

223  XFCNSCAN.  EQ.  4>  GO  TO  800 

224  CALL  CONTOf’ 

AZLAST-AZT 

-  IF (STOP  NE  1>  GO  TO  17  -  - 

800  CALL  STRAK 
STOP 


♦ftSSM 

-  ALIGN  4 

ERROR  SVC  2.  ERRCODE 


D-12 

DECODE  ERROR  BITS 


SVC  2,  fcKKfciCUK  UL'  I  t"'U  I  tK'KUK  PliU  lu  LUNiUUt 

SVC  2, PAUSE 

LM  0, RSRVE  RETURN  TO  FORTRRN 

♦FORT 

IF<K.  EQ,  0>  00  TO  10000 
IFCK.  EQ.  1>  GO  TO  10001 
00  TO  1 

10003  CONTINUE 
♦RSSM 


RLIGN 

4 

ERRCODE 

EQU 

♦ 

DB 

0 

~ *"  •' 

DB  ' 

6 

DC 

H'0' 

DC 

fl<ERCD> 

RLIGN 

4 

ERRBLOK 

EQU 

* 

DB 

0 

DB 

7 

DC 

H'28' 

ERCD 

DC 

C' 

DC 

C'l/O  ERROR 

IN  RTN  RANGE' 

RLIGN 

4 

PRUSE 

EQU 

* 

' 

DB 

0,  1 

RLIGN 

4 

RNCBLK 

DB 

X  ■'  59 ' 

RERD  BLOCK  FOR  LITTLE  RERD 

DB 

9 

LU 

DB 

0.  0 

DC 

R  <  RNC  > 

START  ADDRESS 

DC 

fl<RNC> 

END  ADDRESS 

DSF 

3 

NRCBLK 

DB 

X'59' 

READ  BLOCK  FOR  BIG  RERD 

DB 

9 

LU 

DB 

0,  0 

STRTUS 

DC 

RCRNC? 

START  ADDRESS 

DC 

R  C RNC  > 

END  ADDRESS 

DSF 

3 

♦FORT 


END 


0-0000 


€ 


(  07/13/79  12 : 51 : 51 

♦♦♦LISTING  FOR  ERT1 : INPRRM  FTN 

*N 

BLOCK  DATA 

FOR  PROGRFiM  EXTRAD  ERT  NO.  162 
VERSION  4.  0  LEVEL  780301 
JHW  IBM370 

IMPLICIT  INTEGER+2  <I-N) 

INTEGERS  U,  TL1,  TL2,  HR,  RQUANT,  T,  TM 

COMMON  /CDRflVS/  IC<32,  10),  C<32,  9),  IDC32,  10),  D<32,  9),  IM,  JM 
COMMON  /CONST/  VMISW,  DIV,  VMflG,  VMISWM,  ZD  IV,  RDIV,  Rl,  R2,  R3,  Bl,  B2,  HD  IV 
COMMON  /VPRRM/  VX,  VV 

COMMON  /PNTRS/  NCMX,  NVMIN,  NUMX,  IELSN,  NSCRN,  IESNL,  NVSCN,  NT 
COMMON  /ECONST/  ERRTH,  VMK 

COMMON  /NVLIS/  NVRRM,  NCRRM,  NVO,  ICO,  10,  JO,  JVR,  KTL 
COMMON  /2L00K/  IZOFF,  ZflRYOi) 

COMMON  /INTL/  MHSN, MNSN,  HM,  FNSN 
COMMON  /DVRL/  DELR 
COMMON  /CNTRS/  IRTR<5),  MRTR 
COMMON  /FLGS/  ITVPE,  PR1,  PR2,  PRIN2 
COMMON  /TMRX/  TM 

COMMON  /R2M/  RZMUTH,  R2LRST,  RZSTflR,  NR,  ELEVRT 
COMMON  /DRTR1/  ECL<224), NOFST, KOFST,  ICLRD,  NRN1 
COMMON  /DRTR2/  VCLC736),  MXVC,  NVC1 
COMMON  /DRTR3/  VR<192),  MXVR,  NVR1 
COMMON/F I LTER/TRTRMN,  RRERMN,  DR2M 

COMMON  /PWORK/  KMRX,  T<80), JMXDB, JMRX,  IRMRX,  IR, JR,  IMXJMX 
COMMON  /RZ 2/SINfl,  COSR,  DELTRZ,  ISCfiNF,  NEL,  RI,  SR 
COMMON  /FIXED/  NPR,  IEMRX,  NFC,  I EM,  JEM 
COMMON  /PRSTOR  /  NUP, TRTRC1400),  NUMRX,  IRCTC70), 

+  I DC  <  22 ) ,  I PRNG <  34  > , KRMRX,  MXTR 

COMMON/REi  L/  W<  1025),  HR <258),  NCL,  NID,  NIDP,  INCL, 

X  I MX.  INN, TL1, TL2, RQURNT ,  IDVEL<25S) 

COMMON/THRESH/LDV 
LOGICRL  PR1, PR2, PRIN2 

DRTR  ECL/224*0.  0/, VCL/736+0.  0/, VR/192^0  0/ 

DflTR  TL1/3©/,  TL2/60/ 

DRTR  NOFST/16/,  KOFST/7/,  I CLAD/112/ 

DRTR  DRZM/0  0174533/,  ISCRNF/0/ 

DflTR  RI/0.  /,  SR/900.  /,  RQURNT/1/ 

DRTR  LDV/3/,  MXVR/192/,  NVR1/6/,  MXVC/736/,  NVC1/23/ 

DRTR  NPfl/2/,  IEMAX/22/,  NFC/2/ 

DRTR  KMRX/45/,  JMXDB/80/, JMRX/34/,  IRMRX/2700/,  IR/15/, JR/45/, 

♦  IMXJMX/60/,  NCL/258/,  KRMRX/990/ 

DflTR  N I D/250/, NUP/4/, N I DP/70/, NUMRX/20/,  MXTR/1400/ 

DRTR  AREAMN/0.  4/,  I MX/257/,  IMN/2/ 

DRTR  I2OFF/10/,  DELR/0.  00S7/, FNSN/0  009/ 

DRTR  IRTR/9,  8.  5,  0,  0/ 

DRTR  ITVPE/1/,  PR1/  TRUE  /.  PR2/.  FALSE  /,  PRIN2/.  FALSE.  / 

DRTR  NVRRM/32/,  MNSN/5/,  MHSN/7/,  HM/6  5/ 

DRTR  TM/0/,  NCARM/16/,  VMK/1  E-3/ 

DflTR  I M/32/,  JM/9/, NVO/0/,  ICO/0/,  10/0/,  JO/0/ 

DRTR  VMISW/5  /,  DIV/  2/, VMAG/  01/, ZD IV/.  1/, HD IV/  5/ 

DflTR  AD IV/  04/.  fll/.  4/,  A2/.  3/,  03/  3/,  Bl/.  7/,  B2/.  3/ 

DATA  VX/0  /, VV/0  /,  IESNL/0/,  IELSN/0/ 

DATA  ERRTH/6  4857  E-5/,  NT/0/,  NSCAN/0/,  NVSCN/1/ 


END 


P-14 


F/G  17/9 
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AD-A081  061  ENVIRONMENTAL  RESEARCH  AND  TECHNOLOGY  INC  CONCORD  MA 

AUTOMATIC  WEATHER  RADAR  ECHO  ASSESSMENt  AND  TRACKING. (U) 

MAR  79  R  K  CRANE  F19628-78-C-0076 

AFGL-TR-79-0298  NL 


UNCLASSIFIED 


07/19/79  12:52: 43 

♦♦♦LISTING  FOR  ERT1 : CONTOR.  FTN 

*N -  -  -  . .  - 

SUBROUTINE  CONTOR 

c  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

C  JHW  AFGL  SUDBURV  RFlDflR  SUBROUTINE 

C  VERSION  4.  1  LEVEL  78111? 

C  FIND  EVENTS  ALONG  SINGLE  RADIAL 

m  w ' '  w"  "T*  ■?*  ~  w  t'  ^  4^*4*  *  4^  4*  *  ^  4*  *  4*  *  4^  4*4*  4*  ™  4*  *  4*  4*4*4*  4*  4*  4*  4*  4*  4*  4*4*4*4*  4*  “4*4*  4*  4”  4*  4*4*4*4*4**  4*4*4*4* 

IMPLICIT  INTEGER+2  <I-N) 

INTEGER  *2  W,  TL1, TL2,  HR,  RQUANT ,  TC 
+  ,  T,  IC21C22),  IC22<22) 

COMMON  /SWITCH/  IC1<44), IC2C44), TC<1930), IPCK5400), 

+  IPC2C5400),  IPC2C5400),  IPCNTC1930),  IPTCC44), 

- +  NEMB,  NEMC,  NAC 

COMMON/AZM/  RZMUTH, A2LAST ,  AZSTAR,  NA,  ELEVAT 
COMMON  /AZ2/SINA,  COSA,  DELTAZ,  ISCANF, NEL,  RI,  SA 
COMMON/F I LTER/TATRMN,  AREAMN,  DAZM 
COMMON  /FIXED/  NPA,  IEMAX, NFC,  IEM,  JEM 
COMMON  /PRSTOR  /  NUP, TATRC1400),  NUMRX,  IACT<?0), 

- +  IDC<22),  IPRNG<34),  KAMAX,  MXTR 

COMMON/REFL/  W<1025>,  HRC258),  NCL,  NID,  NIDP,  INCL 
X,  I  MX,  IMN,  TL1,  TL2,  RQUANT,  IDVELC258) 

COMMON  /PWORK/  KMAX, T<80>,  JMXDB,  JMAX,  IRMAX,  IR,  JR,  IMXJMX 
COMMON  /THRESH/LDV 
COMMON  /CONPK/  NCEL,  NMR 

c - - 

DATA  RPD/.  017453/ 

C  - 

C  CALL  CONMSG< 6,  'CONTOR O 

C  IFLAG=1 

C  IFCIFLAG.  EQ.  1)WRITE<3,  7)<W< IX),  IX=1,  102) 

C -  FORMAT < IX,  2015) 

IF  <NR.  NE.  1)  GO  TO  61 
C 

C  INITIALIZE. 

C 

TLS=TL1*RQUANT 

TATRMN=AREAMNd.  E06/SA 

NEMC= IEMAX 

NCEL=1 

JEM=0 

61  CONTINUE 

DO  101  K=l,  IEMAX 
-  IDCOO=0 

101  CONTINUE 

DO  102  J«l,  JMAX 
IPRNGC J)=0 

102  CONTINUE 
NEMB=NEMC 

-  NEMC=NAC* IEMAX  -  -  - 

NEM1=NEMC +1 
IEM»0 
IEM2-0 
IP=0 
IPB=0 

C  FIND  EVENTS 

c — 

DO  281  1=2,  NCL-?rf 

- ~IFTWC I).  GT. "TLS)  130  TO  2311  D  1S 

UO)«0 


2311 


1411 

C 

C 

c 

151 

161 


1913 


1711 


IF<RQUANT  GT  1>W< I >«W< I >/RQUANT 
IF  <W<I>  GT  TL1)  GO  TO  131 
GO  TO  241 

IF  <W<  I-l>.  LE.  TL1  >  GO  TO  141 

GO  TO  151 

IEM-IEM+1 

IEA-IEM+NEMC 

IF<  IEM.  LE.  IEMAX>GO  TO  1411 
WR I  TE  <  3,  1412  >  I  EMAX,  K 

FORMAT < IX.  39HEVENT  COUNTER  EXCEEDED  MAX  VALUE,  1 

IEM=IEMAX 

IC1<IEA>«I-1 

IC2<IEA>-0 

PEAK  DETECTION,  LOCATE  AND  COUNT  PEAKS 

IF  <W<I>-W<I-1>>  171,161,161 

IPB=I-i 

GO  TO  181 

IF  CIPB  EQ  0>  GO  TO  181 
IP- IF+1 

IF(IP  LE  JMAX>GO  TO  1711 

WRITE <3, 1913 > IP,  I EVENT 

FORMAT < IX,  17HN  PEAKS  EXCEEDED,  2I6> 

IP-JMAX 
GO  TO  181 

1 PRNG  <  lP>-<  I  +  IPBV2 

IPB=0 

CONTINUE 

GO  TO  282 

IF  <W<I-1>.  LE  TL1 >  GO  TO  281 
IC2<IEA)-I-1 


16,  5X,  I4> 


KEEP  COUNT  OF  PEAKS  WITH  EVENT 


2412 


C 

C 

1C 

C 

c 

C  1 
C  881 


IF  < IPB  EQ  0 >  GO  TO  251 
IP-IP+1 

IFUP  LE  JMAX>GO  TO  242 
WRITEC6,  1913 > IP,  IEM 
IP-JMAX 
GO  TO  243 

I  PRNG <  IP)-< I  + IPB>/2 
IPB-0 

I DC  < I EM >  -  l  P 

IF<W< I >  LE  TL2>  GO  TO  2412 
IF <W< I-l>  GT  TL2 >  GO  TO  281 
IEM2-IEM2H 

I F  < I EM2  GT  IEMAXJIEM2-IEMAX 
IC21C IEM2>-I-1 
IC22C I EM2  > “0 
GO  TO  281 

IF <W< 1-1)  LE  TL2>  GO  TO  281 
IC22< IEM2  >- 1 -1 
CONTINUE 
CALL  PEAKD 

STORE  PRESENT  PARAMETERS  IN  PREVIOUS  PARAMETERS 
IF <  IEL.  NE.  1>  GO  TO  802 

WRITEC6, 1>  AZMUTH, TL2, IEM2, <IC21<J>, IC22<J>,  J-l, IEM2> 
FORMAT <F8.  2,  218.  /11<2I6,  3X>/il<2I6,  3X>> 

WRITE<8>  AZMUTH,  TL1,  IEM,  <IC1<J>,  IC2<J>,  J-NEMi,  IER> 

JEM- I EM 

RETURN  , , 

END  ,Mft 


oooo  I  onooo 


07/19/79  12:54: 07 

♦♦♦LISTING  FOR  ERT1 :  PERKD.  FTN 

*N 

SUBROUTINE  PERKD 

If  ♦♦♦♦♦♦♦+♦♦♦♦  <♦*♦*♦♦  ♦+♦♦♦*•  If  >f  if*>f  ■+ -flf  if  iflflf  iflf  if  if  if  if  if  if  if  if  if  f  >f  if  if  if  if  if  if  *  if  if  f  f  if  iff  faff 

VERSION  5.  0  LEVEL  7S0616 
JHW  RVCO  IBM260 

DETERMINES  PERK  VALUES  RND  THEIR  ATTRIBUTES. 

*♦♦  If  If  If  If  )f  if  If  if  If  If  if*  If  *  *  ♦  If  +  if  If  If  If  If  >flfif  if  If  if  iff*  :f  if  if  if  if  If  if  if  If  If  if:  if  If  if  If  if  If  If  If  >f  if  if:  If  if  if  if  If  if  if  If  If 

IMPLICIT  INTEGERS  <I-N> 

REAL  UP<6>,  TRTR<1400>,  BUFC8) 

INTEGERS  RSAVE<16> 

INTEGER  *2  TCVL, TBVL, TCVLB,  TCVM, TRTC,  TRTM 
INTEGER  ^2  HB,  IRCT<70>,  IDC<22>, 

+  I PCRNG  <  34  > ,  TM 

INTEGER  *2  W,  RQURNT 
INTEGER+2  TL2,  T,  TC 

COMMON  /DRTR1/  ECL<224>,  NOFST,  KOFST,  I  CL  AD,  NANI 
dOMMON  /SWITCH/  IC1<44>,  IC2<44>,  TC<1980>,  IPC1<5400>, 

+  IFC2<5400>, IPC3<5400>, IPCNT<1980>, IPTC(44>, 

+  NEMB,  NEMC,  NRC 

COMMON  /PWORK/  KMRX,  T<80>,  JMXDB, JMAX,  IRMRX,  IR,  JR,  IMXJMX 
COMMON/REFL/  WC1025>,  HBC258>,  NCL,  NIO, NIDP,  INCL, 

X  I  MX,  IMN,  TM,  TL2,  RQURNT,  I DVEL  C  258  > 

COMMON  /F I XED/NFR,  IEMRX,  NFC,  I EM,  JEM 
COMMON/RZM/  RZMUTH,  R2LRST ,  RZSTRR,  NR,  ELCVRT 
COMMON  /RZ2/SRZ,  CRZ,  DRZ,  ISCfiNF,  NEL,  RI,  SR 
COMMON/F I LTER/TflTRMN,  RRERMN,  DR2M 
COMMON  /THRESH/  LDB 

COMMON  /PRSTOR/  NUP, TRTR(1400>,  NUMRX,  IRCT<70>, 

♦  I DC  <  22  > ,  I PCRNG  <347, KRMRX, MXTR 

COMMON  /CONPK/  NCELL,  NMX 
DRTR  DPR/57.  296/,  A/0.  1/ 

I  EM  IS  NO.  OF  EVENTS  IN  C  RRDIRL. 

INITIALIZE  RND  GENERATE  HC  RRRRV 

NBADR^NCRDR 
NCADR= I RMRX^NRC 
NBKR=NCKR 
NCKfi=KRMRX*NAC 
NRX=NR 
ITV*0 

IF  CNR.  NE.  1>  GO  TO  2109 
•"  LM=6 

LMDP«LM^NIDP 
NBADR^ I RMRX+NRC 
N6KR=KAMRX^NAC 
NRN-0 
NRN1-0 

LMM-LM-1 -  " 

IDX=LM+1 
NIMN»1 

FQUflNT  *RG)URNT 
NCLM*NCL-1 
LDBM=LDB-1 
-  NUMP=2+LM+LDB 
LDX*NUMP-1 
NPDP*LDX^NIDP 
I D2*l*  <  LDB~1 > ♦LM 
DO  17  1*1, KOFST 
17  UP<I>*0. 


1>- 17 


DO  18  1=1,  LIMT 

18  ECL<I>=0. 

DO  19  1=1, NCL 

19  IDVELC I >=0 
C 

NMX=1 

DO  2107  1=1,  NIDP 
2107  I  ACT < I >“0 

DO  2108  J=l,  MXTR 

2108  TRTRCJ>=0. 

2109  NGM=0 

DO  23  K=l,  KMAX 
KAD=K>NCKfi 
23  I PONT  <  KRD ) *0 
C 

1044  IFCIEM.  LE.  0>GO  TO  952 
DO  951  IE=1, I EM 
I  ER=  I  ll+NEMC 
KIE=< IE-1>*KMAX+NCKA 
K I EM=K I E-KMRX 
IPTCC1EA>=0 
ICEST=IC1CIER> 

IC£SP=IC2C IEfl> 

IFC IE.  EQ.  l)GO  TO  232 
DO  233  K«l, KMRX 
KR=K+KIE 
KB=k>KIEM 

233  IPCNT  <KR>=IPCNT  <KB> 

232  IPL=0 

IF<IE.  GT.  1> IPL=IDC< IE-1> 

IP=IDC< IE) 

IF<  IP.  LE.  I  PL  )  GO  TO  951 

IPL=IPL+1 

JE1=0 

JE2=0 

C  FIND  B  EVENTS  ASSOCIATED  WITH  C  EVENTS. 

C  .’EM  IS  NO.  OF  EVENTS  IN  PREVIOUS  RRDIRL. 

C 

IF< JEM  EG'  0>  GO  TO  41 
DO  31  JE=1, JEM 
JEB=JE+NEMB 

IF< IC2C JEB>  .  LT.  ICEST>  GO  TO  31 
IF <  IC1C  JEB>.  GT.  ICESP>  GO  TO  41 
JE2=JE 

IF  C  JE1.  EQ.  0>  JE1=JE 
31  CONTINUE 

C  FIND  THRESHOLDS  FOR  IE  EVENT 

C 

41  DO  51  J=l, JMXDB 

51  T  <  J  >  =0 

NTHRES*1 
DO  71  L-IPL,  IP 
IF  CL.  GT.  JMRX>GO  TO  71 
I Rl* I PCRNG  <  L  > 

IFCIR1.  LT.  ICEST>GO  TO  71 
IFCIRl.  GT.  ICESP>GO  TO  712 
DO  711  K=l,  LDB 
IU=WC IR1> 

IT-IU-TM-K+1 
IFCIT.  LE.  8>G0  TO  711 
IFCIT.  GT.  JMXDB )  I T = JMXDB 
IF<T<IT>.  EQ.  0>NTHRES=NTHRES+1 
T  < IT)=1 


M  711  CONTINUE 
W  71  CONTINUE 


D-18 


IFCNTHRES  GT.  KMAX>WRITE<6,  787SJNA,  NTHRES. KMAX.  IE 

7878  FORMAT < 2X,  33HNUMBER  OF  THRESHOLDS  EXCEEDS  KMAX.  4Ilfl> 

IF<NTHRES.  GT  KMAX  >  IPSRT<*NTHRES-KMAX 

IPT-1 

DO  91  L-l.  JMXDB 
IF  <  T<L  >.  LE.  0>  GO  TO  91 
KA-IPT+KIE 
TC  <  KA  > “L+TM-1 
IPSRT-IPSRT-1 
IF<  IPSRT.  GT.  0)GO  TO  91 
IPT-IPT+1 
91  CONTINUE 
IPT-IPT-1 

IFCIPT  GT.  JRM<JRITE<6<  7879>NA.  IPT.  JR.  IE 

7879  FORMAT < 2X.  31HNUMBER  OF  THRESHOLDS  EXCEEDS  JR. 4I10> 

IF<  IPT.  GE.  JR  MPT*  JR 

IPTC<IEA>*IPT 

IF< IPT  LE  0>GO  TO  951 

C  LOOP  ON  RANGE  IN  IE  EVENT  TO  FIND  CONTOUR 

C 

IBGN-ICEST+i 
IND-lCESF+l 
DO  161  I  =*  I  BGN.<  IND 
C 

C  LOOP  ON  THRESHOLD 

C 

DO  131  K«i.  IPT 
IU*W<I> 

KA»K+KIE 

IF  <1U  GT.  TCCKAM  GO  TO  ill 
GO  TO  141 
111  IU«W<I-1> 

IF  <IU  LE  TC<KA>  >  GO  TO  121 
GO  TO  131 
C 

C  START  RANGE  FOR  SEGMENT  <CONTOUR> 

C 

121  IPCNT  <KR>*IPCNT  <KA>+1 

IF<IPCNT<KR>  LE  IMXJMXJGO  TO  1211 
MRITE<6. 1212 > I TV. K< IE 

1212  FORMAT < 2X.  3GHNUMBER  OF  SEGMENTS  EXCEEDS  I MX.  3110) 

I PCNT <  K A  >  * I MX JMX 
1211  I PE* I PCNT  <KA> 

I REG* I -1 

I RDDR- 1 PE+  <  K- 1 >  * 1 MX J MX+NCADR 
I PCI < I ADDR  >  * I REG 
I PC3  < I ADDR  >  *0 
131  CONTINUE 
GO  TO  161 
C 

C  END  RANGE  FOR  SEGMENT 

C 

141  DO  151  KL-K. IPT 

IF  <WCI-1>.  EQ  -999>  GO  TO  161 
IU-W<I-1> 

KA-KL+KIE 

IF  <IU  LE  TC<KA>  >  GO  TO  161 
IPE»1PCNT<KA> 

I REG* I -1 

IEQL«IPEMKL-1>*IMXJMX+NCRDR 
IPC2< IEOL  >* IREG 
151  CONTINUE 

161  CONTINUE  p 

C 

C  ASSOCIATE  CELLS  LOOP  ON  THRESHOLD  HIGHEST  TO  LOWEST 


940  DO  941  LC*1,  IPT 
KC-IPT-LC+1 
KR-KC+KIE 
KZ-KC+KIEM 
IF<KC.  LE.  0>GO  TO  941 
TCVL=TCCKR> 

TCVM*TCVL+1 
TCVLB*TCVL+LDB 
NFC* IPCNT  CKR) 

NPL-0 

IFC1E  GT  1 ) NFL  * I PCNT  C  KZ  > 

IF CNPC  LE.  NFL) GO  TO  941 
NPL-NPL+1 

C  LOOP  ON  SEGMENTS 

DO  931  IPE-NPL,  NPC 
JRDDR* I PE+  C  KC- 1 > ♦ I MX JMX+NCRDR 
IHBM=IPClv.JRDDR> 

IHB=IHBM+1 
IHD=IPC2C JRDDR) 

K-KC+1 

KV-KR+1 

KX*KZ+1 

NFK=>0 

TRTM=0 

LPE= I PCNT  C  KV  > 

LPL*0 

IFC  IE  GT  1>LPL=IPCNTCKX> 

LPL*LPL+1 

IF<LPE.  LT  LPL.  OR  K  GT  IPT)GO  TO  193 

C 

C  LOOP  SEGMENTS,  NEXT  HIGHER  THRESHOLD 

C 

192  DO  191  L*LPL,  LPE 

I RDDR=L+ C  K- 1 ?  * I MX JMX+NCRDR 
IFCIPC2v IRDDR?  LT  IHBM>GO  TO  191 
IFCIPCIURDDRX  GT.  IHD)GO  TO  193 
NFCEL* I PC3  Cl RDDR ? 

IFCNPCEL  LE  0>GO  TO  1911 
TRTC=TRTRCNPCEL> 

IFCTRTC  LT.  TRTM)  GO  TO  231 

TRTM*TRTC 

NFK=NPCEL 

231  IFCTRTC  LT.  0>TRTC*-TRTC 
C 

C  NPCEL  IS  FOR  NEXT  HIGHER  C ENCLOSED?  THRESHOLD  ON  C  RRDIRL 

C 

IFCTRTC  GT.  TCVLBJGO  TO  932 
191  CONTINUE 
GO  TO  193 
932  NPK=-NPCEl 
GO  TO  193 

1911  NPK— CNIDF+1? 

C 

C  RSSOCIRTE  CELLS  ON  B  RRDIRL,  TOP  DOWN 

C 

193  MPK-0 

IF CNR  EQ  1>G0  TO  361 
TRTM-0 

IF  CJE2  EQ.  0)  GO  TO  371 
DO  261  JE-JE1,  JE2 
JEB-JE+NEMB 
I TRTM-0 

IFC IC2C JEB>.  LT  IHBM)  QO  TO  261 
IFCICKJEB)  GT  IHD>  GO  TO  3661  „  ,, 


I  ««  IV  »  Mi.  4  I  l.'-'UV  A  I  I  I  L.V 


c 

f 

e 

t 


c 

271  IPB«IPTCCJEB> 

1FUPB.  LE.  0>OO  TO  261 

DO  291  LB-1,  IPB 

KB-IPB-LB+1 

KBB- < JE-1 > *KMRX+NBKR 

KBR-KB+KBB 

KBC-KB+KBB-KMRX 

TBVL-TC<KBR> 

JEQL-TBVL+1 

NP2-IPCNT<KBR> 

NP1-0 

IFCJE.  GT  1>NP1— IPCNT <KBC> 

IF<NP2  LE  NP1 >00  TO  291 
NP1-NP1+1 

DO  281  JPE-NP1,  NP2 
I RDDR« JPE-K  KB- 1 > ♦ I MX JMX+NBRDR 
IF< IPC2C IRDDRV  LT  IHBM>GO  TO  281 
I F  <  I PCI v I RQDR > .  QT  IHD>GO  TO  2911 
LPCEL-IPC3C I RDDR  > 

IFv'LPCEL  LE  0>GO  TO  281 
IF<  TCVL  LE.  TBVL >  GO  TO  282 
I EQL- TRTR  <  LPCEL  > 

I F  <  JEQL  LT  I EQL > GO  TO  281 
282  TRTC-TRTR  <  LPCEL  > 

IFCTRTC  LT  TRTM>  GO  TO  281 

TRTM-TRTC 

1TRTM-TRTM 

MFK-LPCEL 

KBM-KB 

JBM-JE 

281  CONTINUE 

2911  IF<  I TRTM.  NE  O’*  GO  TO  261 
291  CONTINUE 
261  CONTINUE 
3661  IF<MPK  EQ  0>GO  TO  371 

IFCTRTM  GT  TCVLB>MPK=-MPK 
GO  TO  421 

371  DO  194  I-IHB,  I HD 

I F  <  HB  < I > .  EQ  -999 > GO  TO  194 
IF< IRBSvHBi 1 >  >  LE  TC<KR>  >G0  TO  194 
IF (NPK  EQ  0>GO  TO  931 

IF <NPK  GT  0>G0  TO  366 

GO  TO  3662 
194  CONTINUE 
C 

C  HRVE  B  COMPRRE  WITHIN  RRNGE 

C 

361  CONTINUE 

IF <  NPK  EQ  0>GO  TO  631 
C 

C  MPK-0  RND  NPK-C  -  NO  COMFRRE 

C  MPK-0  RND  NPK  NE  0  -  NO  B  COMPRRE 

C  NPK-B  RND  MPK.  NE.  B  -  B  COMPRRE 

C  HIGHEST  THIS  RRDIRL 

C 


£  392 


IF<NF*K  LE  0  OR  NPK  GT  NMX>GO  TO  3662 

NO  PRIOR  RRDIRL  FOR  COMPRRISON,  INCREMENT  NPCEL 


NPCEL-NPK 

I NDX- T RTR C NPCEL > - TC VM 
IF < INOX.  GE  LDB  OR  INDX  LE  0>GO  TO  366 
rN-INDX*Lmi 
I NX- I DX+ 1 NDX*LM 


p- :  i 


«  ■*.->«<«  ft  W'l 


lEQL-TRTRvJN' 

IFvlEOL  N£  0  OR  Nfl  EG  l 'GO  TO  3921 
1EQL '■'■NPCEL*  v  IN-LNM-l  '  *NIOF 
JEGL-TFTRv 1EGL ' 

IFvJEGL  LE  O'GO  TO  366 
NPC -NPCEL 
NPCEL-JEQL 

IFvMPC  EG  NPCEL  OR  NPCEL  GT  NMX > GO  TO  366 
GO  TO  339 

39C1  IPC3vJHOOR>-NPCEL 
JN1-JN+NIOP 
JNi-JUl^NIOP 
JN3-JN2*NIDP 
JN4«JN3*NIOP 
JN3-JN4+NI0P 

IF v THTRv  JN1 '  EO  -999  >GO  TO  419 
1ST-1NB 
I SR* I HO 

OO  411  l* 1ST  ISP 
R-RI*SH*vFLOHTvI  1'-  3' 

IU-WV I > 

KU»R*FLOHTv  lU'*OHr 
THTRv JNi'-THIRv JN1 ' *OHT*R 
THTRv  JN-2'«*THTRv JN*'  +  RU 
THTRv  JN3 ' -THTRv  JN3 ' ♦SfC*R*RU 
THTRv  JN4 '“•THTRv  JN4  >*CfC*R*RU 
THTRv  JN3 '“•THTRv  JN3'*  I OVELv  I  > 

411  CONTINUE 
419  KN»*NPCEL*  v  INN-1  '*NlOP 
THTRvkN'-NH 
kNN-NPCEL ♦NPl'F 
THTRvkNN'«IE 

I F  v  I  ST  EQ  2  OR  ISF  EQ  IMN'THTRv  JN1  '  —  999 
GO  TO  366 
366*  NPCEL«-NPk 

3*6  I F  v  NPCEL  GT  NMX  OR  NPCEL  IF  O' GO  TO  931 
I NON- THTR v  NPCEL ' - TCVM 
C 

C  COMBINE  LPCEL  WITH  NPCEL  HT  THIS  LEVEL 

C  COMBINE  BV  SETTING  HREH  HS  POINTER  HNO  ION  TO  NH  -  0 


IF  vLPE  LT  LPL  OR  k  GT  IPT'GO  TO  931 
OO  363  L-IPL  LPE 
l HOOR-L ♦ v k - 1 ' ♦ I MX JMX 
I F  v I PC*  v I MOOR  >  LT  I HBM  > GO  TO  363 
IF v IPClv I MOOR'  GT  IHO'GO  TO  931 
LPCEl  -I PC  V IHOOR ' 

IF vLPCEL  l  E  0  OR  LPCEL  GT  NMX'GO  TO  363 
LNX-LPCEL ♦ LMOF 
IF v THTRv LNX'  EO  0  >GO  TO  363 
IF v NPCEL  EQ  LFHEL'GO  TO  363 
I NDX *  THTR  v  LPCEl  '-TCVM 
IFvlNOX  GE  L  OB  GO  TO  363 
IF v INOX  LE  O'lNDX-O 
INO-LPCEL+v INOX*l '♦LMOP 
IFv  THTRV  INI>>  EQ  0  'GO  TO  363 
I NO- 1 NON -M 
IPG-0 

OO  3663  J - 1 NO -  LOB 
IN**<  J-1>*LM 
I EGL-LF'CEL* J*LMOP 
JEQL-THTRv IEDL ' 

IFvJEQL  EQ  NH' IPG-IPG+i 
OO  3663  I-l.LM 
IEQL-LRCCL*V  I*IN'*N1DP 


0- 


3666 

3667" 

3665 

3664 


365 


EG  >J  OR  IE  LE  l'GO  TO  3664 
DO  3665  l-t,  IE 
I©"I+N£MC 
IPTT-IPTCU©> 

IFvIPTT  LE  ©'GO  TO  3665 
DO  3666  KT«i. IPTT 
KT I »  v  I  -1  > *kM©X>NCk© 

KTfl«kT4kTI 

Kre*kT^KTi~krtax 

NPCT«IPCNT<kT©' 

IEGL*TC<;kT©>+l 

NPCL-© 

IFCI  GT  i  'NPCL*  IPCNT  <,kTB' 

IFvNPCT  LE  NPCL'GO  TO  3666 

NPCL=NPCL*i 

DO  3667  lp-npcl.npct 

I©DDRl=LP-*'<«kT-l  '♦  INXJMN*NC©OR 

IF  ..LPCEL  NE  IPC3>.  I©PDR"GO  TO  3667 

I  NON  T -T6TR «.  NPCEL  '  - 1 EOL 

IFvINOST  LT  LOB'GO  TO  3666 

IPC3v I©OOR'-0 

GO  TO  3667 

IF < l MON  GE  LOG 'GO  TO  366? 

I PC  3  v I ©OOR ' *NPCEL 

CONTINUE 

CONTINUE 

CONTINUE 

IPG-O 

IFvINOX  GE  LOG 'GO  TO  365 

IrtCTv LPCEL '"-NPCEL 

IEQL-LPCEL+v INDN^LM+1 '♦NIOP 

TrtTRv I EOL ' "NPCEL 

IF v INOX  N€  ©'GO  TO  365 

I  ACT  *.  LPCEL  >  *-NIOP-l 

LN"LPCEL*NIPF 

T©TR«.LN'~0 

CONTINUE 

GO  TO  931 


C 

C 

c 

c 

c 

c 

4*1 


COMBINE  NPCEL  ©NO  LPCEL  FE©k  VALUES  EOUftL 


COMBINE  WITH  F  RHDIftL  CELLS 

IF^MPk  LE  ©'GO  TO  4^2 
IF <NPk  LT  ©'GO  TO  366T 
NGM*© 

LPCEL ’MPk 

LNN«LPCEL*LNOF 

IECJL»THTR<LNX' 

KBN©»kBM*  v  JBN-l '  *kM©N>NBk© 

IFvIEOL  EQ  N©  ©NO  NPk  EO  ©  ©NO  TCVL  GT 
♦  TCvkBN©"  GO  TO  465 
I NON-TRTR< LPCEL ' -TCVM 
IMDX-INOX 

IFvNPk  GT  ©'IMOX*T©TR«,NPk'-TCVT1 

IF*. IMDN  LE  INDX'GO  TO  4^12 

NGM*1 

NPCEL -NPk 

INO-INOX 

INOX-IMCX 

INOX*- 1  NO 

GO  TO  4*13 

IF  *  INOX  LT  ©'GO  TO  46M 
NPCEL -LPCFL 


c 

c 

4212 


512 


5214 


5212 


5212 


5311 


531 

8012 


C 

C 

£T - 

C 

C 

C 

481 


COMBINE  WITH  B  -  RADIAL,  C-LEVEL  LOWER 

IFC  INDX.  GE.  LDB>GO  TO  4221 
IN=INDX*LM 
IN1=C IN+1)*NIDP 
ILN=CIN+LM>*NIDP 
I EQL=NPCEL+ I LN 

IFCTATRC IEQL).  NE.  0.  >GO  TO  5311 
I EQL=NPCEL+ I N1 
JEQL=TATRC IEQL> 

IF C  JEQL.  LE.  0.  AND.  NGM.  EQ.  0>GO  TO  4221 

IFCNGM.  NE.  l)GO  TO  5312 

IM=IMDX*LM 

IFCIM.LT.  0>  GO  TO  5211 

IF CLPCEL.  LE.  0  OR.  LPCEL.  GT.  NMX)  GO  TO  422 

I EQL=LPCEL+  C I M+LM ) *N I DP 

IFCTATRC  IEQL>.  NE.  0.  )GO  TO  5311 

IEQL=LPCEL+C IM+1)*NIDP 

IFCTATRC  IEQL >.  GT.  0.  >GO  TO  5313 

LPCEL=NPCEL 

GO  TO  4221 

LPCEL=TRTRCIEQL> 

IFCLPCEL.  EQ.  NPCEL.  OR.  LPCEL.  GT.  NMX ) GO  TO  4221 

I MDX=T ATR  C  LPCEL ) -TCVM 

GO  TO  5314 

IEQL=NPCEL+IN1 

NPCEL=TATRCIEQL) 

IFCNPCEL.  LE.  0.  OR.  NPCEL.  GT.  NMX) GO  TO  4221 

I NDX=T ATR  C  NPCEL  ) -TCVM 

GO  TO  4213 

I PC3 C  JADDR >  =NPCEL 

IEGL=NPCEL+IN1 

IFCTATRC  IEQL).  EQ.  -999.  )  GO  TO  8012 
I ST —  I HB 
I SP= I HD 

I EQL2= I EQL+N I DP 
I EQL3= I EQL2+N I DP 
I EQL4= I EQL2+N I DP 
I EQL5= I EQL4+N I DP 
DO  521  I  =  IST,  ISP 
IU=WCI> 

R=R I +SR*  C  FLOAT  C I -1 ) — .  5) 

RU=R*FLORT  C I U ) *DAZ 
TATRC IEQL)=TRTRC IEQL)+DAZ*R 
T ATR  C I EQL2  )  =T ATR  C I EQL2  )  +RU 
TATRC I EQL3 )  =TRTR  C I EQL3 ) +5 AZ*R*RU 
T ATR  C I EQL4 ) =T ATR  C I EQL4 )  +CRZ*R*RU 
T ATR  C I EQL5 ) =T ATR  C I EQL5 )  + 1 DVEL C I ) 

CONTINUE 
I EQL=NPCEL+ I LN 
T ATR  C I EQL  )  =NR 
NMP=NPCEL+NPDP 
TATRCNMP)=IE 
IEQL=NPCEL+IN1 

IFCIST.  EQ.  2.  OR.  ISP.  EQ.  I  MX  >  TATRC  I  EQL  >=-999. 

LPCEL=NPCEL 
GO  TO  4221 

COMBINE  WITH  B-RADIAL,  C-LEVEL  HIGHER 

IF  FIRST  COMBINE,  AREA=0,  IF  SECOND  OR  HIGHER,  RREA=-1. 
TEST  AREA  TO  ESTABLISH  NEW  NUMBERS 

INDX=-INDX 
IND-NUMP-1 


D-24 


i  l<VT  -L-C'C> 

INS“2 

IPG=0 

TRTRvLPCEL  >=TCVM 

LNP=LPCEL+NPDP 

TRTRvLMP>=IE 

IFvINOX  GE  LDEOGO  TO  4S2 
IN0-LD6-INDX 
DO  4832  I = INOX, LDBM 
IEQL=LPCEL-K  I+1HLMDP 
JEQL-TRTRc I EQL  > 

IF  JEQL  EQ  NO  >  I  PG“  I  PG+1 
4832  CONT I NUE 

DO  483  1=1,  IND 
DO  483  J=i,LN 

IN=LPCEL-KJ-KLDB-I  '"♦‘LM>'*'NIDP 
IM=LPCEL+<  IND-I  >*LMWNIDP 

483  TRTRv.  IN>=TRTRv.  IM> 

I ND= I NDX*LN+1 
INDP=INDN 

482  DO  4S35  I “1, LDB 

IEQL=LPCEL+ I *LMDF 
JEQL=TRTR< I EQL  > 

IF*.  JEQL  EQ  NR> IPG=IPG+1 
4835  CONTINUE 

DO  4S4  I = I NS, IND 
IN-LPCEL-K  I-i>+NIDF 

484  TRTRt  IN>“0 

DO  4841  1  =  1,  INDP 
I EQL=LPCEL+ I ♦LMDP 
4841  TRTRv  IEQL''=Nfl 

IFUPG  EQ  O  OR  IE  LE  i'GO  TO  4SS 
DO  4831  1=1,  IE 
IR=I+NEMC 
IPTT=IPTC< Ifl> 

IF «.  IPTT  LE  0>GO  TO  4S31 
DO  4833  KT  =  1,  IPTT 
KTP=>,  I -1 '"♦■KMRX+NCKfl 
KTR=KT+KTP 
KTB=KT +KTP-KNRX 
NPCT  =  IPCNT<KTR'> 
lEQL=TCa.TR>  +  l 
NPCL=0 

IF  v  I  GT  l>NPCL=lPCNTO;'re'' 

IFvNPCT.  LE  NPCL^GO  TO  4833 
NPCL-NPCL+1 
DO  4834  LP=NPCL.  NPCT 
KRDDR=LP-K  KT -1 > +  1 MX JMN+NCRDR 
IFCLPCEL  NE  IPC3vKRDDR'»  >GO  TO  4S34 
INDXT=TRTR<,LPCEL>-IEQL 
IFvINDXT  LT  LDB  ''GO  TO  4S34 
I PC  3  <  KRDDR  >  =0 
4834  CONTINUE 
4S33  CONTINUE 
4831  CONTINUE 
IPG=0 
488  IN=0 

IFCLPCEL  LE  0  OR  LPCEL  GT  NMX>GO  TO  931 

LND=LPCEL+Lf1DP 

TRTR<LND>“NR 

I PC3  <  JRDOR  ? “LPCEL 

NPCEL-LPCEL 

NGM=0 

GO  TO  512 

485  DO  488  I«NIMN. NIDP  p-:S 

IF<  IRCT «.  I  V  EQ  O >GO  TO  4S" 


m 


i 


WRITE(6,  644 > 

GO  TO  931 
487  LPCEL=I 
N I MN= I  ♦  1 
IRCT  ( I >=1 

IFCNIMN  GT  NIDP>  N I MN=N I DP 
IF(NMX  LT  NIMN)  NMX=NIMN 
TRTR( LPCEL >-TCVM 


C 

C 

C 

502 

851 

852 

8511 

8512 

861 


LNP=LPCEL+NPDP 
TRTR(LNP>=IE 
GO  TO  4S8 
LPCEL=MPK 

IF(LPCEL  LT  0>LPCEL=-LPCEL 

IF (LPCEL  GT  NMX.  OR  LPCEL  LE  0>GO  TO  2662 

I MDX=TRTR ( LPCEL  > - TCVM 

IFCIMDX  LT.  0>GO  TO  632 

DO  441  JE=JE1, JE2 

JEB= JE+NEMB 

IF(  IC2( JEBX  LT  I HBM >  GO  TO  441 
IF< I Cl (JEBX  GT  IHD>  GO  TO  632 
I PB= I PT C  (  JEB  > 

IF( IPB  LE.  0>GO  TO  441 

DO  471  LB=1, IPB 

KB=IPB-LB+1 

KB J  = ( JE-1 > *KMRX+NBKR 

K6R=JE+K6J 

KBC= JE+KBJ-KMRX 

IF(TC(KBRX  NE  TCVL)  GO  TO  471 

MPB= I PCNT  <  KBR  > 

MPL=0 

IFCJE  GT  1>MPL=IPCNT(KBC> 

IF(MP6  LE  MPL  > GO  TO  471 
MPL=MPL+1 

DO  461  JFE=MPL, MPB 
I RDDR= JPE+ (  KB-1  >  + 1 MXJMX+NBRDR 
IF( IPC2( IRDDRX  LT.  I  HBM  >  GO  TO  461 
IF ( IPC1( IRDDRX  GT  IHD)GO  TO  471 
NPCEL=IPC3( I RDDR  > 

I F < NF'CEL  LE  0  OR  NPCEL  GT  NMX) GO  TO  461 
IF (LPCEL  EQ  NPCEL >GO  TO  461 


COMBINE  RT  TB=TC  LEVEL 

I NDX=TRTR (  NPCEL > -TCVM 
I F  ( I NDX.  GE  LDB )GO  TO  461 
IFONDX.  LT.  ©'•GO  TO  8511 
IFCIMDX  LT.  LDB > GO  TO  861 
ND=INDX*LM+1 
DO  852  1=1, LM 
I EQL=NPCEL+  < ND+ I -1 >  *N I DP 
TRTR< IEQL>=0 
IPC3( I RDDR ) =0 
GO  TO  461 
IPC3( IRDDR>=0 
DO  8512  J«l,  LDX 
JN=NPCEL+ J*N I DP 
TRTR<  JN>=0. 

IRCT (NPCEL >=-(NIDP+l> 

GO  TO  461 
LD=IMDX*LM+1 
ND=INDX*LM+1 
LDR=LD*NIDP 
NDR=ND*N I DP 

I EQL-LPCEL+  <  LD+LMM  >  *N I DP 
JEQL-NPCEL+  <  ND+LMM ) *N I DP 


5 


I  M  I  Kl_=  ItllKk  itUL,1 

TRTRN=TATRCJEQL> 

IFCTRTRL,  NE.  0.  .  RND.  TATRN  NE.  0.  >00  TO  ©911 

JEQL=LPCEL+LDA 

TRTRJ=TRTR< JEQL> 

IFCTRTRL.  EQ.  0.  .  RND.  TRTRJ.  LE.  0.  >G0  TO  851 

IFCTRTRL,  GT.  0  >G0  TO  8912 

LPCEL-TRTRJ 

IFCLPCEL.  LE.  0.  OR.  LPCEL.  GT.  NMX>GO  TO  461 
GO  TO  4221 
8912  IEQL=NPCEL+NDA 

TRTRJ=TRTR< IEQL> 

IFCTRTRN.  EQ.  0  .  RND.  TRTRJ.  LE.  0  >G0  TO  8913 
NF'CEL=TRTRJ 

IFCNPCEL.  LE.  0.  OR  NPCEL.  GT.  NMX.  OR  NPCEL.  EQ.  LPCEL>GO  TO  461 
I PCS  < I RDDR  >  =NPCEL 
GO  TO  502 

8913  DO  8914  1=1, LM 

I EQL=LPCEL+  C  LD+ 1 - 1 >  *N I DP 

8914  TRTRCIEQL>=0 

I PC3  < I RDDR  > =0 
GO  TO  4221 
8911  I6NDRV=0 

IEQL=LPCEL+LDA 

JEQL=NPCEL+NDR 

IFCTRTRC  IEQLX  EQ.  -999.  .  OR.  TATRCJEQL).  EQ.  -999.  > 

X  I BNDRV=1 
DO  891  1=1,  LMM 
I EQL=LPCEL+  <  LD+ 1 -1 >  *N I DP 
JEQL =NPCEL+  <  ND+ I -1 >  *N I DP 

IFCIBNDRY.  EQ.  0>TRTRC IEQL >=TATR< JEQL  '+TRTRC IEQL> 

TRTRC JEQL>=0. 

891  CONTINUE 

IEQL=LPCEL+LDfl 

I F  C  I BNDRV.  EQ.  1>TATR<  IEQL>=-999. 

I EQL=NPCEL+ <  ND+LMM  >  *N I DP 
TRTRC  IEQL>=0. 

I EQL=NPCEL+NDR 
TRTR  < I EQL  >  =LPCEL 
I ACT  <  NPCEL  > =-LPCEL 
IPC3C I RDDR  ?  =LPCEL 
461  CONTINUE 
471  CONTINUE 
441  CONTINUE 

632  IFCNPK.  LE  0>GO  TO  3662 
NPCEL=LF'CEL 
GO  TO  366 
C 

C  UNRSSOC I RTED 

C 

631  DO  642  J=NIMN,  NIDP 

IFC IRCT C  J).  EQ  0>GO  TO  643 

642  CONTINUE 

WR I TE  C  6, 644  > 

644  FORMAT CSX.  15H  TOO  MRNV  CELLS > 

GO  TO  931 

643  NPCEL-J 
NIMN=J 

IFCNIMN.  GT  NIDP>NIMN=NIDP 
IFCNMX.  LT.  N I  MN  > NMX=N I MN 
IRCT  C  J)=l 
IPC3C JRDDR > -NPCEL 
DO  671  1=1,  NUMP 
NMP=NPCEL+CI-1>*NIDP  . 

TRTRCNMP:>=0  0 
671  CONTINUE 


NNN-NPCEL+NPDP 
TRTR<NMN)-IE 
IST-IHB 
ISP* I HD 

NP2=NPCEL+NIDP 

NP3=NP2+NIDP 

NP4-NP3+NIDP 

NP5-NP4+NIDP 

NP6=NP5+N I DP 

DO  621  I*IST,  ISP 

R=RI+SR*<FLOAT < I— 1)-.  5) 

IU»W<I> 

RU=R*FLOAT  < I U  >  *DRZ 

TRTR  <  NP2  >  «DAZ*R+TATR  <  NP2 > 

TRTR  <  NP3 >  *RU+TATR  <  NP3 > 

TRTR <  NP4 >  =SRZ*R*RU+TRTR  <  NP4 > 

TRTR <  NP5  >  =TRTR <  NP5 >  +CRZ*R*RU 
TRTR  <  NP6  >  =TRTR  <  NP6 >  + 1 DVE  L  <  I  > 

621  CONTINUE 

NIX=NPCEL+LMDP 

TRTR<NIX)=NA 

IF<  1ST.  EQ  2.  OR.  ISP  EQ.  IMX)TATR<NP2>=-999 
931  CONTINUE 
941  CONTINUE 
951  CONTINUE 
C 

C  CLERN  UP  TRTR  RND  IC  RRRRVS  -  REMOVE  IC  POINTER 

C  TO  DELETED  RRRRVS 

C 

7010  DO  9512  1=1, NMX 

IF<  I  ACT  <  I  >.  EQ  0>GO  TO  9512 

IF< IRCT < I ).  GE  0>GO  TO  9611 

DO  9613  JE=1, I EM 

IER=IE+NE  MC 

KIE=< IE-1) *KMRX+NCKA 

K I EM=K I E-KMAX 

IPT=IPTC<IEfl> 

IF< IPT.  LE  0>GO  TO  9613 

DO  9618  KC=1, IPT 

KR=KC+KIE 

KB=KC+KIEM 

NPC=IPCNT  <KR> 

TCVM=TC<KR>+1 

NPL»0 

IF<  IE.  GT.  1  > NPL=  I PCNT < KB > 

IFCNPC  LE.  NPL>GO  TO  9618 
NPL=NPL+1 

DO  9619  IPE=NPL, NPC 
I RDDR= I PE+  <  KC-1 >  * I MXJMX+NCADR 
IFCI.NE.  IPC3<  I ADDR)  )  GO  TO  9619 
IF<  IRCT  <  I  >.  LT.  -N I  DP)  GO  TO  9614 
INDX*TRTR<  I  >-TCVM 
IEQL=I+<INDX+1)*LMDP 
JEQL=I+<INDX*LM+1)*NIDP 
IF<TRTR< IEQL)  NE.  0  >GO  TO  9619 
IEQL-TATRCJEQL) 

IF<  IEQL.  NE.  -IRCT <  I  )  )GO  TO  9614 
IPC3CI  RDDR>  —  IRCT<  I ) 

GO  TO  9619 
9614  IPC3<IRDDR>-0 
9619  CONTINUE 
9618  CONTINUE 
9613  CONTINUE 

IF< IRCT < I >.  GE  -NIDP>GO  TO  9317  n  ?8 
IRCT  < I >"8 


TRTRU>-0. 

GO  TO  9512 

9517  DO  9512  J-l, LDB 

KEQL-I+<LM*< J-1>+1>*NIDP 
JEQL - TRTR  <  KEQL  > 

IEQL-I+J-LNDP 
IEQL-TRTR<  IEQL> 

IF< JEQL  EQ  -I ACT <  I  >.  AND  IEQL  EG  0> 

♦  GO  TO  9514 
9512  CONTINUE 

GO  TO  9611 
9514  TATR<KEQL>-0 
9611  DO  9612  K-2,  LDB 
IEQL- 1 -MOM-MOP 
JEQL- I  +  <  <  K-l >  *LM+1 >  *N I DP 

IF<TRTR< 1EQL>  NE  O  AND.  TATR< JEQL>  EQ  0  ) 

♦  TATR< IEQL>=0 
CONTINUE 
IRCTCD-1 
CONTINUE 

IF(Nfl  EQ  1 >GO  TO  1020 

END  OF  ASSOCIATION  LOOPS 

DO  991  1-1,  NMX 
I A— l  +  < LDX-1 >  *N I DP 
IFCIACTCI)  EQ  0?GO  TO  991 
IF<TRTR< IR>  EQ  0  >GO  TO  9912 
lEQL-TATRC IA> 

IF<  IEQL.  LE  OMEQL—IEQL 
IFCIEQL  EQ  NAX-1>G0  TO  971 
GO  TO  991 

CHECK  BACKROUND  CONING  DOWN 

INBR=0 
ITERM=1 

DO  9716  J-l, LDBN 
IEQL-I  +  C  C J~l>*LM+i>*NIDP 
JEQL*"IEQL+NIDP 

IF<TATR<IEQL>  LE  0  OR.  TATR<JFOL>  EQ  O  >GO  TO  99 
CONTINUE 
NMP- I +NPDP 
IEQL-TATRcNMP> 

J-l 

IF <  JEM.  LT.  2) GO  TO  96S 
DO  9711  J-l, JEM 
968  JA^J+NENB 

IF< IEQL  NE  J>GO  TO  9711 
9712  IPB=IPTC<JR> 

DO  9713  K-l,  IPB 
KR- <  J-l > -KMRX+NBKR 
KAP-K+KA 
KAM-K+KR-KMAX 
IEQL-TATR< I ) 

IEQL«IEQL-TC<KAP> 

IF< IEQL  NE  LDB>GO  TO  9712 
NP-1PCNT  <.KRP> 

NL-0 

IF< J.  GT.  1>NL— IPCNT <KAM> 

NL-NL+1 

DO  9712  N-NL,  NP 
I EQL-N+ <  K-l >  * I MX JMX+NBADR 
IFd.NE.  IPC2<  IEQL >  )GO  TO  9712 

INBR-INBR+1 

VTTt?M«0 


9612 

9512 

C 

C 

C 

952 

961 


C 

C 

C 

971 


9716 


irusL>K=n-r  «.K-a>*inXjnx+NbHDK 
IST=IPC1<IRDDR) 

I SP= I PC2  < I RDDR  )  +1 
DO  9715  L=IST/ ISP 
IU=W<L) 

IFCIU.  EQ.  -999 > GO  TO  9715 
IF<IABS<IU>.  GT.  TC<KAP))GO  TO  9982 
9715  CONTINUE 
9713  CONTINUE 
9711  CONTINUE 
ITERM=3 

IFCINBR.  EQ.  8>G0  TO  9982 
I TERM=4 

IID=I+ID2*NIDP 

IF<TRTR<  I  ID).  LE.  TRTRMN)GO  TO  9982 

DO  981  J=l,  LMM 

IEQL— I+< ID2+J-1)*NIDP 

981  UP< J)=TRTR< IEQL) 

UP  <  2  >  =UP  <  2  ) /UP  Cl) 

UP<6)=R*UP<2)+UP<5> 

DO  985  M=l,  NOFST 
MG=6+  <  M— 1 ) ♦KOFST+NRN1 
IFCUPCO.  CiT.  ECL<MG)  >GO  TO  986 

985  CONTINUE 
GO  TO  989 

986  NRN=NRN+1 
IFCNRN.  GT.  1  ) NRN=0 
NRB=NRN1 
NRN1=NRN* I CLRD 
LMT=NOFST-l 

DO  988  J=l,  LMT 

KJ=<J-l>*KOFST 

DO  987  K-  1>  LM 

JK=K+KJ 

JK1=JK+NPB 

JK2=JK+NRN1 

IFCJ.  EQ.  ECL •' JK2 >  =UP < K > 

IFCJ.  GE.  M>  JK2=JK2+K0FST 

987  ECL  <  JK2  >  =ECL  <  JK1 > 

988  CONTINUE 

989  CONTINUE 
NCELL*NCELL+1 
I TERM«5 

GO  TO  9982 

9912  DO  9913  J-l<  LDB 
INDP1=I+J  •‘LMDP 

ieql=trtr:indpi> 

IF<IEQL.  LT.  0) IEQL=-IEQL 
IF<IEQL.  EQ.  NR)GO  TO  991 

I F < TRTR <  I NDP1  ) .  LT.  0.  .  AND.  TA'RC  INDP1).  NE.  -999.  )GO  TO  991 

9913  CONTINUE 
ITERM=7 

9982  CONTINUE 

IFCI.  LT.  NII1N)  NIMN=I 
DO  982  J»l,  NUMP 
Jfl= I +<J-1)*N I DP 

982  TRTR<JR)=0. 

IRCT<I>«0 

991  CONTINUE 
1030  CONTINUE 
~1031  CONTINUE 
C037  FORMAT  <3115) 

1040  DO  1  J-l,  I EM 
JR-J+NEMC 
I START- I Cl < JR) 

IST0P-IC2C JR)  D"30 


DO  2  I-ISTRRT, 1STOP 
MH-W< 1-1) 

IF <W<  I  >.  OT 

IF<W< I+l>.  o i  I+l) 

2  HB<I)-MH 
1  CONTINUE 
RETURN 
END 


!'-  M 
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♦♦♦LISTING  FOR  ERTi  TRACK  FTN 

*N 

SUBROUTINE  TRACK 
IMPLICIT  INTEGERS  <I-N> 

INTEGERS  T, TA 

COMMON  /INTL/  MHSN,  MNSN,  HM,  FNSN 

COMMON  /NVLIS/  NVARM,  NCARMj NVO,  ICO,  10, JO, JVR, KTL 
COMMON  /TLIS/  T,  TA,  JDAV,  JHR,  JMIN,  JSEC,  IDAV,  I HR, MIN,  ISEC 
COMMON  /PNTRS/  NCMX,  NVM I N,  NVMX,  IELSN,  NSCAN,  IESNL. NVSCN, NT 
COMMON  /DATA!/  ECL < 224 >,  NOFST.  KOFST,  I CLAD,  NRN1 
COMMON  /ECONST/  EARTH, VMK 
COMMON  /CNT /  COSPH I , S I NPH I ,  C0SPH2 
COMMON  /A22/  SINA,  COSA,  DELTAS,  ISCANF,  NEL,  RI ,  SA 
COMMON  /KTA/NV,  NC, UCN 
C 

C  INITIALIZE 

C 

UCN-0 

VKM-COSPH I ♦VMK 
VKM2“VKM*VMK 
SAVKM2-SA+VKM2 
SAVKM-SA+VKM 
NAN2-NAN1+1 
DO  10  M-l,  NOFST 
MA- <  NOFST-M ) ♦KOFST+NANZ 
IFCECL<MA>.  GT  0  >  GO  TO  22 
10  CONTINUE 
GO  TO  41 

22  NCMX-NOFST-M+1 
DO  30  M-l,  NCMX 
M1-1+  < M-l > ♦KOFST +NAN1 
M2-1+M1 
M3-1+M2 
M4-1+M3 
M5-1+M4 
M6-1+M5 
M7-1+M6 

C  MRITEC6.  50>M1,  M2,  M3,  M4,  M5,  M6,  COSPHI,  VKM.  ECL<Mi>,  ECL<M2>, 

C  ♦  ECL<M3>,  ECLCM4),  ECLCM5>,  ECL<M6>,  NANI 

C  50  FORMAT < IX-  616,  4X,  3HC0S,  F10  5,  4X,  3HVKM,  F10  2,  /4HECL-, 

C  +  6F8  2,  4X,  4HNAN1,  16  > 

ECL  <  M5 ) -ECL  < M5 > ♦SAVKM 
VKME -VKM/ <  ECL  <  M2  > ♦ECL  <  Ml >  > 

ECL  <  M3  > -ECL <  M3 ) +VKME 

ECL  <  M4 > -ECL  <M4 > ♦VKME 

ECL  < Ml > -ECL  < Ml ) ♦SAVKMZ 

R2-ECL  <  M3 ) ♦ECL  <  M3  > +ECL <  M4  > ♦ECL  <  M4  > 

ECL  <  M?  > -ECL  <  M6  > 

30  ECL  <  M6  > - < SQR T  <  R2 > >  *S I NPH I ♦RZ^COSPHZ 
FNSN-FNSN+1. 

NSCAN-NSCAN+1 
IF<NSCAN  Nf  1>G0  TO  41 
KTL-T 
JDAV- I DRV 
JHR- I HR 
JMIN-MIN 
JSEC- I SEC 
DO  40  NC-1.  NCMX 
NV-NC 

CALL  ATRAK 
40  CALL  BTRAK 


l>  -.W 
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♦♦♦LISTING  FOR  ERT1ATRAK.  FTN 

*N 

SUBROUTINE  ATRAK 
IMPLICIT  INTEGERS  <I-N> 

INTEGERS  T,  TA 
LOGICAL  PR1,  PR2,  PRIN2 
INTEGERS  I VCL < 736 > 

COMMON  /2L00K/  I20FF,  2ARV<91 > 

COMMON  /TLIS/  T, TA,  JDAV,  JHR,  JMIN,  JSEC,  IDAV,  IHR,  MIN,  ISEC 

COMMON  /NVLIS/  NVARM,  NCARM,  NVO,  ICO,  10,  JO,  JVR,  KTL 

COMMON  /PNTRS/  NCMX.  NVMIN,  NVMX,  IELSN,  NSCAN,  IESNL,  NVSCN,  NTT 

COMMON  /INTL.  ‘  MHSN, MNSN, HM,  FNSN 

COMMON  /FLGS/  ITVPE,  PR1,  PR2,  PRIN2 

COMMON  /DVAL.  ’  DELA 

COMMON  /DATA2/  VCL<736>, MXVC, NVCi 

COMMON  /0ATA1/  ECLC224>,  NOFST,  KOFST,  I CLAD, NANI 

COMMON  /KTA/NV,  NC, DELW 

EQU I  VALENCE  <  VCL  M),  I VCL  <  1 )  > 

1F<NC  LE  0  OR.  NC  GT  NCMX>  GO  TO  10® 

I F  <  NV.  LE  0  OR.  NV.  GT.  N VARM  >  GO  TO  100 

NCEC*  <  NC  -1 > ♦KOFST  +NAN1 

NVVC“ <NV  -1 ) ♦NVCl 

NCA-2+NCGC 

I2-ECL<NCA> 

NCA-NCA*! 

X-ECL<NCA> 

NCA-NCA+1 

V-ECL<NCA> 

NCA-6+NCEC 
H*ECL  <  NCA  > 

NVA*9*NVVC 

I F < I VCL < NVA >  GT  ®>  GO  TO  10 

C 

C  DEFINE  LOIIEST  ELEVATION  VALUES 

C 

NVA-l+NVVi ' 

VCL  <  NVA  >  ■' ; 

NVA-NVA+1 

VCL  <  NVA  > ■ V 

NVA-NVA+1 

I  VCL  <  NVA  >  i*  1 2 

NVA-NVA+1 

NCA-NCEC+1 

VCL  C  NVA  > -ECL  <  NCA ) 

NVA-NVA+1 

IVCL<NVA>-T 

NVA-NVA+1 

VCL<NVA>«H 

NCA-3+NCEC 

NCA1-4+NCEC 

R2-ECL «.  NCA  )  ♦ECL  <.  NCA  >  +ECL  <  NCA1  > ♦ECL  <  NCA1  > 

R-SQRT<R2> 

NVA— NVA+i 

VCLCNVAJ-R 

NVA-NVA+i 

~  IF<IVCL<NVA>  .HE.  0)  GO  TO  10 
NTT-NTT+1 
I  VCL  <  NVA  > -NTT 


C 

c 


INCREMENT  ATTRIBUTE  ARRAYS 


U.34 


ono  nno  non 


ltd  i^U=l^+l^U»-»- 

IF<  IZL.  LT.  1)  IZL=1 
IF<  IZL.  GT.  91>  IZL=91 
Z=ZflRV< IZL) 

17  NVfl=9+NWC 

I VCL  <  NVfl  )  =  I VCL  <  NVfl  )  +1 

NCfl=7+NCEC 

NVA=23+NVVC 

VCL  <  NVfl  )  =VCL  <  NVfl  )  +ECL <  NCfl  ) 

NVfl=10+NVVC 

VCL  <  NVfl  )  = VCL  <  NVfl  )  +  Z 

NVfl=NVfl+l 

VCL < NVfl )=VCL< NVfl)  +  Z*X 
NVfl=NVfl+l 

VCL < NVfl )=VCL< NVfl)  +  Z*V 

NVfl=14+NVVC 

HL=VCL<NVfl) 

IF<HL  .  GT.  0.  .  OR.  FNSN  LT.  2.  )  GO  TO  11 

HL=H-DELfl*R 
11  NCfl=NCEC+l 
NVfl=13+NVVC 

VCL  <  NVfl ) =VCL  <  NVfl )  +Z*  <  H-HL ) *ECL  <  NCfl ) 

SUMMIT  VALUES 

NVA=NVfl+l 
VCL < NVfl )=H 
NVfl=NVfl+l 
I VCL<NVfl)=IZ 

■  PEAK,  BASE,  AND  TOP  HRRflVS 

NVfl=NVfl+l 
I ZP= I VCL (NVfl) 

IFCIZP.  GT.  IZ)  GO  TO  99 

SET  PEAK 

32  NVfl=16+NVVC 
I VCL  <  NVfl )  =  I Z 
NVfl=NVfl+l 
VCL < NVfl )=H 

99  IF<.  NOT.  PRIN2)  GO  TO  100 
NCfl=NCEC+l 
NVfl=9+NVVC 
100  RETURN 
END 


D-35 


ooo 


* 
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♦♦♦LISTING  FOR  ERT1 : BTRRK.  FTN 
*N 

SUBROUTINE  BTRRK 
IMPLICIT  INTEGERS  <I-N> 

INTEGERS  T,  TR 
INTEGER  ^4  IVR<192> 

COMMON  /DRTR1/  £CL<224>,  NOFST,  KOFST,  I CLAD.  NRN1 
COMMON  /DRTR3/  VR<192>,  MXVR,  NVR1 

COMMON  /TLIS/  T,  TR,  JDRV,  JHR,  JMIN,  JSEC,  IDRV,  IHR,  MIN,  1SEC 
COMMON  /NVLIS/  NVRRM,  NCRRM,  NVO,  ICO,  10,  JO,  JVR,  KTL 
COMMON  /KTR/  NV,  NC,  UCN 
EQU I VRLENCE  <  VR  < 1 > ,  IVR<1>> 

IF<NV.  LE  0.  OR  NV.  GT.  NVRRM >  GO  TO  10 
IF <NC  LE.  0  OR.  NC.  GT.  NCRRM >  GO  TO  10 

DEFINE  LRST  ELEVRTION  VRLUES 


NCEC=  C  NC-1 > ♦K0FST+NRN1 

NVVR-  <  NV'-l  >  ♦NVRl 

NCR-3+NCEC 

NVR-l+NVVR 

VR  <  NVR  > “ECL  <  NCR  > 

NCR-NCR+1 

NVR=NVfl+l 

VR  <  NVR  > “ECL  <  NCR  > 

NCR-2+NCEC 

NVR-NVR+1 

VR  <  NVR ) *ECL  <  NCR  > 

NCR-NCR-1 

NVR-NVR+1 

VR<NVfl>"ECL (NCR) 

NVft-NVR+1 

IVR<NVfl>-T 

NVR-NVR+1 

NCR-6+NCEC 

VR  <  NVR  >  **ECl .  <  NCR  > 

10  RETURN 
END 


n-j<> 
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♦♦♦LISTING  FOR  ERT 1 : COMPAR.  FTN 

*N 

SUBROUTINE  COMPAR 
IMPLICIT  INTEGERS  CI-N) 

INTEGERS  T,  TR 

LOGICAL  PR1, PR2, PRIN2,  RRAV 

INTEGERS  IVCLC73S>,  IVRC192J 

COMMON  /DATA1/  ECL<224>,  NOFST,  k'OFST,  I  CLAD.  NANI 
COMMON  /DATA2/  VCLC736),  MXVC,  NVC1 
COMMON  /DATA3/  VRC192>,  MXVR, NVR1 

COMMON  /CDRAVS/  ICC32,  10  >,  CC32,  9>,  IOC 32,  10  >,  DC  32,  9>,  IM,  JM 
COMMON  /FLGS/  ITVPE,  PR1,  PR2,  PRIN2 
COMMON  /VPARM/  VX,  VV 

-  COMMON  /PNTRSX  NCMX, NVMIN, NVMX, IELSN, NSCRN, IESNL, NVSCN, NT 

COMMON  /CONST.'  VMISW,  DIV,  VMAG,  VMISWM,  ZDIV,  ADIV,  Al,  A2,  A3,  Bl,  B2,  HDIV 

COMMON  /TL I S/  T,  TA,  JDAV,  JHR,  JMIN,  JSEC,  IDAV,  I  HR,  MIN,  I  SEC 

COMMON  /NVLIS/  NVARM,  NCARM,  NVO,  ICO,  10, JO, JYR, KTL 

COMMON  /CNT/  COSPHI,  SINEL,  C0SPI2 

COMMON  /INTL/  MHSN,  MNSN,  HM,  FNSN 

COMMON  /KTA/  NV,  NC,  UCN 

COMMON  /RSLV/  IUV1C32),  IUV2C32),  IUC1C16>,  IUC2C16>, 

+  UVC32>,  UCC16>,  NCR 

EQU I VALENCE C VCL C 1  > ,  IVCL<1>>,  CVRC1>,  IVR<1>> 

IDV»1 

ICV*1 

—  IFCNCMX.  LE  0)  RETURN 
NVMXP=1 

IFCNVMX.  GT.  1)  NVMXP=NVMX 
DO  3  1*1.  NVMXP 
UVCI>=0. 

IUV1CI>=0 
~  '3  IUV2C I >”0 

DO  4  1=1, NCMX 
UCCI>=0 
IUC1CI>=0 
4  IUC2C I >=0 
DO  5  1=1, 1M 
ICCI, 1>=0 
IDC I, 1>=0 
DO  7  J=l,  JM 
J1«J+1 
ICCI, J1>=0 
IDC I,  J1>=0 
CCI,J>=0. 

DC  I,  J>=0. 

7  CONTINUE 
3  CONTINUE 

DO  10  NC*1,  NCMX 
NCEC«CNC-1>+K0FST+NAN1 

-  NC1=1+NCEC 

NC2-NC1+1 
NC3-NC2+1 
NC4-NC3+1 
|  NC6-NC4+2 

NVC-0 

-  DO  40  NV*1,  NVMXP  - 

NVVC* C NV-1 > *NVC1 
NWR«  CNV-l)  ♦NVRl 
MLAST=0 

OCLWL**  p 

NLR*20+NVVC 


Hircx^yHMyve - 

NLR2=5+NVVR 

IF < I VCL<NLR>  LE  0  .  AND.  IVCL<NLR1>  LE.  0>  GO  TO  40 
DTTR=T-IVR<NLR2> 

ATEST* ( VNAG+DTTA>  *  <  VMAG+DTTA> +VM I SUM 

NRR=1+NVVR 

NCA=21+NVVC 

DELX*ECL  <NC3>-v'R<  NRA  )  -  VCL  <  NCA  >  +DTT  A 

delx2-delx*del:< 

IFCDELX2  .  GT.  ATEST)  GO  TO  20 

NRR=2+NVVP 

NCA=21+NWC 

DEL V«ECL  <  NC4  >  -VR  <  NRA  >  -  VCL  <  NC  A  >  +DTT  A 
DELV2-DELV+DELV 
I F  <  DELV2  .  GT  ATEST>  GO  TO  20 
C 

C  ASSOCIATED,  FIND  BEST 

C 

NRA=3+NVVR 

NRA1=NRA+1 

NRA2=NRAl+2 

DELW*RBS <  ECL  <  NC2  > - VR  <  NRA  >  > ♦ZD I V+  <  DELX2+DELV2 > ♦ 

1  DIV+1.  +  ABS  <  ECL  <  NCI  ">  -  VR  <  NRA1  >  >  ♦AD  I V 

2  +ABSCECL<NC6>-VR<NRA2>>^HDIV 
IF<DELWL  GT.  VMISW)  MLAST-1 

GO  TO  20 

30  NVC=NVC+1 

IFvNVC  GT  1)  GO  TO  31 
IUC1<NC>=NV 
UC<NC)=DELW 
GO  TO  33 

31  IF<IUC1<NC>  .  EQ  NV>  GO  TO  33 
NVT=NV 

IFCDELU!  LT  UC<NC>>  GO  TO  32 

35  I F < I UC2 < NC >  LE  0)  GO  TO  36 
I  =  I UC2  <  NC  > 

IF<  I .  GT.  IM>  GO  TO  361 
39  J=ID< I, 1>+1 
IDCI, 1>=J 

IF<  J  LE  JM>  GO  TO  37 
C  PRINT  101,  JM,  NSCAN,  NVT,  NC,  J 

JO=JO+l 
J=JM 

37  IDCI,  J+D-NVT 
D< I,  J  >=DELW 
GO  TO  33 

36  1*1 DV 
IDV=I+1 

IF< I  .  LE  IM)  GO  TO  38 
C  PRINT  102,  IM,  NSCAN,  NVT,  NC,  I 

361  10*10+1 
I  =  IM 

38  IUC2<NC>«I 
GO  TO  39 

32  DX=UC<NC> 

UL'.NC)*DELl4 
DELW«DX 
NVT* I UC1 <  NC  > 

IUC1<NC)-NV 
GO  TO  35 

33  IF < IUV1<NV>  NE  0>  GO  TO  21 
IUV1<NV>-NC 

UV<NV>-DELW 
GO  TO  40 

C  "  n-3S 

C  CLUSTER 

■  -  — 


nor}  I  non  *,  non  on 


21  IF<  IUV1<NV?  EQ  NC)  GO  TO  40 
NCT-NC 


IFCOELW  LT  UV < NV >  >  GO  TO  22 
25  IF C  IUV2<NV>  EQ.  0>  GO  TO  26 
I*IUV2<NV> 


IF<I  GT  IM>  GO  TO  261 
J«ICO,  1>+1 
1CU.  1>-J 

IF<  J  LE  JM>  GO  TO  27 
PRINT  101,  JM,  NSCAN,  NV,  NCT,  J 

FORMAT ('  NO  OF  CELLS  IN  CLUSTER  EXCEEDS  JM  * 

J0*J0+1 

J*JN 


,  SllO1' 


27  Ji*J+l 


ICO, J1>=NCT 
CO,  J>*DELW 
GO  TO  40 


26  I*ICV 


ICV*I+1 

IF< I  LE  IM>  GO  TO  28 
PRINT  102,  IM,  NSC  AN,  NV,  NCT.  I 

102  FORMAT < '  NO  OF  ENTRIES  IN  CLUSTER  ARRAY  EXCEEDS  IM 
1  5I10> 


261  10*10+1 
I*  IN 

28  IUV2<NV>«I 
GO  TO  29 
22  DX*UV  <  NV  > 

UV  <  NV  >  *OELW 
DELW*DX 
NCT  * I UV1 <  NV  > 
IUVi<NV>=NC 
GO  TO  25 


NO  COMPAR  V,  TRY  VCL 

20  NCA=5+NWC 

DEL  T*  T- 1 VCL  <  NC A  > 

ATEST-VM I SWN+  C  VNAG+DELT  >  +  < VMAG+DELT > 

NCA*1+NVVC 

NCA1-21+NVVC 

DELX*ECL  <  NC3  > -VCL  <  NCA > -VCL  <  NCR1 > *DELT 

DELX2*DELX+DELX 

IF  <DELX2  GT  ATEST>  GO  TO  40 

NCA*2+NVVC 

NCA1*22+NVVC 

OELV-ECL  <  NC4 > -VCL < NCA > -VCL  <  NCA1 >  +DELT 

DELV2-DELV+DELY 

I F  <  DEL V2  .  GT  ATEST>  GO  TO  40 

NCA*3+NWC 

NCA1-NCA+1 

NCA2*NCRl+2 

DELW»A8S  t  < ECL C NC2 > -FLOAT  < I VCL < NCA  >  >  >  +2D I V  > ♦ <  <  DELX2 
'  1  +DELV2>*DI V>+1  +  A8S  <  ECL  <  NCI  >  -  VCL  v  NCA1  >  >  +  AD I V 

2  ♦  ABS  <  ECL  <  NC6  > -VCL  <  NCA2 ) >  *HP 1 V 

IF < ML AST  NE  0  AND  DELWL  LT.  DELW > OELW*PELWL 
I F  <  DEL  W  LE.  VM I SW >  GO  TO  20 
4©  CONTINUE 

IF <NVC  .  GT  0?  GO  TO  10 


ISOLATED  CELL.  NO  COMPAR 


I F < NVMX.  LT.  NVMIN?  GO  TO  501 
00  50  NV*NVM1N,  NVMX 
NVVC*  C  NV-i  >  *NVC1 

J  _ 


non  o o  o 


NLH=^tfflJWC - 

NCAi=9+NS'VC 

IF<I\'CL<NCA>  EQ  0  AND  IVCLvNCAl  >  CO  0>  GO  TO  55 

50  CONTINUE 
501  NV-NVMX+1 

IF<NV  .  LE.  NVARM >  GO  TO  51 
PRINT  103»  NVARM, NV 

103  FORMAT < '  <C>  NO  OF  ACTIVE  TRACKS  EXCEEDS  ARRAV  MAXIMUM  NVMX 
1  2 1 10  > 

NVO=NVO+l 
NV=NVARM 

51  NVMX=NV 
55  NVMIN=NV 

I F  <  NVMX.  LE  0>  NVMX=NV 
UCN=UC<NL  > 

CALL  ATRAK 
CALL  BTRAK 
I UV1 <  NV  >  =  -NC 
IUC1<NC>*-NV 
UV<NV)=0  0 
UC<NC>=0  0 
10  CONTINUE 

HAVE  LIST  OF  COMPARISONS,  NOW  RESOL V  CONFLICTS 

DO  60  NV=i,  NVMXP 
I F < I UV1 < NV >  LT  0?  GO  TO  60 
IFUUVUNV).  EQ  0>  GO  TO  61 
NC=IUV1<NV> 

IF <NC  GT  NCMX >  GO  TO  61 
IF< IUC1<NC>  LE  0)  GO  TO  61 

I F  <  I UV2 < NV >  EQ  0  AND  IUC2CNO  EQ  0>  GO  TO  70 
NCR=NC 
CALL  RESOL V 
NC=NCR 
GO  TO  60 
70  UCN=UCCNC) 

CALL  ATRAK 
CALL  BTRAK 
I U VI <  NV  > = - 1 UV1 <  NV  > 

I UC1 < NC  >  =- 1 UC1 < NC  > 

UVCNV>=0  0 
UC<NC>=0  0 
GO  TO  60 

C  NO  NC  COMPAP,  FIX  HEIGHT  STATISTICS 
61  IFCFNSN  LT  1  1)  GO  TO  60 
NVCR*7+ < NV-1 >  +NVC1 
CA=VCLCNVCA) 

HTC=CR*S INEL+CA+CA*COSP  1 2 
60  CONTINUE 
RETURN 
END 


non 


07/19/79  12:04: 00 

♦♦♦LISTING  FOR  ERT1 : RESOL V.  FTN 

#N 

SUBROUTINE  RESOLV 
LOGICAL  PRIN2,  PR1,  PR2 
IMPLICIT  INTEGERS  CI-N> 

INTEGERS  IVCLC736X  NVT,  IVS,  IVT,  KV,  J,  IDI J,  JJ,  NCT,  ICS,  ICT, 

+KC, ICIJ 

DIMENSION  VC284> 

COMMON  /NVLIS/  NVARM,  NCARM,  NVO,  ICO,  10,  JO,  JVR,  KTL 
COMMON  /COMB/  I VC 22,  7>,  IVMX 

COMMON  /PNTRS/  NCMX,  NVM  IN,  NVMX,  IELSN,  NSCAN,  IESNL, NVSCN, NT 
COMMON  /CDRAVS/  ICC32,  10>,  C<32,  9>,  1DC32,  10>,  DC32,  9>,  IM,  JM 
COMMON  /FLGS/  ITVPE,  PR1,  PR2,  PRIN2 

-  COMMON  /DATA2/  VCLC73€>,  MXVC, NVC1 

COMMON  /CONST/  VMISW,  DIV,  VMAG,  VMISWM,  2D  IV,  AD  IV,  Al,  A2,  A2,  Bl,  B2,  HD  IV 
COMMON  /CNT/  COSPHI,  SINEL,  C0SPI2 
COMMON  /INTL/  MHSN,  MNSN,  HM,  FNSN 
COMMON  /KTA/  NV,  NC,  UCN 

COMMON  /RSLV/  IUV1<32>,  IUV2C32>,  IUC1C16),  IUC2C16>, 

-  +  UVC32),  UCC16X  NCA 

EQU I VALENCE  C  VCL  C 1  > ,  IVCLC1>> 

I VMX=NVARM 


HAVE  CLUSTER,  ORDER  LISTS 

c - 

DO  4  1=1,6 
4  IVC1,  I)=0 
I VT =4 
IVS=5 
ICT =1 

-  ICS=2 

JV=1 

JC=1 

LV=0 

LC=0 

KV=1 

-  KC=1 

NCT =NCA 
-  K0F2=128 
K0F3=256 

IFCNCT.  LE.  0.  OR.  NCT.  GT.  NCMX)  GO  TO  100 
C  PROCESS  NCT 

- 65  IFC  IUC1CNCT>  .  LE.  0.  OR.  IUC1CNCTX  GT.  NVMX>  GO  TO  66 

IFCUCCNCTX  LE.  0.  >G0  TO  66 
NVT  = I UC1 <  NCT  > 

IFCUCCNCTX  GT.  0.  >UCCNCT>=-UCCNCT> 

CALL  COMB INC  NVT ,  I VS,  IVT, KV, J> 

IFC  IUC2CNCT)  .  LE.  0.  OR.  IUC2CNCTX  GT.  IM>  GO  TO  62 

- I_IUC2CNCT5 - 

JX=IDC I,  1> 

IFC  JX  .  LE.  0>  GO  TO  62 
IFC  JX.  GT.  JM>  JX=JM 
IDCI,  1>=-IDCI, 1> 

JX1=JX+1 

- DO  611  J=2,  JX1 - - - 

IDIJ=l£>CI,  J> 

-  CALL  COMBINC IDI J,  IVS,  IVT,  KV, JJ) 

611  CONTINUE 

C — “  "PROCESS  NVT  —  "  -i_4l 

62  IFC IUV1CNVT)  .  LE.  0.  OR.  IUV1CNVTX  GT.  NCMX>  GO  TO  62 


o  o  o 


IF <UV<NVT>.  LE.  0.  >G0  TO  62 
NCT  —  I UV1 <  NVT  > 

I F  <  UV  <  NVT  > .  GT.  0.  )  UV  <  NVT  5  =-U V  <  NVT  ) 

CALL  COMBTN<NCT,  ICS,  ICT, KC, J> 

IF<IUV2<NVT)  .  LE.  0.  OR.  IUV2<NVT>.  GT.  IM>  GO  TO  63 
I=IUV2<NVT> 

JX=IC<I,  1> 

IF<  JX  .  LE.  0>  GO  TO  63 
IF<JX.  GT.  JM>  JX=JM 
IC<I,  1>=-IC<I,  1> 

JX1=JX+1 
DO  621  J=2,  JX1 
ICIJ=IC<I,  J> 

CALL  COMBINC  ICI J,  ICS,  ICT,  KC,  JJ> 

621  CONTINUE 
C 

C  RUN  COMPAR  LIST  TO  FLUSH  OUT  FULL  SET 

62  DO  631  K=JV,  KV 
NVT  =  I V  <  K,  IVS> 

IF  (NVT.  LE.  0.  OR.  NVT.  GT.  NVARM)  GO  TO  631 
IF<UV<NVT>.  LE.  0.  >  GO  TO  631 

I F C I UV1  < NVT ) .  GT.  0.  AND.  IUV1CNVTX  LE.  NCMX>  GO  TO  64 
631  CONTINUE 
GO  TO  66 
64  JV=K 
LC=LC+1 
GO  TO  62 

66  DO  661  K=JC,  KC 
NCT-I V<K,  ICS) 

IFCNCT.  LE.  0.  OR.  NCT.  GT.  NCMX>  GO  TO  661 
IF <UC<NCT).  LE.  0.  )GG  TO  661 

IF<  IUC1CNCT).  GT.  0.  AND.  IUCKNCT).  LE.  NVARM >  GO  TO  6? 
661  CONTINUE 
GO  TO  68 

67  JC=K 
LV=LV+1 
GO  TO  65 

68  IFCLC  .  EQ.  0)  GQ  TO  69 
LC=0 
JV=1 
JC=1 
LV=0 
GO  TO  63 

69  IFCLV  .  EQ.  0)  GO  TO  70 
LV=0 
JV=1 
LC=0 
JC=1 

GO  TO  66 

HAVE  ORDERED  LIST, NOW  FIND  BEST  MATCH 

70  IF<KC.  LE.  1  .  OR.  KV.  LE.  1)  GO  TO  100 
KV=KV-1 

IF<KV.  GT.  IVMX>  GO  TO  100 
KC«KC-1 

IFCKC.  GT.  IVMX)  GO  TO  100 
IMSM=0 

DO  701  K=l,  KV 
NV»IV<K,  IVS) 

IF<NV.  LE.  0.  OR.  NV.  GT.  NVMXJGO  TO  701 
IF <UV<NV).  LT.  0.  )UV < NV > *»-UV < NV > 

701  CONTINUE 

DO  71  K*1,KC 
KA2-K+K0F2 


D-42 


ooo 


Nnj-NTKnj 

VCK>=0. 

VCKfl2>=0. 

VCKfl2>=0. 

IVCJG  ICT)=0 
IVCK,  IVT>=0 
IVCK,  3>=0 
IVCK,  6>=0 
IVCK,  7>=0 
NC»IVCK,  ICS? 

IFCNC,  LE.  0.  OR.  NC.  GT.  NCMX)  GO  TO  71 
IFCUCCNC).  LT.  0.  >UCCNC)=-UCCNC) 

NV«IUC1CNC> 

IFCNV.  LE.  0.  OR.  NV.  GT.  NVMX)  GO  TO  71 
IFCIUV1CNV).  LT.  0>GO  TO  71 
IFCIUV1CNV).  NE.  NDGO  TO  711 
V  C  K  >  =UC  C  NC  > 

UV  C  NV ) =-UV  <  NV  > 

UC  C  NC ) =-UC  C  NC  > 

GO  TO  71 
711  IMSM® IMSM+1 
71  CONTINUE 

I F  C KV.  LE.  CKC-IMSM)  .  OR.  IMSM.  £Q.  0)  GO  TO  75 
FIRST  ROUND  MIN  WEIGHT  SELECTION 

KNV=0 

KNC=0 

DO  72  K=l, KC 

IFCIVCKj  3  >  .  NE.  0)  GO  TO  72 
NC=IVCK, ICS) 

IFCNC.  LE.  0.  OR.  NC.  GT.  NCMX)  GO  TO  72 
IFCUCCNC).  LE.  0.  l)GO  TO  72 
NV= I UC1 <  NC ) 

IFCNV.  LE  0.  OR.  NV.  GT.  NVMX)GO  TO  724 
IFCUVCNV).  GT.  0.  l)GO  TO  725 
IVCK,  7 ) =NV 
KNV=KNV+1 

724  I=IUC2CNC) 

IF  Cl.  LE.  0.  OR.  1.  GT.  1M)  GO  TO  721 
JX=IDCI,  1) 

IFCJX.  LT.  0)JX=-JX 

IFCJX.  LE.  0.  OR.  JX.  GT.  JM)  GO  TO  721 

NV=0 

DWT=999. 

DO  723  J=l,  JX 
J1=J+1 
NVT»=IDC  It  Jl> 

IFCNVT.  LE.  0.  OR.  NVT.  GT.  NVP1X)  GO  TO  723 
IFCUVCNVT)  .  LE.  0.  1)  GO  TO  723 
DELW=DCI,  J) 

IFCDELW.  LE.  0.  1)  GO  TO  723 
IFCDELW.  LT.  DWT ) DWT =DEL  W 
IFCDWT  ,  EQ.  DELW)  NV=NVT 
723  CONTINUE  - 

IFCNV  .  LE.  0  OR.  NV.  GT  NVMX.  OR  DWT.  GT.  VMISW.  OR  DWT.  LT  . 
1  GO  TO  721 
GO  TO  726 

725  DWT-UCCNC) 

726  KR2=K+K0F2 
'  VCKR2)»DWT 

IFCUVCNV).  GT  0  >UVCNV>— UVCNV) 

IFCUCCNC).  GT  0  )UCCNC)»-UCCNC> 

IVCK,  ICT)-NV 
GO  TO  -72 
721  KNC“KNC+1 


n  -  4  3 


n  n  n 


irVMXU.  Ul.  KU  UU  IU  (  d 

IVCKNC,  I  VT)*=K 
72  CONTI HUE 

I F <. KNV.  LE.  0.  AND  KNC  LE  Q)GO  TO  75 
IFCKNC  EQ.  0.  OR  KNC.  GT  KC)  GO  TO  80 

CASCADE  REORDER  OF  COMPAR  LIST 


J«0 

731  J=J+1 

IFCJ  .  GT  KNC)  GO  TO  80 
K=IVCJ, IVT) 

IF <K.  LE.  0  OR  K  GT  KC)  GO  TO  731 
NC*IV<K,  ICS) 

IFvNC  LE.  0.  OR.  NC  GT  NCMX)  GO  TO  739 
NV=IUC1<NC) 

IFCNV  LE  0  OR.  NV  GT  NVMX)  GO  TO  739 
DO  738  L=l>  KC 

IF<IVvL,  3).  EQ  NV ) GO  TO  73S1 
IFCIVCL,  ICT)  EQ.  NV)GO  TO  7382 
738  CONTINUE 
GO  TO  739 
7381  NCT«IV<L,  ICS) 

DELT=VMISW+V<L> 

GO  TO  7383 
73S2  NCT  = I V<L<  ICS) 

KA2=L+K0F2 
DELT =VM I SW  * VCKA2 ) 

7383  KT=L 

IF <NCT.  LE  0  OR  NCT  GT  NCMX)GO  TO  739 
IFvIUCl<NCT)  LE  0  OR  IUC2CNCT).  LE.  0>PO  TO  739 
DELM=UC  < NC ) 

IFvDELW  LT  0  )DELW—DELW 
I=IUC2<NCT ) 

IFU.  LT  0)1=-I 

IFU  .  LE  0  OR  I  GT  IM)  GO  TO  739 
JX=ID< 1/ 1) 

IF < JX.  LT  0)JX=-JX 

IFCJX  LE  0  OR  JX.  GT  JM)  GO  TO  739 
DWT“999. 

NVB-0 

DWT1=999 

NV1=0 

DO  732  L=l> JX 

IF<D<I,L).  LE  0  1)  GO  TO  732 
IF<D<I,L)  LT  DWT)DWT-D< I> L) 

Ll=L+i 

NVT1  =  1ABS<  IDO, Ll>) 

I F <. DWT  EQ  Dv  I  <  L  >  >  NVB-NVTi 

IFCNVTl  LE  0  OR.  NVT1.  GT.  NVMX)  GO  TO  732 

IF <UVCNVT1)  GT  0  1  AND  D<I,L>.  LT  DWT1)DWT1-D< I,  L) 

IF<DWTi  EQ  D<I,L>)  NV1-NVT1 

732  CONTINUE 

IFCNVB  LE  0  OR.  NVB.  GT  NVMX  OR  DWT  LE  .  1  OR  DWT  GT  VMISW) 

1  GO  TO  739 

734  IFvNVl  LE  0  OR  NV1.  GT  NVMX  OR  DWT1  GT.  VMISW  OR  DWT1.  LE.  .  1) 
1  GO  TO  735 
DELW1-DELW+DWT1 
DELW2-DELW+DWT 
IFCDELWl  C.T  DELT)  GO  TO  735 
IV<K,  6>-N'‘ 

KA3-K+K0F  ; 

V<KA3)-DELW 
I  V<KT<  6>-NVTl 

KA3-KT+KA3  D-44 

V<KA3)-DWT1 

i 


non 


I F  C  DELU12.  GE.  DELHI)  GO  TO  729 
I VCKT, 7 >  =NVB 
GO  TO  739 

735  DO  736  1=1,  KC 

IF  C  I  VC  I,  3>.  EQ.  NVB)  GO  TO  739 
IFCIVCI,  ICT>.  EQ.  NVB)  GO  TO  737 

736  CONTINUE 
GO  TO  739 

737  KR2=I+K0F2 

DELT =DELT+ VC  Kfl2 ) 

DELW2=DELW+DWT+ VM I SW 

IFCDELW2.  GT.  DELT)  GO  TO  739 

IVCK,  6>=NV 

KA3=K+K0F3 

VCKA3)=DELW 

IVCKT,  6>=NVB 

KA3=KT+K0F3 

KA2=I+K0F2 

VCKA3)=DWT 

IVCI,  ICT>=0 

VCKA2>=0. 

739  IVCK, 7>=0 
KNV=KNV-1 
GO  TO  731 

EXCHANGE  PAIRS  FOR  MIN  MEASURE 

80  IFCKNV.  LE.  0.  OR.  KNV.  GT.  KC>  GO  TO  75 
DO  801  K=l,  KC 
NVB=I VCK, 7) 

IFCNVB.  LE.  0.  OR.  NVB.  GT.  NVMX>  GO  TO  801 
IFC IUV2CNVB).  LE.  0>  GO  TO  801 
NC=I VCK,  I CS> 

IF<NC.  LE.  0.  OR.  NC.  GT.  NCMX)  GO  TO  801 
NV=  IVCK,  2) 

DO  802  L«l, KC 

IFCNVB.  EQ.  IVCL,  ICT).  OR.  NVB  EQ  IVCL,  3>)  GO  TO  803 

802  CONTINUE 
GO  TO  801 

803  NCB=IVCL,  ICS) 

I  =  I UV2  C  NVB ) 

IF  Cl.  LT.  0>I=-I 
JX=ICCI, 1) 

IFCJX.  LE.  0.  OR.  JX.  GT.  JM>  GO  TO  801 
DO  807  J=l,  JX 
J1=J+1 

IFCICCI,  Jl).  EQ.  NCB)  GO  TO  808 

807  CONTINUE 
GO  TO  801 

808  DSET=C< I,  J) 

DUC=UC  C  NC  > 

IFCDUC.  LT.  0.  >DUC=-DUC 
DELWB=DSET  +DUC 
KR2*K+K0F2 
KR3=K+KOF3 
DELW1=VCKA3> 

KR=KR-KOFST 

IFCDELW1.  LE.  .  1)  DELW1=VCKR2) 

IFCDELW1.  LE.  .  1)  DELW1=VCK> 

KR2-L+K0F2 

-  KR3=L+K0F3  "  - 

DELW2»VCKR3> 

IFCDELW2.  LE.  .  1>  DELW2-VCKR2) 

IFCDELM2.  LE ..  1)  DELW2-VCL) 

DEL W=DELW1 +DELM2 
TF<T>F»  U  I  F  OFI  UR)  00  TO  QO-i 


P-45 


o  on 


V<KA3>-UC<NC> 

1F<V<KA3>  LT  0.  >V<KR3>=-V<KA3> 

KA3=L+KQF3 
V<KA3)**DSET 
I  V<K, 6>*NVB 
IV<L,  6>«NV 
801  IV<K, 7>»0 
75  CONTINUE 

UPDATE  ATTRIBUTES 

DO  78  K*l,  KC 
KA3=K+K0F3 
KA2-K+K0F2 
I F  <  PR I N2 ) 

1  WRITE <6,  788)  IV<K.  ICS),  IV<K,  6),  IV<K,  ICT), 

2  IV<K,  3),  V<KA3>,  V<KA2>,  VCK) 

788  FORMAT  <  IX,  415.  3F8  2) 

NC=IV<K, ICS> 

IFCNC  LE  0  OR.  NC.  OT.  NCMX>  GO  TO  ?S 
IF<IUC1<NC>  LE.  0)  GO  TO  78 
NV-IVCK, 6> 

IFCNV.  LE  0  OR  NV  GT  NVMX>  GO  TO  810 
DWT-VCKA3) 

GO  TO  820 

810  NV-IVCK, ICT> 

IF<NV.  LE  0  OR  NV  GT  NVMX)  GO  TO  811 
DWT**V<KA2) 

GO  TO  820 

811  NV»IV<K, 3) 

I F  <  NV.  LE  0  OR  NV.  GT  NVMX  >  GO  TO  79 
DWT-VCK) 

820  IFCDWT  LE  0  1  OR  DUT.  GT.  VMISW>  GO  TO  79 
IF< IUV1<NV>.  EQ  0)  GO  TO  79 
UCN-DMT 

call  atra; 

CALL  BTRA.: 

GO  TO  77 

79  I F < NVMX.  LT  NVMIN>  GO  TO  7911 
DO  791  I-NVMIN,  NVMX 
NVVC*  < I-1)*NVC1 
IVA«9+NVVC 
IVA1-IVA+11 

I F  <  I VCL <  I VA1 ) .  EQ  0  AND.  IVCLCIVA).  EQ  0)  GO  TO  792 

791  CONTINUE 
7911  I-NVMX+1 

IF<  I.  LE.  NV  ARM  >  GO  TO  7921 
C  PRINT  103,  NVARM,  I 

C  103  FORMAT < y  <R)  NO  OF  ACTIVE  TRACKS  EXCEEDS  ARRAV  MAXIMUM  NVMX 
C  1  2110) 

NV0“NV0+1 
I “NVARM 
GO  TO  77 
7921  NVMX* I 

792  NV-I 
NVMIN-I 
UCN-0 
CALL  ATRAK 
CALL  BTRAK 

77  IUVi<NV>—  NC 
IUCKNC)—  NV 
UV<NV)-0 
UC<NC>-0 
7»  CONTINUE 
00  99  K-l, KV 


n-4t> 


NV»1V<.K,  lVi»> 

IF <NV.  LE  0.  OR.  NV.  GT.  NVMX>  GO  TO  99 
IF<IUV1<NV>.  LE.  0)  GO  TO  99 
IUV1<NV>  —  IUV1CNV) 

IFCFNSN.  LT.  1.  1>  GO  TO  99 
NVfl«?+  <  NV-1 >  *NVC1 
Vft«VCL<NVfl> 

HTC«Vfl*S I NEL+Vfl*Vfl*COSP I 2 
99  CONTINUE 
100  RETURN 
END 


o  o  o  o 


r  7 

C 

* 
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♦♦♦LISTING  FOR  ERT1 : COMB IN.  FTN 

SUBROUTINE  COMBINCN,  IS,  IT,  K,  J> 

IMPLICIT  INTEGERS  <I-N) 

INTEGERS  N, IS,  IT,  K,  J 
COMMON  /'COMB/  IV<32,  7),  IVMX 

INSERT  N  INTO  ORDERED  RRRRV  IV<K, IT> 

RETURN  NEW  RRRRV  RS  IVCK, IS) 

I  =  IS 
IS=IT 
I  T=  I 
L=0 

DO  10  J=1  K 
IVJ=IV<J,  IT) 

IF< I VJ.  LT  0)IVJ=-IVJ 
IF<  I VJ-N)  20,20,40 
IFCIVJ.  EQ  0)  GO  TO  40 
IV<J, IS)=IV<J,  IT) 

J=K 
L=1 

IV<J,  IS)=N 
DO  50  I=sJ,  K 
IL=I+L 

IV<IL,  IS)=IV<I,  IT) 

K=K+L 

IF<K.  GE.  IVMX)  GO  TO  70 
IVCK,  IS)=0 
GO  TO  80 

WRITE <6, 100)  K,  IVMX 

FORMRT  <  '  ERROR  IN  COMB  IN,  13,  K,  IVMX',  2X,  2110) 
K=I VMX-1 
RETURN 
END 


» 


I 


K 


B 


I 

c 


I 


T 


o  o 
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♦♦♦LISTING  FOR  ERT1 : STRflK  FTN 

*N 

SUBROUTINE  STRflK 
IMPLICIT  INTEGERS  <I-N> 

INTEGERS  T,  TM,  Tfl 

INTEGERS  IVCLC736),  ITCLC21),  IVRC192) 

DIMENSION  flTST <62>, NUMC62) 

LOGICAL  PR1,  PR2,  PRIN2 
DIMENSION  TCL<21>,  DUM<6> 

COMMON  /TMflX/  TM 

COMMON  /FLGS/  I  TYPE,  PR1,  PR2,  PRIN2 

COMMON  /PNTRS/  NCMX,  NVMIN,  NVMX,  IELSN,  NSCAN,  IESNL,  NVSCN,  NT 
COMMON  /TLIS/  T,  Tfl,  JDAV,  JHR,  JMIN,  JSEC,  IDAV,  I  HR,  MIN,  I  SEC 
COMMON  /NVLIS/  NVARM,  NCARM,  NVO,  ICO,  10,  JO,  JVR,  KTL 

COMMON  /CONST/  VMISW,  DIV,  VMflG,  VMISWM,  ZD  IV,  flDIV,  fll,  A2,  A3,  Bl,  B2,  HD  IV 
COMMON  /INTL/  MHSN, MNSN,  HM, FNSN 
COMMON  /VPflRM/  VX,  VV 

COMMON  /DflTfll/  ECL<224>,  NOFST, KOFST,  I CLAD,  NflNl 
COMMON  /DATA2/  VCL<736>,  MXVC, NVCi 
COMMON  /DATA3/  VR<192>,  MXVR,  NVR1 
EQU I VALENCE <  TCL  < 1 > ,  ITCL<1)> 

EQU I VALENCE < VCL <  1 ) ,  IVCL<1>>,  <VR<1>.  IVR<1>> 

DATA  I ZERO/0/,  I OUT/31/,  LAN/62/ 

MAN-0 

NVOUT =NVMX-1 
DO  19  N»l,  LAN 
ATST  <N)°*0. 

19  NUM<N>-0 


VXC=0 

-  VVC=0 

NSN-0 

NACT=0 

NFN=FNSN-1. 

I F < NFN.  LE.  0>  GO  TO  35 

C  WR I TE  <  2  >  KTL,  NVSCN,  JDAV,  JHR,  JM I N,  JSEC 

C  1  ,  NFN,  NSCAN,  NT,  NVMX,  JVR,  NVO 

C  2  ,  10,  JO,  I  ZERO,  VX,  VV,  DUM 

IF<PR1.  OR.  PR2) 

1WRITEC6, 1000 >  JDAV, JHR, JMIN, JSEC,  NVSCN 
1000  FORMAT < 1H1,  3X,  'PROGRAM  TRACK4  <781207)  ', 5X,  '  TRACK  TIME', 414, 
15X,  'VOL  SCAN',  16//,  IX,  'TRK  LOCATION  DBZ  CELL  PEAK  ', 

2'  SURFACE  SUMMIT  VELOCITY  A'/,  IX 

3'  ID  EAST  NORTH  AV  VOL  DBZ  HT  DBZ  HT  DBZ  HT', 

4'  EAST  NORTH  L'/,  IX, 

3'  KM  KM  KM3  KM  KM  KM', 

S'  M/S  M/S  L'> 

DO  20  NV-1,  NVOUT 
~  NVA*< NV-1 > ♦NVCI 

NV23-23+NVA 
A-VCLCNV23) 

DO  21  N-l,  I OUT 
NN-N+MAN 

IF<A.  GT.  ATST<NN>>  GO  TO  22 

- "21  CONTINUE 

NN-NN+1 
22  LAN-MAN 

MAN-MAN+IOUT 
IFCMAN.  GT  I OUT)  MAN*® 

IFML-0 


n-i  9 


>/u  *.->  l  l-l<  i  UU  t 

ML-M+LflN 

IF <ML.  NE  NN  OR  IFML  EQ.  1>  GO  TO  24 

8TST<MM>-fl 

NUM<MM>*NV 

LRN»LftN-l 

IFML-1 

GO  TO  23 

24  flTST <MM)»flTST <ML> 

NUM<MM>*NUM<ML> 

23  CONTINUE 
20  CONTINUE 
50  DO  100  N*l.  I OUT 
NN=N+MfiN 
NV=NUMCNN> 

I F < N V.  EQ.  0>  GO  TO  100 

NVft*  <  NV— 1 >  *NVC1 

NVl*l+NVfl 

NV2=NV1+1 

NV3=NV2+1 

NV4«NV3+1 

NV5*NV4+1 

NV6*NV5+1 

NV7=NV6+1 

NV8=NV7+1 

NV9=NV8+1 

NV10=NV9+1 

NV11=NV10+1 

NV12=NV11+1 

NV13=NV12+1 

NV14=NV13+1 

NV15=NV14+1 

NV16=NV15+1 

NV17=*NV16+1 

NV1S=NV17+1 

NV19=NV18+1 

NV20=NV19+1 

NV21=NV20+1 

N V22=N V21 + 1 

NV23=NV22+1 

NR1=1+  <  NV— 1 >  *NVR1 

NR2=NR1+1 

NR3=NR2+1 

NR4=NR3+1 

NR5*NR4+1 

NR6=NR5+1 

I F < I VCL C N V9 >  .  LE.  0>  GO  TO  200 
22DIV-1.  0/VCL<NV10> 

ITCL<1>-IVCL<NV5> 

TCL < 2 > «VCL < NV11 > *2ZD I V 
TCL < 3 > *VCL C NV12 > *22D I V 

TCL<4>*10.  *RL0G10 < VCL < NV10 > /FLOAT < IVCL<NV9>>> 

TCL < 5 ) -VCL < NV13 > *22D I V 

ITCL<7>*IVCL<NV16> 

TCL<8>»VCL<NV1?> 

ITCL<18>»IVCL<NV3> 

TCL<11>*VCL<NV4> 

TCL<12>-VCL<NV€> 

I M3* I VCL  <  NV16 ) -3 
30  ITCL<13>-IVCL<NV15> 

TCL < 16 > *VCL < NV14 > 

1VL8*IVCL<NV8) 

IF<  IVL8.  LT.  8)IVL8—  IVL8 
ITCL<17>*IVL8  D-SO 

ITCL<18)-0 


ITCL<20>=0 
ITCL<21>=0 
VXT = VCL  <  NV21 > 

VVT =VCL  <  NV22  > 

ITIM=IVCL<NV20> 

JTIM=ITCL<1> 

IF<  ITIM.  EQ.  0  .  OR.  ITIM.  EQ.  JTIM>  GO  TO  40 
DELTM= JT I M- I T I M 
DELTM=1.  0/DELTM 

TCL <  20  >  =  <  TCL C  2  > - VCL  <  NV18  > ) *DELTM 
TCL  <  21 >  = < TCL  <  3  > -VCL  <  NV19  > ) *DELTM 
VXT=TCL<20> 

VVT=TCL<21> 

40  I  VCL  <  N  V20  )  =  I T  CL  <  1 ) 

VCL  <  NV19  )  =TCL  <  3  > 

VCL  <  NV18  >  =TCL  <  2  > 

VCL  <  NV21  >  =fll*VXT*R2*VCL  <  NV21  > +R3*VX 
VCL  <  NV22  >  =A1*VVT+R2*VCL  <  NV22 > +R3*VV 
VRCNR1>=TCLC2> 

VR<NR2>=TCL<3> 

VRCNR3>=TCL<4> 

VR  <  NR4  >  =VCL  <  NV4 ) 

I VR  <  NR5  >  =  I TCL  <  1  > 

VR<NR6>=TCL<8> 

IZVRL=VR<NR3> 

IDTC=FLORT <  IVCLCNV9)  )/<FNSN-l.  >*10.  +.  5 

NSN=NSN+1 

VXC=VXC+TCL<20> 

VVC=VYC+TCL  <  21 > 

VXP=VCL  <  NV21 >  *1000. 

VVP=VCL  <  NV22  >  *1000. 

NRCT =N  ACT +1 

C  WR I TE  <  2  >  TCL.  VCL  < NV39 > .  VCL  < NV40 > .  IDTC 

IF<PR1>  GO  TO  59 
IFC.  NOT.  PR2>  GO  TO  58 

IFC  <  IDTC.  GT.  MNSN.  OR.  IZVRL.  GT.  35>.  AND.  IVCL<NV8>.  GT.  0> 

1  I VCL  <  NV8  >  =— I VCL  <  NVS  > 

IFCIVCLCNVS).  GE.  0>  GO  TO  58 

59  WRITE <6.  1005 >  ITCLC17).  VRCNR1).  VRCNR2).  IZVRL.  TCL<5>.  ITCLC7). 
1TCLC8>.  ITCL<10>. TCL<12>.  ITCL<15>.  TCL<16>. 

2VXP. VVP. IDTC 

1005  FORMAT < IX.  14.  2F5.  0.  13.  F6.  1.  13.  F5  1.  15.  F5.  1.  14. 

+  3F5.  1.  2X.  I2> 

58  DO  41  1=9.  14 
I R= I +NVA 

41  IVCL< IR>=0 

I VCL  <  NV23  >  =0 
GO  TO  100 

"  200  IF < IVCL<NV20>  .  LE.  0>  GO  TO  102 

IF < <T-IVCL<NV20>  >  LE.  TM>  GO  TO  100 
DO  101  1=1.20 
IA=I+NVfl 

101  IVCL< IR>=0 
IA=23+NVA 
VCL<  lfl>=0. 

I F  <  NV.  LT.  NVM  IN)  NVM I  N=N V 

102  IVCL<NV23>=0 
100  CONTINUE 

IF<NSN  .  EQ.  0>  GO  TO  56 
VN=NSN 

-  VX=B1*VXC/VN+B2*VX 

VV=B1*WC/VN+B2*VV 
DO  43  JI=1.  I OUT 
II-IJ+MAN 

- I*NUPKID 

IFU.  EQ.  0)  GO  TO  43 


D-51 


MV11-V  4 

IF<IVCL<NVA>  GT  0>  GO  TO  43 
NA=21-KI-1>*NVC1 
VCL<NA>-VX 
NA=NA+1 
VCL<NA>=VV 
43  CONTINUE 
VXP=VX*1000 
VVP=W*1000. 

56  IF <NVSCN.  GT.  1>  GO  TO  55 
WRITE<6,  1003  > 

1003  FORMAT  </7VlX,  "  VOL  START  TIME  NO  EL  LAST  TRACK', 

1  '  OVERFLOW  AV  VELOCITY  V,  IX, 

2  SCAN  DAV  HHMM  SS  SCANS  SCAN  TOTAL  ACT  ', 

3  '  NV  IC  I  J  EAST  NORTH '/> 

55  WRITE <6,  1004  >  NVSCN,  JDAV,  JHR,  JMIN,  JSEC,  NFN,  NSCAN,  NT,  NACT 
1,  NVO,  ICO,  10,  JO,  VXP,  VVP 

1004  FORMAT  <  IX,  15,  14,  13,  12,  13,  15,  216,  15,  17,  214,  IX,  2F6  i> 
NVSCN=N VSCN ♦ 1 

NVO=0 

ICO=0 

10=0 

JO=0 

JDAV=IDAV 
JHR=1HR 
JMIN=MIN 
JSEC= I SEC 
FNSN=1.  009 
KTL=T 
RETURN 
END 


•\m\’Pi\  i 
S'VMt'l output 


SAMPLE  OUTPUT  FI  Li! 


Two  types  of  output  are  produced  following  each  volume  scan: 

1.  A  list  of  attributes  describing  each  of  the  31*  most  significant 
cells  detected  and  tracked  within  the  volume  scan 

a.  ITCL  (17)  -  the  ID  number  of  this  cell  track 

b.  VR  (1)  -  cell  location  to  east  {(-)west]  of  radar  (km) 

c.  VR  (2)  -  cell  location  to  north  f(-)south]  of  radar  (km) 

d.  IZVAL  =  VR  (3)  -  average  ref lectivity  of  cell  (dBZ) 

e.  TCL  (5)  -  cell  volume  (km^) 

f.  ITCL  (7)  -  peak  reflectivity  of  cell  (dBZ) 

g.  TCL  (S)  -  height  of  peak  reflectivity  (km) 

h.  ITCL  (10)  -  reflectivity  at  base  of  cell  (dBZ) 

i.  TCL  (12)  -  height  of  cell  base  (km) 

j.  ITCI.  (15)  -  reflectivity  at  summit  of  cell  (dBZ) 

k.  TCI  (16)  -  height  of  cell  summit  (km) 

l.  VXP  =  VCI,  (391*10;)  -  cell  velocity  toward  east  [(-)westj  (m/s) 

m.  VYI  =  VCL  (40)*103  -  cell  velocity  toward  north  [(-)south]  (m/s) 

n.  ID'  C  =  I VCL  (9)/(FNSN-l)*10 

number  of  scans  cell  detected/number  of  scans  per  volume  scan 
percent  of  elevation  scans  cell  detected  at 

2.  A  suim try  list  of  statistics  for  the  entire  volume  scan 

a.  NVSCN  -  volume  scan  number 

b.  JDAY,  JHR,  JMIN,  JSEC  -  start  time  of  volume  scan 

c.  NF.\  =  FNSN-1  -  number  of  elevation  scans  processed  in  this  volume  scan 

d.  NSCAN  -  number  of  last  elevation  scan  in  the  volume  scan 

e.  NT  -  number  of  cell  tracks  updated  this  scan 

f.  NACT  -  number  of  possible  cell  tracks  stored  from  current  and 

previous  scans 

g.  NVO  -  number  of  significant  cells  detected  but  not  outputted 

h.  ICO  -  number  of  internally  paired  and  clustered  cells  over  the 

dimensions  of  CI.UST  array  in  subroutine  RF.SOLV 

i.  10  -  number  of  detected  cells  over  array  dimensions  of  cluster  array 

n  subroutine  COMPAR 

j.  .10  -  number  of  detected  cells  over  array  dimensions  of  IC  array  in 

subroutine  COMPAR 

k.  VXP  -  an  estimate  of  the  average  velocity  east  [(-)west]  any  new 

cells  will  have  (m/s),  set  to  0  on  first  scan 

l.  VYP  -  an  estimate  of  the  average  velocity  north  [(-)south]  any  new 

cells  will  have  (m/s),  set  to  0  on  first  scan 


*see  Appenuix  F  -  Option  to  Increase  Number  of  Significant  Cells 


F.-2 


SAMPLE  OUTPUT  OF  FIRST  VOLUME  SCAN 


PROGRAM  TRACK 4  <7S1207>  TRACK  TIME  120  16  46  59  VOL  SCAN 


TRK  LOCATION 

DB2 

CELL 

PEAK 

SURFACE 

SUMMIT  1 

VELOCITY 

A 

IOEAST 

NORTH 

AV 

VOL  DBZ 

HT 

DB2 

HT 

DBZ 

HT  EAST  NORTH 

L 

KM 

KM 

KM3 

KM 

KM 

KM 

M/S 

M/S 

L 

1?  -90. 

40. 

45 

1.  5 

45 

4.  3 

45 

4.  3 

45 

4.  2 

0.  0 

0.  0 

1  -60. 

-63. 

34 

2.  7 

34 

1.  8 

34 

1.  8 

34 

1.  8 

0.  0 

0.  0 

3 

2  -96. 

37. 

41 

4.  0 

41 

2.  3 

41 

2.  3 

41 

2.  2 

0.  0 

0.  0 

3 

18  -89. 

32. 

47 

2.  7 

47 

4.  1 

47 

4.  1 

47 

4.  1 

0.  0 

0  0 

3 

3-56 

—67. 

42 

2  8 

42 

'6.  7 

42 

1.  8 

42 

6  7 

0.  0 

0.  0 

■“* 

T 

19  -80. 

38. 

36 

5.  4 

36 

6.  8 

36 

2.  7 

36 

6  8 

0.  0 

0  0 

7 

20  -86. 

30. 

37 

2.  6 

38 

3.  9 

38 

3.  9 

36 

7.  0 

0.  0 

0.  0 

7 

5  -94. 

34. 

46 

2.  8 

47 

2.  1 

47 

2.  1 

46 

4.  4 

0.  0 

0  0 

7 

-  21  -97. 

53. 

39 

1.  8 

39 

4.  8 

39 

4.  8 

39 

4.  S 

0.  0 

0.  0 

2 

4  -50. 

-66. 

36 

2.  4 

36 

1.  7 

36 

1.  7 

36 

1.  7 

0.  0 

0.  0 

3 

- 31=1817 

- 35 

'41 

177 ' 

41 

—8.  2 

—  41 

2.2 

“41 

-8.2 

"0."0 

0.  0 

3 

22  -91. 

50. 

42 

6.  3 

43 

4.  6 

43 

4.  6 

42 

8  0 

0.  0 

0.  0 

7 

30  -93. 

37 

46 

1.  5 

46 

4.  3 

46 

4.  3 

46 

4.  3 

0.  0 

0.  0 

3 

6  -43. 

63. 

37 

2.  0 

37 

1.  5 

37 

1.  5 

37 

1.  5 

0.  0 

0.  0 

3 

“  23  -37. 

54. 

52 

2.  5 

52 

2.  7 

52 

2.  7 

52 

2.  7 

0.  0 

0.  0 

3 

7  -4. 

25. 

IBIQ 

0.  2 

50 

0.  4 

50 

0.  4 

50 

0.  4 

0.  0 

0.  0 

3 

374“= 727 

53.' 

44 

370 

"45' 

3.  9 

- 

T9 

'“457 

“379 

“  “0.70 

"  070“ 

7 

8  -50. 

60 

41 

4.  3 

41 

1.  6 

41 

1.  6 

41 

1.  6 

0.  0 

0.  0 

3 

24  -75. 

83. 

36 

3.  7 

36 

4.  9 

36 

4.  9 

36 

4.  9 

0.  0 

0.  0 

3 

25  -73. 

56. 

39 

6.  3 

40 

10.  0 

38 

4.  0 

40 

10.  0 

0.  0 

0.  0 

7 

9  -27. 

2. 

43 

0.  4 

44 

0.  5 

44 

0.  5 

44 

0  5 

0.  0 

0  0 

3 

10  -14. 

25. 

41 

0.  2 

41 

0.  5 

41 

0.  5 

41 

0.  5 

0.  0 

0.  0 

3 

ll  -6. 

15. 

42  071 

42 

0.  2 

42 

“0.3 

42 

0.3 

0.  0 

0  0 

3 

12  -80. 

61. 

40 

3.  8 

40 

2.  2 

40 

2.  2 

40 

2.  2 

0.  0 

0.  0 

3 

13  -17. 

10. 

39 

0.2 

39 

8.  3 

39 

0.3 

29 

0  3 

0.  0 

0.  0 

3 

15  -4. 

53 

49 

0.  9 

49 

1.  0 

49 

1.  0 

49 

1.  0 

0.  0 

0.  0 

3 

26  -57. 

82. 

43 

1.  5 

44 

4.  3 

44 

4.  3 

44 

4.  3 

0.  0 

0.  0 

3 

16  -6. 

1. 

41 

0  1 

41 

0.  1 

41 

0.  1 

41 

0.  1 

0.  0 

0.  0 

3 

“27  “-547 

80. 

'42 

17  4 

42 

4r2 

42“ 

472 

“42 

4.2 

0  0 

0.0 

3 

28  -62 

55. 

36 

1.  0 

36 

3.  5 

36 

3  5 

36 

3.  5 

0.  0 

0.  0 

3 

29  -52. 

-70 

45 

1.  1 

45 

3.  7 

45 

3.  7 

45 

3.  7 

0.  0 

0.  0 

3 

~ VOL  5TRRT  TIME  TTO  -EL — LAST - TRACK 

SCAN  DRV  HHMM  SS  SCANS  SCAN  TOTAL  ACT. 
~  1  120' 1646 ”59  '  3  4  32  31 


- OVERFLOW 

NV  IC  I 
25  8  0 


- RV  VELOCITY 

J  EAST  NORTH 
0  0  0  0.0 


APPENDIX  F 

Options  to  Increase: 

Number  of  Significant  Cells 
Number  of  Active  Tracks 


Option:  Increase  Number  ot'  Significant  ('ells  Processed  bach  Scan 


Both  the  real-time  program  and  the  post -mission  program  are  set  up  to 
process  lb  significant  cells  each  scan.  To  increase  the  number  of  cells 
the  dimensions  of  several  arrays  and  address  offsets  must  be  increased. 

In  block  common  I>ATA1  -  for  each  addtional  cell 

increase  dimension  of  IT'l.  by  14 
increa  ;e  NOFST*  by  1 
increase  10I.A0*  by  7 

In  block  common  RSI.V  -  for  each  additional  cell 
increase  dimension  of:  IIIC1,  HJC28UC  by  1 

In  block  common  NV1.IS  -  for  each  additional  cell 
increase  NCARM*  by  1 

In  subroutine  R1-S0I.V  -  for  each  additional  cell 

increase  dimension  of  V  hv  24 
increase  KOT:2tiKOI:3  by  8 

Note:  to  process  more  than  32  cells  each  scan  the  number  of  active  tracks 
must  be  increased  to  at  least  that  number  (see  option  to  increase 
active  tracks'). 


•variable  set  in  block  DATA 


Option:  Increase  Number  of  Active  Tracks  Updated  l.ach  Scan  and  Outputted 

Hath  Voliune  Scan 

Both  the  real-time  program  and  the  post-mission  program  are  set  up  to 
update  and  output  31  active  tracks  at  any  one  scan.  To  inert  ise  the 
number  of  tracks,  the  dimensions  of  several  arrays  and  add res  .  offsets 
must  be  increased. 

In  block  common  DATA 2  -  for  each  additional  track 

increase  dimension  of  VCI,  and  IVOR  bv  23 
increase  MX VC*  by  23 

In  block  common  DATA3  -  for  each  additional  track 

increase  dimension  of  VR  and  IVR  by  u 
increase  MXVR*  by  7 

In  block  common  RS1.V  -  for  each  additional  track 
increase  dimension  of:  HIV  1 ,  IllV2tillY  by  1 

In  block  common  CPRAYS  -  for  each  additional  track 

increase  the  first  dimension  of:  IC.C.lIUiU  by  1 
increase  IM*  by  1 

In  block  common  NV1TS  -  for  each  additional  track 
increase  NVARM*  by  l 

In  block  common  COMB  -  for  each  additional  track 
increase  the  first  dimension  of  IV  by  1 

In  subroutine  STRAK  -  for  each  additional  track 

increase  dimension  of  ATST  t,  NUM  bv  7 
increase  I OUT  by  1 
increase  (.AN  by  2 


‘variable  set  in  block  PAT A 


